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PREFACE 


The  development  of  experimental  facilities  for  rail  vehicle 
testing  at  the  DOT  High  Speed  Ground  Test  Center  is  being  com- 
plemented by  analytic  studies  being  conducted  by  Transportation 
Systems  Center  under  the  UM204  UMTA  Rail  Supporting  Technology 
Program  for  the  Urban  Mass  Transportation  Administration,  Office 
of  Research,  Development,  & Demonstrations.  The  purpose  of 
this  effort  has  been  to  gain  insight  into  the  dynamics  of  rail 
vehicles  in  order  to  guide  the  development  of  the  wheel/rail 
simulators  and  to  establish  an  analytic  framework  for  the  design 
and  interpretation  of  tests  to  be  conducted  at  these  facilities. 
Continuation  of  these  efforts  are  expected  to  result  in  defini- 
tion of  the  interrelations  between  track  construction  and  main- 
tenance requirements  and  the  vehicle  design  parameters,  required 
for  meeting  ride  vibrations  and  noise  transmission  standards  at 
minimum  costs . 

The  work  described  here  represents  an  initial  effort  to- 
wards meeting  these  objectives.  This  report  describes  work 
currently  in  progress  and  subject  to  revision.  It  is  intended 
primarily  to  provide  information  on  the  current  status  of  in- 
house  TSC  analytic  efforts  conducted  from  November  1971  to  May 
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1.0  INTRODUCTION 


A major  limitation  on  rail  transportation  is  the  tracking 
error  produced  by  the  lateral  vibrations  of  rail  vehicles  and 
wheel  assemblies  in  response  to  track  irregularities.  In  addi- 
tion to  producing  undesirable  vibration  in  the  passenger  compart- 
ment, these  vibrations  are  capable  of  producing  severe  wear  of 
both  the  rails  and  vehicle  wheels.  Under  extreme  conditions, 
these  oscillations  may  result  in  derailment. 

The  most  dramatic  effect  of  the  lateral  dynamics  of  rail 
vehicles  is  the  hunting  phenomenon.  At  low  speeds  (25  to  50  mph) 
in  vehicles  with  lightly  damped  secondary  suspensions  the  phe- 
nomenon is  observed  as  a large  (possibly  violent)  oscillation  of 
the  vehicle  body  which  occurs  at  a characteristic  speed.  At 
speeds  above  and  below  this  characteristic  speed  the  oscilla- 
tions are  reduced.  This  low  speed  body  hunting  normally  does 
not  result  in  a catastrophic  situation  but  does  contribute  to 
excessive  wheel  and  rail  wear  and  passenger  discomfort.  At  very 
high  speeds  (above  100  mph)  the  hunting  phenomenon  appears  as  a 
violent  oscillation  of  the  wheel  assemblies  which  is  limited 
only  by  flange  impact  and  eventual  derailment.  This  oscillation 
begins  at  a critical  vehicle  speed  and  the  motion  is  unstable 
at  all  velocities  above  the  critical  velocity. 


1-1 


Although  the  hunting  instabilities  have  been  observed  for 
many  years,  the  mechanism  of  the  instability  was  first 
fully  understood  only  in  the  past  decade.  Analytical  and  ex- 
perimental studies  conducted  by  Wickens  at  British  Railways 
(Ref.  1)  and  Matsuidaira  at  Japanese  National  Railways  (Ref.  2) 
resulted  in  the  development  of  linearized  analytic  models  which 
could  be  applied  to  determine  the  stability  of  two  axled  rail 
vehicles  operating  on  straight  track.  These  studies  also  re- 
sulted in  the  development  of  roller  rigs  which  were  employed  to 
simulate  straight  track  for  full  scale  stability  investigations 
of  prototype  truck  and  vehicle  designs. 

Although  stability  is  a necessary  condition  for  a new  rail 
vehicle  design  it  is  not  sufficient  to  assure  safe  and  reliable 
performance  on  less  than  a perfectly  straight  track.  In  order 
to  assure  safety,  limit  wheel  and  rail  wear  and  provide  satis- 
factory ride  vibration  characteristics  it  is  necessary  to 
evaluate  the  response  of  the  vehicle  and  wheel  assemblies  to 
rail  irregularities.  The  U.S.  Department  of  Transportation  is 
currently  designing  and  fabricating  a Wheel/Rail  Dynamics 
Research  facility  which  will  include  a dynamic  track  simulator 
capable  of  full  scale  vehicle  tests  to  fully  evaluate  the  coupled 
response  of  rail  vehicles  to  track  irregularities  and  limiting 
rail  conditions. 

This  development  of  experimental  facilities  for  rail  vehicle 
testing  is  being  complemented  by  analytic  studies  being  conducted 
by  the  Transportation  Systems  Center  (TSC)  on  behalf  of  the  Urban 
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Mass  Transportation  Administration,  Office  of  Research,  Develop- 
ment, & Deomonstrations . 

The  analytic  models  developed  as  a result  of  these  studies 
will  provide  a framework  for  interpretation  of  the  results  of 
experimental  investigations  of  the  dynamics  of  rail  vehicles 
and  permit  extrapolation  of  the  test  results  to  new  vehicle 
design  concepts  and  to  prediction  of  the  dynamic  performance  of 
vehicles  under  real  track  conditions . The  analytic  efforts 
will  also  support  the  design  of  the  simulation  equipment  by 
providing  a comparison  of  vehicle  performance  on  track  with  that 
expected  in  the  simulation . The  analyses  will  also  identify  the 
range  of  test  parameters  required  for  vehicle  performance  evalu- 
ation and  identification  of  critical  test  conditions  that  might 
approach  limits  of  the  simulation  capabilities.,  The  results 
obtained  from  simulator  and  track  tests  will  be  used  to  validate 
and  improve  the  analytic  models  to  permit  their  application  to 
new  vehicle  and  component  designs  and  provide  improved  specifi- 
cations for  track  alignment. 

This  document  is  intended  to  provide  a summary  of  the  work 
currently  in  progress  at  TSC  towards  establishing  analytic 
models  of  the  lateral  response  of  rail  vehicles  to  track  irreg- 
ularities and  preliminary  results  of  these  analyses.  Section  2 
develops  the  mechanism  of  the  hunting  instability  through  the 
use  of  simplified  wheelset,  truck  and  vehicle  models  to  develop 
an  intuitive  understanding  of  the  hunting  phenomenon.  It  is 
shown  that  at  low  speeds  the  energy  associated  with  the  lateral 
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vibrations  of  rail  vehicles  is  dissipated  by  a complex  inter- 
action between  the  creep  forces  developed  by  the  different  axles 
of  the  vehicle  through  the  compliance  of  the  vehicle  suspension 
and  wheel  assembly  structures.  The  linearized  models  of  rail 
vehicle  lateral  dynamics  for  two  axled  vehicles  developed  by 
Wickens  (Ref.  1)  and  Cooperrider  (Ref.  3)  are  extended  to  in- 
clude track  irregularities.  It  is  shown  that  by  appropriate 
substitutions  this  model  can  also  be  applied  to  four  axled 
vehicles  for  the  important  special  cases  of  a rigid  truck  frame 
and  of  flexible  truck  frames  with  either  very  low  lateral  stiff- 
ness or  low  torsional  stiffnesses  between  axles.  A heuristic 
formulation  is  presented  for  calculating  the  critical  speeds  of 
rail  vehicles.  This  argument  predicts  low  speed  hunting  will 
occur  when  the  kinematic  frequency  associated  with  the  geometry 
and  compliance  of  the  wheel  assembly  is  in  the  neighborhood  of  a 
vehicle  body  lateral,  yaw  or  roll  natural  frequency.  High  speed 
hunting  will  occur  when  the  kinematic  frequency  approaches  the 
larger  of  the  wheel  assembly  yaw  and  lateral  natural  frequencies. 
This  heuristic  argument  yields  similar  results  for  the  critical 
speed  to  those  obtained  from  the  "resonance"  theory  but  is  based 
upon  stability  arguments  rather  than  treatment  of  the  kinematic 
oscillations  of  the  wheelset  as  a forcing  sinusoid.  The  basis  of 
the  "resonance"  theory  of  lateral  motions  and  the  limitations 
associated  with  this  theory  are  discussed  in  Reference  4.  This 
result  is  essentially  in  agreement  with  the  numerical  results 
obtained  by  Wickens  and  Cooperrider  and  those  obtained  for  a 
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typical  rapid  transit  vehicle  in  Section  2. 

In  Section  3 the  linear  models  developed  in  Section  2 are 
extended  to  include  the  non-linearities  included  by  Cooperrider 
in  Ref.  3 of  flange  impact,  wheel  slip  and  suspension  friction 
and  those  implied  by  profiled  wheels.  The  model  is  extended  to 
include  track  irregularities  and  a simulation  program  is  being 
planned  for  exercising  this  model.  Section  3 also  discusses  the 
extension  of  the  analyses  to  include  the  coupling  of  forward, 
vertical  and  lateral  oscillations  which  are  believed  to  limit 
wheel  adhesion  at  high  speeds. 
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2,0  LINEARIZED  MODELS  OF  LATERAL 
RAIL  VEHICLE  DYNAMICS 


2 . 1 SUMMARY 

The  mechanism  of  lateral  guidance  of  rail  vehicles  by  the 
wheel  conicity  is  reviewed  based  on  linearized  models.  These 
models  are  developed  towards  demonstration  of  the  mechanics  of 
the  hunting  instability  and  the  influence  of  the  compliance 
between  axles  and  secondary  suspension  damping  on  controlling 
this  instability.  The  equations  of  motion  for  predicting  the 
characteristic  roots  describing  transient  performance  and  for 
predicting  lateral  response  to  track  irregularities  are  described. 
The  requirements  of  a relatively  stiff  compliance  between  axles 
is  in  conflict  with  the  guidance  requirements  in  negotiating 
curves.  The  relationships  between  compliance  between  axles  and 
curve  negotiations  are  developed. 

Approximate  relationships  are  developed  for  predicting 
wheel  slip  in  curve  negotiation.  The  lateral  dynamics  programs 
are  checked  against  published  results.  The  implications  of 
these  results  are  applied  for  the  parameters  of  a typical  tran- 
sit car.  Section  3 discusses  the  extension  of  these  models 
to  include  the  non-linearities  associated  with  the  wheel/rail 
interaction  and  real  suspension  system  designs. 
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2 . 2 ISOLATED  WHEELSET 


Rail  vehicles  are  guided  along  a track  by  the  creep  forces 
developed  by  the  interaction  of  the  conicity  of  the  wheels  with 
the  rail  for  small  tracking  errors  and  by  the  forces  developed 
by  the  contact  of  the  flange  and  the  rail  for  large  lateral 
displacement  relative  to  the  track  centerline. 

As  shown  in  Figure  2-1,  when  the  wheel  set  is  given  a lat- 
eral displacement  y the  radius  r^  of  the  right  wheel  increases 
by  ay^  where  a is  the  cone  angle,  and  the  radius  of  the  left 
wheel  decreases  by  ay 


r,  = r + ay 
1 o J 


r„  = r - a 
2 o y 


so  that  the  velocities  of  points  a and  b on  the  axles  are 


(2-1) 

(2-2) 


V 

a 


V 


b 


V 

r 


(ro  + 

ay) 

(2-3) 

(r  - 
o 

ay) 

(2-4) 

This  implies  an  angular  velocity, 

• _ Va  - Vb  _ 2V  ay  _ - V ay 
V 21  ~ 2 rQ£  “ r0£ 


(2-5) 


If  the  vehicle  is  turning  around  a constant  radius  curve,  R, 
and  is  in  a steady  state  condition,  with  no  oscillations  this 
angular  velocity  is: 

i = ^ (2-6) 
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Figure  2-1.  Steering  Action  of  Wheelset  with  Coned  Wheels 
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so  that  in  rounding  a curve  the  lateral  tracking  error  is 

+rQi 

y = -jnf  (2-7 

For  a wheelset  with  a radius  of  17"  a cone  angle  of  0.05  and 
flange  clearance  of  13/16"  on  a track  with  a gauge  of  56.5"  the 
minimum  curve  radius  that  can  be  negotiated  without  flange 
contact  is: 


R 


rc£ 

ay 


17  x 28.25 
0.05  x 13/32 


23,600" 


So  that  the  sharpest  curve  that  can  be  negotiated  without  flange 
contact  is  1970  feet  which  is  about  3/8  of  a mile. 

On  straight  track  the  wheel  set  corrects  for  a lateral 
position  error  by  turning  according  to  Equation  2-5  to  produce 
a lateral  velocity  component: 

y = V ^ (2-8) 


Combining  Equations  2-5  and  2-8  we  obtain: 


y + 5b r Y = 0 (2_9) 

o 

This  is  the  equation  of  motion  of  an  undamped  second  order 
system  with  a sinusoidal  solution.  The  frequency  of  oscillation 
is : 


fk  = b ^7  <2-10> 

Noting  that  dx  = Vdt,  Equation  2-9  can  be  rewritten  as 

y"  + br  Y = 0 (2-11) 

ro*' 
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The  solution  of  this  equation  is  a sinusoid  in  space  with  a 
characteristic  wavelength 


This  undamped  oscillation  of  a single  unrestrained  wheelset  is 
known  as  "kinematic  hunting"  since  the  wheelset  traces  the  same 
path  independent  of  velocity.  f^  and  are  known  as  the 
kinematic  frequency  and  kinematic  wavelength  respectively.  For 
rQ  = 17".  i = 28.25",  a = 0.05,  the  kinematic  wavelength  is: 


for  a speed  of  60  mph,  this  wavelength  corresponds  to  a kine- 
matic frequency  of  1.71  Hz. 

When  the  wheelset  travels  along  an  irregular  track  as 
shown  in  Figure  2-2  the  velocities  of  points  a and  b for  a pure 
rolling  condition  on  the  axles  are: 


(2-12) 


51.3  ft. 


+ a 


(2-14a) 


( 2-14b) 


This  implies  an  angular  velocity 


o 


(2-15) 


For  pure  rolling  of  the  wheelset  on  the  track 


y = v ip 


(2-16) 


2-5 


2-6 


Figure  2-2.  Wheelset  Travelling  on  Irregular  Track 


We  define  the  location  of  the  track  centerline  as 


and  Equation  2-15  becomes 


Combining  Equations 


Y + 


ro£ 

cTv 


\jj  = 6 


2-16  and  2-18  yields: 


(2-17) 


(2-18) 


(2-19) 


Noting  that  dx  = Vdt  equation  2-19  can  be  rewritten  as 

y + 3k2  y"  = 5 (2-20) 

where 

2 = / Ak  \ 2 = ro£  (2-21) 

^k  y 2tt  J a 

For  straight  track,  6 = constant,  the  solution  is  an  oscillation 
at  the  kinematic  frequency  centered  about  y = 6 as  obtained 
above . 

For  a vehicle  travelling 

6 

For  small  angles. 

The  solution  to  Equation 

y 


on  curved  track  of  constant  radius  R, 
2 


= R 


(-5) 


1/2 


(2-22) 


R - 


X 

2R 


2-20  for  this  input  is: 

= (R  + Bk2)  - Is 


(2-23) 


(2-24) 


2-7 


so  that 


(y  - 6) 


ro* 


aR 


(2-25) 


as  obtained  in  Equation  2-7. 

For  sinusoidal  track  irregularities  the  steady  state  solution 
to  Equation  2-20  is 


y = 


(2-26) 


where  A is  the  wavelength  of  the  irregularity  and  f = V/A  is  the 

frequency  corresponding  to  that  wavelength.  For  wavelengths  of 

39  feet  and  A.  = 51.3  ft.  the  center  of  the  wheelset  will  trace 
k 

a path  that  has  an  amplitude  of  about  135%  of  the  amplitude  of 
the  irregularity.  For  irregularity  wavelengths  of  78  ft.  the 
response  amplitude  will  be  1.77  times  the  irregularity  amplitude. 
The  tracking  error  in  following  irregularities  is 


(y  - 6) 


(2-27) 


For  wavelengths  of  39  feet  the  tracking  error  is  235%  of  the 
irregularity  amplitude  and  for  wavelengths  of  78  ft.  the  error 
is  77%  of  the  irregularity  amplitude.  Therefore  with  a flange 
clearance  of  13/16  inch,  the  maximum  irregularity  in  the  track 
centerline  that  can  be  accommodated  without  flange  impact  at  a 
39  foot  wavelength  is  0.346  inch  peak  to  peak. 
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Forces  and  torques  applied  to  the  wheelset  are  resisted  by 
forces  generated  at  the  wheel  rail  contact  which  are  approxi- 
mately proportional  to  the  relative  velocities  between  wheel  and 
rail  at  the  contact  for  relative  velocities  which  are  small 
(less  than  0.1%)  compared  to  the  forward  vehicle  velocity.  For 
the  wheelset  shown  in  Figure  2-2,  the  difference  between  the 
actual  velocities  and  those  predicted  for  pure  rolling  contact 
in  the  forward  and  lateral  directions  are: 


(2-28) 


(2-29) 


(2-30a) 


vbL  = y - vi|j 


(2-30b) 


The  corresponding  creep  forces  are: 


F = -f  — 
aT  T V 


F = -f  — — = + f 
bT  TV  T 


f 


T 


(2-32) 


(2-31) 


(2-33) 


The  net  force  in  the  y direction  is : 


(2-34) 


The  net  torque  in  the  ip  direction  is: 
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M 


- -2ft  (^y  + ^)+  2£t(^) 


(2-35) 


which  for  zero  force  and  moment  yield  Equations  2-16  and  2-18. 

If  the  wheelset  is  given  mass  and  inertia, 

Fy  = my  (2-36) 

M = Cip  (2-37) 

Equations  2-34  and  2-35  become 


my  + 2f 


cijj  + 2fr 


($-*)-» 


^ y + V)=  2f 


aZS 

T r 


Equations  2-38  and  2-39  can  be  written  as 


y ♦ ^ 


o£ 

T r 

o 


” 

— “ 

r* 

~2f  L 

y (s) 

2f  Z2s\ 

+ ' ) 

V (s) 

L J 

u 

2fr 


aZS  (S) 


(2-38) 


(2-39a) 


( 2-3  9b) 


where  S is  the  laplace  transform  variable.  The  transient  solu- 
tion is  defined  by  the  roots  of  the  determinant  of  the  coeffi- 
cient matrix.  This  determinant  is: 


0 4f  f Z _ 

D ■ «l£t  F1  + s +l 

O V 


r2MfT£‘ 
V 


2Cf , 


V 


3 4 

S + MCS 


(2-40 


The  roots  of  this  equation  lie  in  the  right  half  plane  for  all 
values  of  velocity  indicating  that  the  simple  wheelset  is  un- 
stable under  all  conditions  unless  some  additional  restraint 
is  provided. 
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2.3  ISOLATED  RIGID  TRUCK 


For  the  rigid  two  axle  assembly  shown  in  Figure  2-3,  the 
lateral  velocity  of  the  leading  wheelset  is 


y±  = y + hftiji  (2-41) 

and  for  the  trailing  wheetset 

y2  = y - hfti  (2-42) 


The  lateral  creep  forces  are  from  Equation  2-34. 


Fyl 

- -2fL  ( 

} - 

) -2£lm 

i 

V 

fy2 

The  moments  generated 

- -2fL  ($  - * 
at  each  axle 

^ +2fLhft 
are 

i 

V 

Ml  = -2£t 

/aft  (y  + 

V ro 

hftijj ) 

*t) 

2f  aftfi, 

+ 

r 

o 

M2  ~ -2f T 

/aft  (y  - 

' ro 

hfti^ ) 

2f  aft62 

+ 

r 

o 

The  net  lateral  force  is 


Fy  = 4f7 


- (S  - •) 

The  net  moment  acting  on  the  truck  is: 

M = M-j-  + M2  + (F1  - F2)  hft 


M = /-  4£t  ( 


aft  , ft2ii 

— y + — — 
r y V 
o 


+ 4f, 


+ 4f 


aft 


61  + 62 


T r 


h2£2iA 

v J 


(2-43) 

(2-44) 

(2-45) 

(2-46) 

(2-47) 

(2-48) 

(2-49) 


For  no  applied  net  force  or  moment  these  equations  reduce  to: 


y 


Figure  2-3.  Two  Axle  Rigid  Truck  Assembly 
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(2-50) 


Z 

V 


- ip  = 0 


y + 


r l 
o 

aV 


1 + h" 


(2-51) 


If  we  define 


!kt  = (l  + h2  y 


(2-52) 


-T  / 

Equations  2-50  and  2-51  are  identical  to  equation  2-16  and  2-18 
from  a single  wheelset.  The  equation  of  lateral  motion  of  the 
rigid  truck  is: 


y + 


KT 


y = 


61  + S2 


(2-53) 


This  results  in  a kinematic  wavelength 


AKT 


1 + h“ 


K 


(2-54) 


for  h = 1.41  and  A = 51.3  ft,  the  kinematic  wavelength  for  a 
rigid  truck  is  88.8  feet.  The  sharpest  curve  that  can  be 
negotiated  without  flange  contact  for  the  wheelset  used  as  an 
example  in  Section  2.3  is  5,910  ft. 

Since  the  rigid  truck  does  not  provide  any  dissipation 
forces  to  damp  the  kinematic  oscillation  and  has  the  same 
equation  form  as  the  single  wheelset,  the  addition  of  either 
mass  or  inertia  will  result  in  unstable  behavior.  This  apparent 
instability  of  a rigid  truck  was  noted  by  Langer  in  Reference  5 
and  caused  him  to  conclude  that  guidance  through  the  action  of 
the  wheel  conicity  is  impossible  and  that  only  flange  guidance 
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acted  to  keep  a truck  on  a track.  Fortunately,  as  shown  below, 
most  truck  designs  have  finite  stiffnesses  between  axles  which 
results  in  an  energy  dissipation  that  permits  stable  behavior. 

The  response  of  the  rigid  truck  to  irregularities  can  be  obtained 
from  Equation  2-26  by  substituting  A ^ for  AK  and  the  average  of 
6^  and  6^  for  6. 

2.4  FLEXIBLE  TRUCK 

The  simplified  models  given  in  the  previous  sections  do  not 
provide  any  means  of  dissipation  of  the  kinematic  hunting 
oscillation.  As  shown  in  the  following  paragraphs,  a real  truck 
does  provide  dissipation  forces  for  the  hunting  mechanism  by  a 
complex  interaction  between  the  truck  axles  through  the  elastic 
members  connecting  the  two  axles.  The  rigid  truck  and  simple 
wheelset  represent  the  limiting  cases  where  the  compliance  of 
these  elastic  members  is  extremely  small  or  extremely  large. 

For  the  model  shown  in  Figure  2-4,  the  creep  forces  and 
moments  generated  by  each  wheelset  motion  are: 


(2-55) 


(2-56) 


(2-57) 


(2-58) 
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Figure  2-4.  Flexible  Two  Axle  Truck  Assembly 
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If  no  external  forces  are  applied,  the  requirements  of  equilib- 
rium are: 


M1  + M2  + h£ (F1  “ F2)  = ° 

F! 

yx  - y2  = ^ + K~ 

- V ~ Flh£  + M1 

The  substitutions: 


yA  = 


ip  A = 


, y 
, $ 


yl  + y 2 
2 

*1  + ip2 

2 


(2-59) 

(2-60) 

(2-61) 

(2-62) 

(2-63) 

(2-64) 


and  some  algebraic  manipulations  permit  these  requirements  to  be 
rewritten  as: 


$ - f - o 


o£_  - h££L  yA_  _ h££L 

ro  Y fT  V J"T 


TT**  + £t-  ^ + *2 


fL  • fL 

yA  + ^ yA  ^ ^A  - h£^  = 0 

y y 


fT°£  A , , * , f T&2  , A “T 
yA  + i|»A  + — — ip  A = 


fma£  / 6^  - 6 ^ 


K r 
ip  o 


K.  V 

<P 


K , r 

ip  o 


(2-65) 


(2-66) 


(2-67) 


(2-68) 


For  large  values  of  K , yA  = hl<|> . For  large  values  of  , ipA  = 0 
so  that  Equations  2-65  and  2-66  for  a rigid  truck  reduce  to: 

A 

$ - £ = 0 (2-65) 


2-16 


y + 


r a 
o 

aV 


1 + h2  _L  ] ^ = 


61  + 62 


(2-69) 


which  agree  with  Equations  2-50  and  2-51  derived  for  a rigid 


truck.  If  -*  0,  from  Equation  2-67 


ZA  _ 


V 


If  K -*  o from  Equation  2-68, 


ipA  = 0 


a i . , £ * a£ 

— yA  + ipA  = — 
r J V y r 

o o 


h - {2 


(2-70) 


(2-71) 


Substitution  of  Equation  2-70  into  Equation  2-66)  yields 


ct£  - , l2  j a £ 
— Y + v-  <P  = — 
o o 


61  + 62 


(2-72) 


addition  and  subtraction  of  Equations  2-65  and  2-70  and  of  2-71 
and  2-72  yield  the  simple  wheelset  equations  obtained  in 
Section  2.2. 

The  transient  response  of  the  flexible  truck  described  by 
Equations  2-65  through  2-68  is  defined  by  the  roots  of  the 
determinant : 


D = 


S 

V 


, fL  S 

*7^ 


1 + 


fLS 
K V 

y 


o 

hfL 


K 


-1 


i S 
V 


-h£ 


fTai 

Kipr 


1 + 


f £^ 

T S 


:2-73) 
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which  provides  the  characteristic  Equation: 


Some  simplification  results  from  defining 


and  regrouping  terms  to  obtain: 


For  the  limiting  cases  of  K -»■  the  characteristic 

equation  predicts  the  rigid  truck  kinematic  hunting.  For 
K -*■  0 , Ky  0 , the  equation  reduces  to  the  pair  of  single 
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wheelset  kinematic  hunting  equations.  We  also  obtain  the  result 
that  for  small  values  of  either  the  lateral  or  rotational  stiff- 


ness, the  solution  is  the  kinematic  oscillation  of  a single 
wheelset . 

Equation  2-76  may  also  be  written  in  the  form: 


where  for  oscillatory  solutions  of  the  equations  of  motion 


order  system.  3a  and  g^  are  the  associated  damping  ratios.  The 
roots  of  the  characteristic  equation  (Eq.  2-76)  have  been 
calculated  as  functions  of  the  dimensionless  stiffnesses  for 
axle  distances  of  h = 1.41  times  the  gauge  distance  {21)  which 
is  typical  of  truck  designs  and  for  axle  distances  of  h = 10 
which  is  typical  of  truck  center  separations.  These  results 
are  plotted  in  Figures  2-5  through  2-10. 

It  is  seen  that  the  interaction  between  the  two  axles 
through  the  elastic  members  does  result  in  a dissipation  of 
energy  and  that  flexible  trucks  can  be  designed  with  a damping 
ratio  of  0.35  for  the  kinematic  truck  oscillation.  The  kinematic 
wavelength  of  the  damped  truck  would  be  about  midway  between 
the  wavelength  calculated  for  the  rigid  truck  and  that  calculated 
for  a simple  wheelset.  A change  in  the  creep  coefficient  can 
result  in  a dramatic  change  in  truck  behavior.  An  increase 
in  the  creep  coefficient  will  increase  the  kinematic  frequency 
at  a given  speed.  Depending  on  the  truck  design  a change  in 


0 


(2-76a) 


w^  and  w^  are  the  natural  frequencies  of  the  equivalent  second 
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Figure  2-5.  Kinematic  Frequency  of  Two  Axled  Vehicle  Vs.  Non- 
Dimensional  Yaw  Stiffness  (h=1.4) 
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Figure  2-6.  Damping  of  Kinematic  Modes  of  Two  Axled  Vehicle  Vs. 
Non-Dimensional  Yaw  Stiffness  (h=1.4) 
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Figure  2-8.  Damping  of  Kinematic  Modes  of  Two  Axled  Vehicle  Vs. 
Lateral  Stiffness  (h=1.4) 
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Figure  2-9.  Kinematic  Frequencies  of  Two  Axled  Vehicle  Vs. 

Non-Dimensional  Lateral  Stiffness  (h=10) 
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Figure  2-10.  Damping  Ratio  of  Kinematic  Modes  of  Two  Axled 
Vehicle  Vs.  Lateral  Stiffness  (h=10) 


creep  coefficient  can  result  in  either  an  increase  or  a decrease 
in  the  effective  damping  of  the  transient  oscillations. 

As  seen  by  comparison  of  Figures  2-6  and  2-10  the  dissipa- 
tion available  for  damping  transient  oscillations  increase  with 
vehicle  length.  However,  a high  stiffness  between  axles  for  a 
long  vehicle  results  in  a limitation  on  the  ability  of  the 
vehicle  to  negotiate  curves.  For  truck  designs  where  the  axle 
suspensions  are  connected  to  a rigid  truck  body  as  shown  in 
Figure  2-11,  the  effective  stiffnesses  between  axles  are: 


K 


* 


K 

Y 


\ 

+ k yh2H 2 J 


For  the  stiffest  primary  suspension  designs  of  LIM  motor  test 
vehicle : 


k = 52,000  lb/in 
Y 

k.  = 20  x 10^  lb-in/rad 


resulting  in  effective  stiffnesses  for  102"  distance  between 

axles  of  K,  = 10  x 10^  ~ — Y and  K = 3,350  i-^-.  The  dimension- 
ip  rad  y in 

less  stiffnesses  for  a creep  coefficient  of  10^  lb  and  a cone 
angle  of  0.05  are: 


Ky 


Ky3K  _ 3,350  x 100 


10 


0.335 


I 


(*k\2  % ( 100 

\£  / fT6K  \28 . 5 


2 10  x 106 

106  x 100 


1.23 
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Figure  2-11.  Model  of  Two  Axle  Flexible  Truck 
Including  Truck  Body 
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For  a cone  angle  of  0.025  these  non-dimensionless  stiffnesses 
would  be  0.475  and  1.75. 

For  a flexible  truck  travelling  in  a constant  radius  curve 
the  Equations  2-65  to  2-68  can  be  simplified  by  the  substitutions 


* , 

ip/S  0,  ip  = 

~y  (tq)  = 0 = lb  (ti) 

1 S- 

yA  = -hZip 

5 i (tj)  -<$  2 (tf)  =0 


yield  a tracking  error 


y 


ipA 


r Z 

o*1 

aR 


1 + 


h -t 
fT 


1+ 


(£) 


fLfT 

KyKi^ 


h£ 

R 


X + 

VrX^TJ 


(2-77) 


(2-78) 


f_L  /h£\ 

Ky  \ R ) 

yA  = (2-79) 

i j.  a Z fLfT 
rQ  KyK\p 
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The  tracking  error  of  the  forward  wheelset  is 


Yl  = y + yA  = 


r0£ 

aR 


1 + 


h2r 

rT 


1 + 


a£ 


rQ  Kyh£ 


1 + 


a£ 


fLfT 

KyKip 


(2-80) 


for  the  case  of  a rigid  truck  frame  the  tracking  error  is 


?-  fai 

y aR 


1 + h 


2 f L 


(2-8i: 


as  obtained  in  Section  2.3  and  for  an  infinitely  flexible  truck 
frame  we  obtain  the  simple  wheelset  solution  of  Section  2.2. 

For  low  values  of  yaw  stiffness  the  tracking  error  for  the  lead- 
ing wheelset  decreases  with  decreasing  lateral  stiffness.  How- 
ever, for  large  values  of  yaw  stiffness,  a decrease  in  lateral 
compliance  results  in  a larger  tracking  error  than  that  for  a 
rigid  truck. 

The  above  result  has  also  been  obtained  by  Newland  in 
Reference  6 from  the  static  case  of  the  model  shown  in  Figure 
2-11,  with  the  substitutions: 


Tf  _ V 
KiP  ~ 2~ 


k 

K = JL 


y 2 \kijj  + kyh^£  y 

Newland  also  obtains  the  tracking  error  due  to  application  of  a 
lateral  force  P as 
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(2-82) 


,fL/  fT  V 

fh£ 

r0£ 

h_fT\-*-  + rQ  Kyh£/ 

Ph£ 

KiJj  " 1 

aR 

a£  f Lf T 

' 4f 

a i fTifT 

1 rQ  KyKtJj 

1 + rQ  KyKij; 

Slip  of  the  wheels  will  occur  when  the  creep  forces,  required 
to  maintain  the  geometric  relations  of  the  axles  implied  by 
the  elastic  restraint  of  the  assemblies,  exceeds  a critical 
value.  This  force  is  defined  by  the  coefficient  of  sliding 
friction  y as  y W/2  where  W is  the  axle  load.  This  results  in 
a minimum  radius  curve  that  can  be  negotiated  by  the  flexible 
truck  without  wheel  slip.  This  radius  is  obtained  by  Newland 
as : 


, fL 

a£  ^fT 

2 

1/2 

2fTh£ 

1 + h^  fT 

.1  + rQ  VKyh£  /_ 

- p W f Lf  T 

■*-  + ro  KyKi|j 
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For  an  infinitely  rigid  truck  this  minimum  radius  is 

R - ^[l  + h=]V2  <2-8 

Measurements  of  the  limiting  adhesion  coefficient  range  from 
U = 0.15  to  y = 0.3.  The  creep  coefficient  is  of  the  order  of 
150  times  the  normal  contact  force.  For  a rigid  truck  with 
1 = 28.25  inch,  h = 1.41,  the  minimum  radius  that  can  be 
negotiated  without  wheel  slip  is  2,870  feet  for  y = 0.3  and 
5,740  feet  for  y = 0.15.  As  obtained  above,  flange  contact 
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for  this  rigid  truck  will  occur  at  a radius  of  3,410  feet. 
Sharper  curve  radii  may  be  negotiated  by  flexible  trucks  without 
slipping.  For  the  case  of  zero  yaw  stiffness  any  curve  can  be 
negotiated.  For  zero  lateral  stiffness  the  minimum  curve  radius 
that  can  be  negotiated  without  slipping  is: 


20,000  lb  this  minimum  radius  is  565  feet  which  is  considerably 
sharper  than  the  1970  foot  radius  at  which  flange  contact  will 
begin  for  a single  isolated  wheelset.  The  above  equations 
assume  that  the  vehicle  is  travelling  at  zero  speed  around  a 
curve  with  no  superelevation.  When  the  vehicle  is  travelling 
at  speed,  a centrifugal  force  will  be  generated  that  may  or  may 
not  be  compensated  by  the  superelevation  of  the  track.  If  there 
is  no  slip,  the  tracking  error  can  be  calculated  from  Equation 
2-82  by  substituting  for  P the  difference  between  the  centrifu- 
gal force  and  the  lateral  gravity  force  component  produced  by 
the  superelevation.  The  minimum  radius  that  can  be  negotiated 
without  slip  is  calculated  by  reducing  the  effective  adhesion 
force  yW  by  the  same  quantity.  This  is  equivalent  to  reducing 
y by  the  lateral  acceleration  sensed  by  a passenger  in  g units. 
Therefore,  the  smallest  radius  curve  that  can  be  negotiated 
without  slip  for  the  rigid  truck  example  above  with  lateral 
accelerations  of  0 . lg  is  4,300  feet  for  y = 0.3  and  17,300  feet 
for  y = 0.15. 


(2-85) 


Q 

For  a yaw  stiffness  of  10  x 10  lb-in/rad.  and  an  axle  load  of 
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Misalignments  between  axles  will  cause  tracking  errors  and 
loss  of  adhesion  on  straight  track  similar  to  those  that  would 
occur  if  the  vehicle  was  travelling  around  a curve  whose  center 
was  located  at  the  intersection  of  the  center  lines  of  the  two 
axles.  For  a truck  having  an  80  inch  wheelbase  and  an  axle 
misalignment  of  one  milliradian,  this  equivalent  radius  would 
be  6,670  feet  or  1.27  miles.  A misalignment  of  0.1°  would  be 
equivalent  to  travelling  around  a curve  having  a radius  of 
3,830  feet. 

2.5  PARALLELOGRAM  WHEELSET  ASSEMBLY 

In  practice,  many  rail  vehicle  trucks  are  designed  so  that 
the  journal  boxes  of  the  wheel  axles  are  connected  by  rigid 
structural  members,  (which  serve  to  equalize  the  load  between 
axles)  which  rest  on  a rubber  type  pad  mounted  on  the  journal 
boxes  as  shown  in  Figure  2-12  . The  truck  frame  which  houses 
the  drive  motors  and  provides  suspension  for  the  vehicle  body 
is  connected  to  these  equalizing  members  through  the  primary 
suspension  springs.  This  type  of  assembly  suggests  the  model 
shown  in  Figure  2-13.  For  a hard  rubber  pad  (about  70  durometer) 
having  a shear  modulus  of  about  150  psi , an  effective  diameter 
of  3.6  inch  and  a 1/8  inch  thickness,  the  stiffnesses  of  the 
journal  box  connection  to  the  equalizer  bar  on  side  frame  are 
about : 

k = k = 12,000  lb/in. 
x y 

k = 37,000  lb-in/rad. 

a 
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PRIMARY  SUSPENSION 


Ficmre  2-12.  Typical  Assembly  of  Wheel  Axles  and  Equalizer  Bars 
'g  in  Transit  Truck  Designs 
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Figure  2-13.  Model  of  Equalizer  Bar-Wheel  Axle  Assembly 


This  results  in  effective  stiffness  between  axles  referring  to 
model  of  Figure  2-4,  of: 


„ 37,000  00/ 

Ky  ~ — — - = — = 22.5  lb/ 

h £ ( 40 ) 


in 


K ip  = k^£  = 12,000  x 28  = 9.4  x 10  lb-in/rad. 

The  very  small  effective  lateral  stiffness  of  this  type  of  de- 
sign permits  the  model  of  Figure  2-13  to  be  redrawn  as  shown  in 
Figure  2-14.  Since  no  moments  can  be  resisted  at  the  link 
connections  of  the  parallelogram, 

M = M + M2 

M = M 
1 2 


The  equations  of  motion  of  the  wheelset  assembly  are 

Y2 


p = Fl  + F2  = ~2fL^y- 

/a£yi 

M = Mg  + M2  = -2fT(— 


*l/2:  2fL\V_  - ^2 


V 


■)  -2£t( 


a£y2  £ ip2 


+ 


V 


M = -2K^(2 ipa  - - \p2) 

which  may  be  rewritten  as: 


/ a £ 

+ 2fT ( ?o  ^ ^1  + ^ 2 ) 


M = -4Kip(i/ja  - ip) 

/a£  £^g\ 

M = _4fT\^  y + ~v  j 


a£ 

+ 4fT  zr~ 
L o 


(2-86) 

(2-87) 


(2-88) 

(2-89) 

(2-90) 

(2-91) 


These  equations  are  recognized  as  having  exactly  the  same  form 
as  those  of  a simple  wheelset  suspended  from  a frame  having  a 
yaw  stiffness  of  4 Ki/j  as  indicated  in  Figure  2-15. 
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Figure  2-14.  Parallelogram  Model  of  Equalizer 
Bar-Wheelset  Assembly 
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Figure  2 15.  Equivalent  Single  Axle  Model  of  Equalizer  Bar-Wheelset 
Assembly  for  Small  Yaw  Stiffness  Between  Equalizer  Bar 


2.6  SPRING  SUSPENDED  WHEELSET 

For  high  frequency  vibrations,  the  suspension  system  acts 
to  isolate  the  vehicle  body  from  motions  of  the  truck  assemblies. 
As  a result  of  this  isolation  the  vehicle  body  appears  to  be  a 
rigid  frame  moving  at  a constant  forward  velocity  along  the 
track  for  vibrations  at  frequencies  well  above  the  car  body 
suspension  natural  frequency.  Under  these  conditions  the  truck 
and  suspension  can  be  modelled  as  shown  in  Figure  2-16.  For 
truck  designs  in  which  the  wheelset  assembly  can  be  modelled  as 
a single  wheelset,  the  equations  of  motion  are: 

My  + 2fL(f  - + ky(y  - ya)  = 0 (2-92) 

— / ay  Zip  \ — a£6 

cip  + 2fTJ^—  + y-J  + k^dp  - ipa)  = 2f T — (2-93) 


For  high  frequency  motions  the  car  body  displacements  y and  ip 

cl  cl 

approach  zero.  The  characteristic  equation  governing  the 
transient  response  of  this  system  is: 


+ 2S 


3/fT*  f L \ 2/ 

\VC  + mvj  + s 


+ w 


, 2 fL 

+ 2S(wl|J  mV  + Wy  cv 


2 fTb' 


+ IWy  Wg 


4 f uf  T&‘ 

mCV2 


+ 


4fTfL 


mC 


= 0 (2-94) 


The  substitutions: 


2 rQZ 


2fTZ‘ 


, wk  = 


V 

3^ 


Wr 


cs  ' WL  m3 


2f  L 


S - wkSx 
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Figure  2-16.  Wheelset  Suspended  from  Vehicle 
by  Simple  Suspension 


permit  rewriting  this  equation  as: 


1 + 


? O ? 

wip  WL  + wy  WT 


wk2(wT2  + wL2) 

7 5 

Wy  ^ Wip*  + WT  WL 


wlC2Cwy2  + 

2~~T~ 

Wyz  W^Z 


W, 


*)  ■+  W-, 


Wr 


— 5 T 

+ wt  wl 


Si' 


2 SlJ  + 


wk 

~1 


WT  WL  + Wy  Wijj' 


IS}  = 


(2-95) 


Rouths  criteria  gives  the  kinematic  frequency  at  which  hunting 
instability  will  occur  as: 


2 

? WT  ? 

2 _ wtp  + wl2  wy 

~/l  wt2W1  . (V  - wy2)2 

\ wl2/ V + w*2)2 

The  above  criteria  can  be  simplified  by  noting  that 


(2-96) 


wy2  = KyB 
wL2  2 f L 

w ty2  = K^B 

Wip2  2fip£2 


and  that  the  creep  coefficient  is  about  150  times  the  normal 
force , 


wy2  _ KyB 
wL2  150W 


2-40 


and  that  the  lateral  vehicle  natural  frequency  is: 


and  3 ~ 100  inch 


1_  iKy  g 

fl  ~ 2 tt  V W 

wy2  _ ( 2 tt ) 2f126  __  f]_23 

wl2  386  x 150  1500 

and  3 ~ 100  inch 


wL2  15 


is  normally  less  than  2 HZ  and  w^  and  w^  are  normally  of 
similar  magnitudes.  The  magnitude  of  the  term 

» . f »:  - 

\WL  + wiJj 


is  of  the  order 


A < 


225 


225 


.07 


This  term  can,  for  many  designs,  therefore  be  neglected  in 

ematic 
(£  t \ /miL' 


Equation  2-96  and  the  critical  kinematic  frequency  is  given  by: 

2 


wk2  s 


w 


w. 


1 + — i 


m£' 


fL^  C 
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wk 


2 2 
2 _ wip  + wy 


( 2 -97a ) 


The  high  speed  hunting  instability  is  therefore  expected  to  occur 
when  the  kinematic  frequency  of  the  wheelset  is  70%  of  the  larger 
natural  frequency  of  the  wheelset  on  the  vehicle  suspension,  if 
the  two  natural  frequencies  are  different.  The  instability  will 


occur  at  the  natural  frequency  when  the  two  natural  frequencies 
are  equal.  At  low  frequencies  the  terms  k^.  (y  - y a1'  V*  - V 
represent  the  inertia  forces  associated  with  the  car  body  and  a 
portion  of  the  creep  forces  associated  with  the  other  axles.  If 
there  is  no  damping,  these  inertia  forces  tend  to  destabilize 
the  lateral  guidance  of  the  vehicle.  These  destabilizing  forces 
are  reacted  by  the  dissipation  provided  by  the  interaction 
between  vehicle  axles.  The  effective  inertia  forces  are  a maxi- 
mum at  the  vehicle  body  natural  frequencies.  It  is  therefore 
expected  that  lateral  vibrations  will  increase  with  speed  until 
these  natural  frequencies  are  reached  by  the  kinematic  frequency, 
At  speeds  well  above  the  natural  frequencies  the  simplifications 
used  above  are  valid  and  the  suspension  acts  to  stabilize  the 
vehicle.  A characteristic  speed  whose  kinematic  frequency  is  at 
the  vehicle  natural  frequency  evidenced  by  large  body  oscilla- 
tions is  therefore  not  unlikely. 

Predictions  of  the  critical  speeds  based  upon  the  above 
arguments  yield  results  which  are  similar  to  those  that  would 
be  obtained  from  the  "resonance"  theory  described  in  Reference 
4 . This  theory  considered  the  kinematic  oscillation  as  a forc- 
ing function  to  the  vehicle  suspension.  The  "resonance"  theory, 
however,  does  not  account  for  the  unstable  nature  of  the  hunting 
problem  and  can  in  some  cases  result  in  grossly  incorrect 
results.  However,  for  many  vehicle  designs,  these  approxima- 
tions permit  a good  first  estimate  of  the  critical  speeds.  In 
order  to  fully  describe  the  vehicle  behavior  a more  complete 


2-42 


model  is  required.  The  simplest  credible  model  for  predicting 
vehicle  response  and  critical  speeds  is  the  seven  degree  of 
freedom  two-axle  vehicle  model  described  below. 


2.7  TWO  AXLE  VEHICLE  MODEL 

The  two  axle  vehicle  discussed  by  Wickens  (Reference  1) 
and  Cooperrider  (Reference  3)  shown  in  Figure  2-17  can  be  used 
to  model  vehicles  with  either  rigid  trucks  or  flexible  trucks 
whose  behavior  is  equivalent  to  the  single  axle  as  discussed 
above.  This  model  can  also  be  applied  directly  to  study  of  more 
complex  trucks  when  the  natural  frequencies  associated  with  the 
truck  assembly  are  much  larger  than  natural  frequencies  of  the 
car  body  mounted  on  the  secondary  suspension  system.  For  more 
complex  vehicle  assemblies,  the  model  of  Figure  2-17  serves  as 
a basic  building  block  for  construction  of  the  more  complete 
model . 

The  equations  of  motion  of  this  model  are: 


(2-98) 


6 


1 


M 


1 


(2-99) 
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Figure  2-17.  Lateral  Dynamics  Model  for  Two  Axled  Vehicle 
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(2-101) 


(2-102) 


(2-103) 


(2-104) 


where : 

Ky  = 2ky 

cy  = 2cy 

Kip  = 2kip  + Kyh2£2 

cipa  = 2cip  + cyh2£2 

Kq  = Kye2  + Kzb2  = Kye2  + 2kQ 
2 2 2 

c0  = eye  + czb  = Cye  + 2c0^ 

kg  = Gravitational  stiffness  resulting  from  wheel/rail 

normal  forces  (discussed  in  more  detail  in  Section  3). 
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These  equations  can  be  non-dimensionalized  and  cast  into  a 
simpler  form  by  the  substitutions: 
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With  these  substitutions  Equations  2-98  through  2-104  are 
reduced  to: 
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where  the  dot  represents  differentiation  with  respect  to 
dimensionless  time  (t) . 

This  normalized  form  serves  to  reduce  the  magnitude  of  the 
coefficients  in  the  equations  of  motion  with  a resulting  simpli- 
fication in  the  computation  requirements  to  obtain  characteristic 
roots  and  frequency  response  characteristics.  Numerical  solu- 
tions of  the  equations  of  motion  indicate  that  stability  problems 
normally  represent  motions  at  the  wheelset  kinematic  frequency. 
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In  the  normalized  set  of  equations,  displacements  and  their 
derivatives  (e.g.,  ip , ip , ip)  at  the  kinematic  frequency  have  the 
same  magnitude.  This  permits  a rapid  qualitative  evaluation  of 
the  significance  of  the  various  design  parameters  in  the  deter- 
mination of  critical  speeds.  The  normalized  form  also  defines 
the  scaling  laws  required  for  an  experimental  scale  model  of 
rail  vehicle  behavior.  That  is,  for  conducting  a scaled  ex- 
periment, the  coefficients  of  Equations  2-105  through  2-111  must 
be  the  same  for  both  the  model  and  the  full  sized  vehicle. 

Computer  programs  are  being  prepared  to  obtain  the  character- 
istic roots  of  the  above  equations  in  order  to  describe  the 
vehicle  transient  response  and  stability.  A program  is  also  being 
prepared  to  calculate  the  response  of  the  vehicle  to  sinusoidal 
track  irregularities  and  stationary  random  track  irregularities. 

In  order  to  check  the  results,  the  programs  are  being 
applied  to  the  example  vehicle  used  by  Wickens  in  Reference  1 
with  the  following  parameters. 


m = 90  slug  (2,810  lb) 

M = 400  slug  (12,900  lb) 
C = 360  lb-ft-sec2 
J = 12,000  lb-ft-sec2 


K = 2 x 104  lb/ ft 

y 

K,  = 7.11  x 10^  lb-ft/rad 

ip 

Kq  = 2.43  x 106  lb-ft/rad 

O 

k = 9.5  x 104  lb/ft 

g 
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I = 4,000  lb-ft-sec2 
f = 3 x 106  lb 
f_  = 3 x 106  lb 
a =0.4 
H = 2.5  ft 


e = 4 .0  ft 
& = 2.5  ft 

h =2.0 
k.  = 1.06  x 105 
3k  = 3.31  ft 


The  results  of  the  calculations  for  this  vehicle  are  in  agreement 
with  those  obtained  by  Wickens.  The  oscillatory  components  of 
the  solution  are  represented  as  the  natural  frequency  and  damping 
ratio  of  an  equivalent  second  order  system.  Negative  damping 
ratios  represent  unstable  occillations . 

The  results  of  the  computations  for  a vehicle  with  no 
damping  are  shown  in  Figures  2-18  and  2-19.  As  shown  in  Figure 
2-18  the  frequencies  of  oscillation  are  the  vehicle  lateral 
natural  frequencies  and  the  frequency  corresponding  to  the  wheel- 
set  motion.  At  very  low  speeds,  damping  of  the  wheelset  motion 
exists  due  to  the  restraint  provided  by  the  gravitational  stiff- 
ness and  the  stiffness  between  axles.  When  the  wheelset  fre- 
quency approaches  the  lateral  natural  frequency  there  are  large 
changes  in  the  damping  ratio  associated  with  that  frequency.  If 
there  is  no  damping  in  the  secondary  suspension,  the  oscillations 
become  unstable  as  shown  in  Figure  2-19.  The  introduction  of  a 
small  amount  of  damping  (as  shown  in  Figure  2-20)  into  the 
suspension  (3%)  results  in  the  suppression  of  the  unstable 
oscillations  at  speeds  that  do  not  coincide  with  the  character- 
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FREQUENCY (HERTZ 


Figure  2-18 . Characteristic  Frequencies  of 
Two-Axled  Vehicle 
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DAMPING  RATIO 
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Figure  2-19. 


Damping  Ratio  of  Undamped  Two- 
Axled  Vehicle  Oscillations 
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DAMPING  RATIO 
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Damping  Ratio  of  Lightly  Damped 
Two-Axled  Vehicle  Oscillations 
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istic  speeds  at  which  the  wheel  set  frequency  is  equal  to  a body 
lateral  natural  frequency.  As  shown  in  Figure  2-21  the  low 
speed  unstable  behaviour  can  be  completely  eliminated  by  the 
introduction  of  additional  damping.  However,  we  still  note  a 
characteristic  speed  at  which  there  is  a sharp  reduction  in  the 
effective  damping  ratio  which  occurs  when  the  wheelset  frequency 
is  equal  to  the  lower  body  mode  natural  frequencies.  This  is 
evidenced  by  an  increase  in  the  lateral  oscillations  of  the 
vehicle  until  the  characteristic  speed  is  exceeded.  Once  this 
speed  is  passed,  however,  the  lateral  oscillations  decrease  and 
the  damping  is  about  the  same  as  that  under  crawl  conditions. 

At  very  high  speeds,  the  vehicle  behaves  according  to  the 
behaviour  predicted  in  Section  2.5  for  the  spring  suspended 
wheelset.  If  the  gravitational  stiffness  is  included  in  the 
lateral  natural  frequency  w^  Equation  2-96  correctly  predicts 
unstable  motion  of  the  wheelsets  at  all  speeds  above  95.1  ft/sec. 

The  low  speed  hunting  phenomenon  has  been  observed  during 
tests  of  two  R-42  cars  borrowed  from  the  New  York  City  Transit 
Authority  at  the  DOT  High  Speed  Ground  Test  Center.  Figure 
2-23  shows  the  lateral  accelerations  of  the  wheelset  at  speeds 
of  15,  30  and  50  mph,  maximum  accelerations  of  both  the  wheelsets 
and  the  car  body  lateral  motions  were  observed  at  about  30  mph 
at  a frequency  of  about  0.9  Hz.  The  dominant  frequency  of  car 
body  vibration  at  all  speeds  was  about  0.9  Hz.  This  would  imply 
a kinematic  wavelength  of  about  50  feet  which  agrees  with  the 
51.3  feet  predicted  for  conical  wheels  with  a conicity  of  0.05 


2-54 


DAMPING  RATIO 


H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 h 

-20  -10  O 10  20 , 30  40  50  60  70  80  90  lOO  110  120  130  140  150  160  170  180  190  200  210  220  230  240 


- -0.  llOOOO 
-0. 105000 

A 

: A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

-0.1000000 
- -0.0950000 

G 

G 

G 

G 

O 

O 

G 

O 

G 

G 

G 

- -0.0900000 
-0.  0850000 
- -0.0800000 

V 

V 

^7 

^7 

V 

S7 

^7 

V 

^7 

^7 

^7 

.0750000 
.0700000 
.0650000 

- -O.  0600000 
i.  0550000 
.0500000 

- O.  0450000 
. 0400000 

- -O.  0350000 
.0300000 

- -O.  0250000 
. 0200000 
. 0150000 
.OlOOOOOO 
.00500000 


- -O. ' 

- -o. 


o 

+ 


6 


0 


9 


-j-  Wheelset  frequency 
Qwheelset  frequency 
A Roll  frequency 
O Yaw  frequency 
V Sway  frequency 

STABLE 


o. oosooooo 
0.  OlOOOOOO 
0.0150000 


UNSTABLE 


Figure  2-21.  Damping  Ratio  of  Damped  Two- 
Axle  Vehicle  Oscillations 


2-55 


above.  The  measurements  of  lateral  displacements  of  the  truck 
body  indicate  a dominant  wavelength  of  about  50  feet  as  shown  in 
Figure  2-22  (plots  have  common  time  scale)  at  all  speeds.  The 
large  wheelset  accelerations  at  30  mph  imply  that  there  is  an 
increase  in  the  lateral  wear  forces  at  the  wheel/rail  interface 
at  this  speed.  A reduction  in  wheel  and  rail  maintenance  costs 
could  probably  be  achieved  if  this  characteristic  speed  could  be 
avoided . 
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Figure  2-22.  Measurements  of  Lateral  Journal  Box  Acceleration 
and  Truck  to  Rail  Displacements  on  R-42  Cars 
on  Pueblo  Test  Track,  Nov.  1971 


3.0  EFFECTS  OF  NON-LINEARITIES  ON  LATERAL 
RAIL  VEHICLE  DYNAMICS 


3 . 1 SUMMARY 

Section  2 develops  the  mechanics  of  lateral  guidance  and 
lateral  oscillations  of  rail  vehicles  based  on  linearized 
models  of  the  vehicle.  These  analyses  provide  good  results  for 
well  maintained  vehicles  and  roadbeds  on  straight  routes.  How- 
ever, for  rail  vehicles  travelling  on  track  with  large  lateral 
and  vertical  irregularities  and  negotiating  sharp  curves,  the 
non-linearities  that  result  from  flange  contact,  wheelslip, 
suspension  friction  and  mechanical  stops  must  be  taken  into 
account.  In  addition  it  is  necessary  to  consider  the  effects 
of  vertical  wheel  motions  which  may  result  in  a decrease  in 
creep  coefficient  and  a loss  of  adhesion.  The  mechanics  of  these 
non-linearities  and  the  modifications  in  the  equations  of  motion 
resulting  from  them  are  discussed  in  the  following  paragraphs. 

3.2  PROFILED  WHEELSET 

In  Section  2 it  was  assumed  that  the  wheels  were  perfectly 
conical  so  that  the  rolling  radius  was  linearly  proportional  to 
the  wheelset  lateral  displacement.  Although  this  may  be  approxi- 
mately accurate  for  newly  ground  wheels,  wear  of  the  wheel  sur- 
face will  act  to  hollow  out  the  wheels  to  produce  curvature. 

New  rail  vehicle  designs  have  been  using  profiled  wheels  in  order 
to  take  advantage  of  the  increase  in  gravitational  stiffness 


that  results  from  the  action  of  the  normal  forces  at  the  wheel/ 

rail  contact.  The  gravitational  stiffness  is  defined  as  the 

force  per  unit  lateral  displacement  that  is  required  to  move  a 

loaded  wheelset  laterally  in  the  absence  of  friction. 

If  a horizontal  force  Fyg  is  applied  to  the  wheelset  shown 

in  Figure  3-1.  the  wheelset  will  translate  laterally  a distance 

y,  and  tilt  through  an  angle  0g  . This  will  result  in  raising 

the  axle  load  W against  gravity  through  a vertical  distance  6 . 

z 

This  height  change  is: 
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The  work  performed  by  the  force  Fyg  to  produce  a virtual  dis- 
placement 6y^  is  equal  to  the  change  in  potential  energy  of  the 
system 
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6y1  is  equal  to  the  change  in  potential  energy  of  the  system 
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The  gravitational  stiffness  may  be  defined  as: 
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Figure  3 1.  Gravitational  Stiffness  Produced  by  Lateral  Curvature 
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For  conical  wheels 
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we  obtain: 


For  profiled  wheels 
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The  gravitational  stiffness  is  then  for  small  displacements 
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For  large  displacements  the  lateral  force  is  given  by: 
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It  is  seen  from  Equation  3-13  that  the  effective  gravitational 
stiffness  can  be  made  quite  large  without  affecting  the  norminal 
conicity.  This  gravitational  stiffness  adds  to  the  lateral 
suspension  system  stiffness  in  the  analysis  of  Section  2.6  to 
produce  an  increase  in  the  vehicle  and  wheelset  critical  speeds. 
Profiled  wheel  designs  are  currently  being  employed  by  the 
British  and  the  Swiss  to  take  advantage  of  this  effect.  The 
total  force  applied  by  the  rails  on  the  wheelsets  in  the  lateral 
direction  is 
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where  y^  = y - 6 = displacement  from  track  centerline 

is  defined  by  the  wheel  gauge,  track  gauge  and  the  wheel  rail 
profiles.  Changes  in  track  gauge  will  produce  corresponding 
changes  in  6 . 

The  moment  applied  by  the  rails  on  the  wheelset  due  to 
creep  forces  is 
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The  normal  forces  also  produce  a destabilizing  torque  that  is 

given  approximately  by 
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The  net  torque  acting  on  the  wheelset  is 
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For  a profile  that  can  be  represented  by  Equation  3-10, 
Equation  3-18  becomes 
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3.3  TRACK  COMPLIANCE 

The  analyses  given  above  have  assumed  that  the  rails  and 
roadbed  are  infinitely  stiff  compared  to  the  gravitational 
stiffness  and  have  negligible  deflection  as  a result  of  the 
creep  forces.  These  assumptions  are  essentially  correct  for 
small  displacements  from  the  track  centerline.  However,  for 
large  displacements  where  the  flanges  do  come  into  contact,  the 
gravitational  stiffness  becomes  quite  large  and  significient 
track  deflections  may  take  place  Assuming  the  track  has  a linear 
load  deflection  characteristic,  the  deflection  of  the  track 
centerline  is  given  by 
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the  change  in  track  gauge  is  given  approximately  by: 
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where  K is  the  overall  track  lateral  stiffness  and  K is  the 
stiffness  between  opposite  rails  to  a gauge  spreading  force. 
Since  the  track  is  now  capable  of  an  instantaneous  velocity. 
Equation  3-15  must  be  modified  to  include  the  track  velocity  in 
the  creep  forces.  The  creep  force  in  the  lateral  direction  is 
now  given  by: 
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making  note  of 
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and  assuming  the  track  deflects  parallel  to  its  initial  geometry 
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we  obtain 
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The  form  of  Equation  3-23  permits  interchanging  the  space  vari- 
able x and  the  time  variable  Vt  in  the  computation.  Otherwise 
it  would  be  necessary  in  computations  to  maintain  the  variables 
as  functions  of  both  x and  t. 
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3.4  CREEP  FORCES 


The  formulations  given  above  have  assumed  linear  relations 
between  the  creep  velocities  and  the  creep  forces.  The  Johnson 
and  Vermuellen  Analysis  (Reference  5)  which  is  in  good  agreement 
with  laboratory  data  gives  the  following  relations  between  creep 
velocities  and  creep  forces: 
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FR  = V FX  + Fy  (3-26> 

for  F < yN 

K. 

where : 

y = adhesion  constant  (coefficient  of  friction) 

= constant  depending  on  wheel/rail  curvatures 
<j>  = constant  depending  on  wheel/rail  curvatures 
a,b  = major  and  minor  semi-axes  of  contact  ellipse  (pro- 
portional to  N.^^) 

G = Shear  modulus  of  material 
F = longitudinal  creep  force 
Fy  = lateral  creep  force 
v^  = longitudinal  creep  velocity 
Vy  = lateral  creep  velocity 
V = forward  velocity 
N = Normal  force. 
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Equations  3-24  and  3-25  can  be  rewritten  as 
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for  F < yN 

K 


where  and  C ^ are  functions  of  the  geometry  and  material  of 


the  wheel  and  rail. 

The  analyses  of  Section  2 have  assumed  a constant  ratio 
between  creep  force  and  creep  velocity  until  the  adhesion  limit. 
This  approximation  is  valid  only  for  very  small  values  of  the 
creep  forces  and  creep  velocities.  Use  of  a constant  creep 
coefficient  results  in  the  prediction  of  larger  creep  forces 
than  those  which  will  actually  exist.  For  designs  having  low 
primary  suspension  stiffness,  the  actual  value  of  the  creep 
coefficient  has  little  influence  on  the  dynamic  characteristics 
of  the  vehicle  as  long  as  the  dimensionless  stiffness  is  small. 
When  the  dimensionless  stiffnesses  are  of  the  order  of  one, 
however,  the  actual  magnitude  of  the  creep  coefficient  is 
significant . 

The  existence  of  lateral  creep  forces  resulting  from  axle  - 
misalignment  on  rigid  truck  designs  results  in  a decrease  in  the 
creep  coefficients  in  both  the  lateral  and  longitudinal  direc- 
tions . 


3-9 


Field  experiments  on  the  adhesion  coefficients  as  a function 
of  speed  indicate  a decrease  in  adhesion  and  creep  coefficients 
as  a function  of  speed.  Recent  work  by  Paul  (Reference  5) 
indicates  that  this  decrease  in  apparent  adhesion  is  due  to 


oscillations  in  the  magnitude  of  the  contact  stresses  resulting 
from  wheel/rail  vibrations.  This  effect  can  be  approximated  by 
setting 

N = Nq  (l  + a sin  wt) 

where  N is  the  nominal  normal  load  and  "a"  is  the  ratio  of  the 
o 

oscillatory  component  of  the  normal  load  to  the  nominal  normal 
load.  Equations  3-27  and  3-28  become: 
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Since  the  vibrations  associated  with  variations  of  the  normal 
load  will  be  at  high  frequency  (above  40  Hz)  compared  to  those 
of  interest  for  the  investigation  of  lateral  dynamics  (below 
10  Hz),  an  approximate  relationship  between  creep  force  and 
average  creep  velocity  is  obtained  by: 
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The  creep  force,  creep  velocity  relations  implied  by  the  above 
expressions  resulting  from  variations  in  normal  force  are  shown 
in  Figure  3-2. 

The  fluctuations  in  normal  force  result  from  the  vertical 
interaction  of  the  rails  and  vehicle  wheels  in  response  to  track 
and  wheel  and  wheel  irregularities.  A computer  program  has 
been  prepared  for  calculation  of  the  response  of  rail  vehicles 
and  prediction  of  track  deflection  due  to  vertical  rail  irregu- 
larities. This  program  is  described  in  Appendix  B.  The  program 
assumes  the  half  car  model  shown  in  Figure  3-3  which  is  a valid 
approximation  at  short  wavelengths.  Typical  results  of  the 
calculations  are  shown  in  Figures  3-4  and  3-5.  at  relatively 
short  wavelengths  corresponding  to  irregularity  frequencies  in 
the  neighborhood  of  the  wheel/rail  natural  frequencies  large 
amplifications  of  the  irregularity  motions  occur.  The  accelera- 
tions associated  with  these  motions  result  in  fluctuations  of 
the  contact  stress  which  produce  a decrease  in  both  the  adhesion 
limit  and  the  effective  creep  coefficient  as  discussed  above. 
Under  some  conditions  these  changes  may  actually  act  to  stabiliz- 
ing the  hunting  behaviors  of  the  vehicle  which  under  other 
conditions  the  adhesion  may  increase  the  likehood  of  derailment. 
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Figure  3-2.  Half-Car  Model  to  Evaluate  Influence  of  Track  Compliance 
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Figure  3-3.  Creep  Force  Vs  Creep  Velocity  in  the  Presence 
of  Normal  Load  Oscillations. 
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Figure  3-4.  Displacements  of  Car  Body  (y]_)  , Truck  Body  (y2) 
and  Wheels  (y3 , y4 ) of  a Rail  Vehicle  Travelling 
Over  Track  with  Vertical  Irregularities 
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Figure  3-5.  Force  (lbs)  at  Wheel/Rail  Interface  Due  to  1" 
Rail  Irregularity  vs  Frequency  (Hz)  and 
Irregularity  Wavelength 
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3.5  SUSPENSION  SYSTEM  NON-LINEARITIES 

Figure  3-6  schematically  indicates  the  assembly  of  typical 
transit  vehicle  trucks.  For  small  longitudinal  lateral  and  yaw 
motions  the  clearance  between  the  journal  boxes  and  the  truck 
frame  permit  the  wheelsets  to  move  relative  to  the  truck  frame 
to  compensate  for  track  misalignments.  However,  as  the  motions 
become  large  the  journal  boxes  will  come  into  contact  with  the 
truck  frame  resulting  in  a sudden  stiffening  of  the  effective 
primary  suspension.  Large  acceleration  on  braking  forces  will 
produce  a similar  effect.  When  the  relative  motions  between 
the  car  body  and  truck  are  large,  sliding  will  occur  on  rubbing 
surfaces.  There  will  also  be  impacts  with  mechanical  stops  and 
snubbers  which  have  been  designed  to  prevent  excessive  relative 
motions . 
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Figure  3-6.  Typical  Truck  Assembly 
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APPENDIX  A 

COMPUTER  PROGRAM  FOR  PREDICTING  CHARACTERISTIC 
ROOTS  DESCRIBING  TRANSIENT  LATERAL  DYNAMICS 
OF  TWO  AXLED  RAIL  VEHICLES 


Prepared  by 
Leonard  Somers 
Service  Technology  Corp. 
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KENTRON  HAWAII,  LTD 
DOT-TSC  ADP  Support  Services  Project 

PROGRAM  "POLY552"  - AN  INTRODUCTION 

The  computer  program  "POLY552"  described  herein  was  devel- 
oped to  obtain  the  characteristic  roots  of  a set  of  seven  dif- 
ferential equations  which  define  a relatively  simple  system, 
representing  a four-wheeled  railway  vehicle  such  as  a wagon  or  a 
bogie.  This  simple  system  contains  the  essential  physical  con- 
tent of  the  larger  systems  which  must  be  considered  in  investi- 
gations on  stability.  Further  details  on  the  equations  may  be 
found  in  the  appendix. 

POLY552  is  implemented  in  the  MAC-360  language,  an  algebric 
compiler  developed  at  the  MIT  Instrumentation  Laboratory  for  use 
in  fields  such  as  dynamics  and  control  theory.  MAC-360  features 
a three-line  format  corresponding  to  the  3 levels  of  an  ordinary 
algebric  equation,  thus  inputing  equations  in  a form  as  close  to 
the  original  as  possible. 

The  resultant  matrix  of  equations  is  solved  by  POLY552 
through  use  of  a modified  MAC  routine  " AU 232  . PLO YMATRX " , renamed 
"POLY551",  which  evaluates  the  determinant  of  a matrix  with 
polynomial  elements  using  the  MAC  determinant  function.  The 
output  of  the  determinant  function  are  the  polynomial  coeffi- 
cients of  the  system  polynomial. 

A non-MAC  routine  DPRQD  (from  the  IBM  SSP  library)  is 
linked  to  POLY551  to  provide  rooting  of  the  system  polynomial. 
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Routine  DPRQD  uses  a double  precision  ratio  quotient  algorithm 
to  obtain  the  roots  of  a polynomial.  The  generated  complex  roots 
are  then  manipulated  to  yield  the  frequencies  and  damping  ratios. 

The  complete  set  of  frequencies  and  damping  factors  are 
stored  and  then  plotted  using  the  MAC  routine  DCGPLOT . Labelling 
and  lettering  of  the  plots  are  also  included  in  the  plot  portion 
of  the  program. 

POLY552  is  catalogued  in  the  MAC  symbolic -program  storage 
file  at  the  MIT  Draper  Laboratory.  It  can  be  called  into  use  by 
specific  control  cards.  The  control  cards  to  summon  POLY552 
the  necessary  input,  and  the  output  and  plots  expected  are  dis- 
cussed in  the  next  few  pages. 


MO  Is 200 

COO  = SQOMFG 

/ SQOMFG 

SOI F300 

K 

L 

MO  lc400 

S’*  =c.oo 

SCI c 500 

2 

MO  1 8600 

SM  = R 0 0 , S M = 

B 1 5 , SM  =815 

SOI c7CC 

1 41 

61 

R01 °FCO 

r NIC  CF 

FIRST  ROW,  ST  A 

FT  OF  SECOND 

ROW. 

EG  1 8 50 0 

2 

”015000 

SM  = SM  , 

SM  = SM  , 

SM  = 815  H L 

S01E100 

161  1 

162  2 

201  PR1MF 

RO 1 c200 

F NIC  CF 

SEC OWC  ROW,  START  ^F  TFIRQ 

ROW  . 

MO  IS  300 

B 36  = - 2 BETA  R OMEGA 

/OMEGA 

SO  1 c 400 

TWO  T K 

T 

MO lc  500 

C 33  = SQOMFG 

/ SQOM  EG 

SOI c 600 

K 

T 

MO  1 c 700 

SM  =C  33 

SU15«00 

322 

*01co00  SM  = l-e36,Svl  = B36 

S020C00  321  351 

*020100  F.ML  CF  THIRD  ROW,  START  OF  FOURTH  ROW. 

MO  2 C 200  SM  = SM  , SM  = SM 

S02C300  431  321  4?2  322 

RQ2C400  END  CF  P0UPTH  FOR,  START  CF  p I F TH  ROW. 

M020500  SM  = -2  BETA  OMFGA  /OMcGA  , SM  = -SM 

S02C600  601  Y K Y 641  601 


M02C700 

m 

"jC 

II 

SQOMFG 

/SGCMFG  , 

SM 

SO 2 3 300 

6 42 

K Y 

661 

R02CS00 

END  OF  FIFTH  ROW 

, STAR. 

T OP  SIXTH 

ROW. 

M02 1 000 

SM 

-2  8FT 

A OMFGA 

OMFGA 

SO  2 1 IOC 

761 

F I VP 

CT  VF 

M021200 

SM 

-2  BET 

A OMEGA 

OMEGA 

SO  2 1 300 

771 

cOUR 

FOUR 

M02 1400 

SM 

2 °FT  A 

OMFGA  /OMFGA  , 

SO  2 1 500 

SOI 

PSI  K 

PSI 

r.02 160C  rND  L F SIXTH  POw,  START  OF  SEVENTH  ROW. 


SM  F 

641  PRIMR 


/ { SQCMb  (3  H L ) 

K PST  PRIME 

/SQOMFG 
K PS  I 

SM  = SQOMFG  /S00M<=G 
P02  K 


PS  I 


MO  14700 
SOI  4800 

MO  1 4C0G 
SO  1 5000 

MO  15100 
SO  1 5 200 

MO  15300 
SOI 54C0 

MO  15500 
S01560C 

MO  15700 
SO  1 5 500 

M015900 

SOlfcOOO 

MOlfclOC 
SO  1 6200 

BO  1 6 300 

WOl 6400 

MO  16 500 

MO  16600 

MO  16700 

M01 6 tOO 
SO  1 6 900 

E017000 
MO  1 7 100 
SO 1 7200 

MO  17 300 
SOI  7400 

F017500 
MO  1 7 600 
SO  17700 

M017800 
SO  1 7900 

P01 3000 

MO  16100 


A 33  = 

SQR 

- SQOM 

EG  / SQOMEG 

T 

A 

T 

SM 

- 1 t 

SM 

- A3  3, 

SM  = -RQR 

3 00 

320 

350  T 

SM 

= 1, 

SM 

= A 3 3 

460 

430 

SM 

= -It 

SM 

= It 

SM  = E 

600 

SM 

760 

S M 

770 


SM 

RQO 

Sw 


640 

660 

p 

RIME 

-SQOMFG 

/( SQOMFG 

H L 

) 

FIVE 

PS  I 

PRIME 

-SQOMF  G 

/SQOMEG  , 

SM 

= 1 

FOUR 

PSI 

p 00 

-S  QO MEG  /(SQOMFG  F ) 

THRrF  THETA  PRIME 

-SM  , SV  =1 


c 40  °00  c 60 


count  = o 

DU  TO  12  FOR  V = VST ART ( V I NC )VMAX 
PRINT  HCG » TIMEOFDAY,  SP2 
TT  ME 

OMEGA  = V/ BETA 
K K. 


2 

SQOMFG  = OMEGA 
K K 

PRINT  BCG , V,OBFGA  ,SQCN'EG  ,SP2 
K K 


2 

VELOCITY  CMEGA  OMEGA 

K K 

B15  = -2  BETA  R CMEGA  /OMEGA 
ONE  L K L 

START  OF  FIRST  ROW  OF  MATRIX  ELEMENTS. 


BOO  = 1 - 615 
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HJ  I K 

0001 

0002 

00C3 

0004 

0 0 C 5 

0006 

00C7 

0008 

0009 

0010 

0011 

0012 

0013 

0014 


IV  G L t Vp  L 20 


FM  AC 


P AT  F 


72250 


SUBROUTINE  FMAC(CALFIL) 

00U8LE  PRECISION  CALFIL(l),  Q(15),  LE(15)t 
IC  = CALFIH  1) 

CALL  PPRUDUALF  I L ( 2 ),  IC,Q,EE,PPL  , I X,  IFP  ) 

CALFIL(l)  =IER 

CALFIL(2)=IX 

CALFIL(2+IX+1)=Q(IX+1 ) 

no  3 i=  l, i x 

C AL F I L ( 1+2 ) = Q ( I ) 

CALFILI  I +2+  I X+l ) =EE  ( I ) 

3 C AL  F I L ( 2+ I X *2  +1  + 1)  =POL(I) 
f A L F I L ( 4+  IX*3)  =PGL ( I X + l ) 

RETURN 

END 


POL(  15) 
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RO  31300 

S rT 

TWO 

MO  3 1 400 

MODF Y ,1 1 NE ,MAk  K,NMARK , NPOI NTS , Pc9  10^, 

MO  3 1500 

VVM  IN,  VVMAX,  WINCH,  VVT  ICK  , VVLOC  , 

MO  3 1 600 

YAMIN, YAM AX, YAINCH, YAT ICK ,YBLOf , 

R03L700 

SF  T 

T HPER 

MO  3 1800 

MOD  F Y , LINE ,MAP  K ,NMARK , NPOI NTS , PEP IOD, 

MO  3 1500 

VVMIN, VVMAX, W INCH,  VVT  ICK,  VVLOC, 

MO  3 2000 

YAM  IN, YAMAX, YAINCH , YATI CK , YCLQC  , 

R03210C 

SE  T 

pQUR 

MO  3 2200 

MOOFY,  L INF  ,MARK  ,NM  ARK  , NPOI  NTS  , PEP.  IOD 

M032300 

WMIN, WMAX, WINCH, VVTICK, VVLOC, 

M032400 

YAM  IN , YAMAX, YA I NCH, YAT I CK , YDLOC , 

R032500 

SET 

F I V F 

M032600 

MODEY , LINE , MARK, NM ARK , NPOI NTS , PFR  IPO, 

MO  3 2700 

WM  IN,  VVM  AX,  WINCH  , VVT  ICK  , VVLOC  , 

M032R00 

YAMIN, YAMAX, YAINCH, YATICK,YFL PC, 

RO  3 2°00 

S CT 

S IX 

MO  3 3000 

MODEY, LINE ,MAKK , NMAftK , N POI NTS IOC, 

M033100 

VVM  IN, VVMAX,VV INCH, VVT ICK, VVLOC, 

MO  3 3200 

YAM  IN  , YAMAX ,YA I NCH , YAT ICK, YFLOC  , 

R033300 

St  T 

SEVEN 

MO  33400 

MODEY,  LINE , MA R K , NMARK , N PCI  NTS,  PERIOD, 

MO  33500 

WM  IN,  VVM  AX,  W INCH  , VVT  ICK  , VVLOC  , 

MO  33600 

YAM  IN, YAMAX, YA INCH, YAT ICK,YGLOC 

MO  3 3 70  0 

CALL  DCGPLOT (LABL) , 

MO  3 3 8 0 0 

00, 3100, 3100,0.06 

MG  3 3500 

CALL  DCGPLOT { LETR ) , 

M034000 

1.1,  (PI/2), (-1.2), 3. 5, 
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M034  1GG 


22,41  ,21,40,52 ,21 ,37,  19,56,60, 24,  21,41 ,51, 57,2  8, 9<?c 


”034200 

MJ3430G 

R0  3 4400 
P.034500 
R034600 

MO  34 70 C 

MG34800 
MC  34Q00 
>*035000 
MQ35100 
M035200 
3 3300 
MO  3 5400 
mu  3 5 500 
MO  3 5600 
RC3570C 
M .) 35900 
Md 35900 
”036000 
P036  100 
MO 3 6200 
MO  3 6300 
”036400 
HO  36500 
M036600 
”03  6700 
”036800 
R036900 


CALL  OCGPLOKLETR  ) , 

1.1 ,0 ,3, (-1.21 ,53,21 ,35,39, 19, 25, 51,56 ,coq 

START  OF  PLnT  GENF RATION  OF  CAPPING  RATIOS 

NS<-'TS  = 7,L  INE  = C,  VARK=1  ,KMARK=1  , NPO I NT  S=  LOOP,  PFR  I0P=}.5 

VVMIN=G,VWAX  = 220,  VVINCH=5.5,  VVTICK=  10  , VV  LOC=  8 000  0 

YAM  I N=— G . 10,YAMAX  = C.10, YAINCH  = 7,Y ATICK =0. 00 5 , Y AL OC =V VL OC  + 1 

MODE X =0 » MCDEY  =2  » Y 3L0C  = V VLOC+3  , YCL OC  = VVLOC  +5,YDL0f.  = VVLnC  + ^ 

YFLOC=VVLCC+c»  YFLOC= V VLOC+1  1,  YGL3C  =VVLOC+  1 3 

CALI  OCGPLOTt PLOT) ,NSFTS, 

SET  r NIF 

MODE  X ,LINF  ,MA K K , NMARK  , KPOI  NTS  , PFR  I CD , 

VVMIN,  WMAX.VV  INCH,VVTICK,VVLOC, 

YAM  IN, YAM AX, Y A INCH, YA  TICK ,YAL°r , 

SET  TWO 

MODFY , L IMF ,MA°K, NMARK, MPOI NTS, PE P I OD, 

VVMI N , VVM  AX , V V I NCH , VVT ICK , VVL  OC , 

Y A” I N, YAM AX, Y A INCH, YATICK, YBLOC, 

SET  T HP  FF 

MOD  - Y ,LINE, MARK, NMARK, NPOINTS , PFR  100, 
VVMIN,VVMAX,VVINCH,VVTICK, VVLQC , 

YAMIN, YAM AX, YAINCH, YATICK, YCL Of, 

Sri  FOUR 

MODE Y , LI NE , MA R K , NMC RK , N PCI  NTS , PEP IOD, 

VV”  IN,WMAX,VVINCH,  VVT  I C-K  , V VLDC  , 

YAM  I N , YAMA  X , YA  I NCH  , Y AT  I f K , YDLOC. , 

SET  FIVE 
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MO  3 7 GOO 

MODEY  , LI  NE  ,MAE.K,  NMARK , N PCI  NTS , PER IOD, 

MO  3 7 100 

VV«  IN, VVM AX, VV INCH, VVT ICK .VVLHC , 

MOB  / 200 

YA  M I N , YAM A X , Y A I NCH , YAT I CK , Y FLOG , 

R037300  SFT 

SIX 

MO  3 7 40  0 

MODE  Y , L 1 NE ,MAP  K, NMARK ,NPOI NTS , PER  I CO, 

MO  3 7 500 

VVMIN,VVMAX,VV TNCH,VVT ICK,VVLGC, 

M037600 

YAMIN,YAMAX,YAINCH, YATICK, YFLOC , 

P03770G  Sr  T 

Sc  VFN 

MO  3 7 R 0 0 

MGDEY , L I N fc , M A P K, NMARK, NPCI NTS , PFP I CO , 

MO  37900 

VVMIN  ,VVMAX,VVINCH,VVTICK, VVLOC, 

M03Q000 

YAM  IN  , YAM AX, Y A INCH, YATICK, YGLOf 

M 0 3 8 1 0 0 

CALL  DC.  OP  LOT  (LABL)  , 

M038200 

00,3100,3100,0.06 

MO  3 P 300 

CALL  DCGPLCT (l ETR ) , 

MO  3 PACO 

1. 1 , ( Pl/2 ) , (-1.2 ) ,3.5  , 

M03P500 

20,  l7, 36, 3°, 25, 37, 2 3, 13, 22,17,  IP, 51, 38, 41, 9CE 

M03P600 

CALL  DCGPLCT ( L°TP } , 

MO 3 F 700 

1 .1 ,0 ,3, ( -1.2 ), 33, 21,  35, 3 3,  19, 25, El, 56,°°° 

MO3R80O 

CALL  DCGPLCT, (-1) 

R 03  P 900  c NO 

0<=  PLOT  GENFP  AT  I ON 

MO  3° 000 

CAL  L {NON-MAC  ) GC  , (-1  ) 

M03S1GG 

START  AT  BEGIN 
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MG2  ° 300 

PRI NT  HDG  » LC  OP , S PI 

V02P400 

RUN  NO. 

MQ2  E 500 

UCHIT  = 2000  + 5C  LOOP 

M028600 

SET  FILE  READ  OOUT 

VO?  c 700 

PILE  ? FAD  V,X 

MO 2 3 BOO 

PRINT  HDG » V » X 

MO  2 890  0 

VELOCITY  NO.  ROOTS 

M029000 
SO  2C  100 

FILE  READ  CMG A TO  CNG A 
1 X 

Mu29200 
SO  2 G 300 

PRINT  MSG  t UMGA  TC  CMGA  tSP2 
1 X 

M02  c400 

P F A L ROOTS  IN  HERTZ 

M02c500 

S029600 

FILF  READ  BFTA  TO  BETA 

1 X 

M029700 

S029PCC 

PRINT  MSG.BETA  TO  BETA  ,SP2 
1 X 

MO  2 3G00 

IMAGINARY  ROOTS  IN  HERTZ: 

MO  3000 J 

100 

DUMMP  = 0 

P03C100 
RG  30  200 

START  0 

c PLOT  GENERATION  OF  PEAl  ROOTS 

MO 3 0300 

NSE  TS  = 7,L I NF=0 ,MARK  = 1 , NMARK=1 » NPU I NT S = LOOP , PEP  100= 15 

MO  3 0900 

VVM  IN=0,VVMAX=220, VV INCH=5. 5, VVT I f K = 10 , VVLOC =7 5000 

M03C500 

YAM  I N = 0 » Y A MA  X = 12  ,YAIFCH=7  » Y AT  I CK=  1 ,Y AL 0C  = VVLOCFl 

M03G600 

M0ncX=0, M0DtY=2 , Y3LOC=VVLOC+3, YCL0C=VVL0C+5» YOLnC=VVLI 

MO  3070  0 

YFL  0C  = VVL  0C+  Q » YF  LOC  = V VLOC+ 11»YGL0C=VVL  0C+ 1 3 

M 0 3 0 P 0 0 

CALL  DCGPLGT (PLOT) »NSFTS» 

R030900 

SET  ONE 

MO  3 100  0 

MUDE:X  ,LI  Nfc  ,MARK,NMARK  ,N POINTS  , PER  IOD, 

v03 1 1 00 

VVM  IN, VVMAX.VV  INCH, VVTICK, VVLOC , 

M03  1200 

YAM  1 N , YAM  A X',  YA  I NCH  , YA  T I C K , Y ALCC  , 
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M02  5000 

WO’E  KING  VECTOR  IS: 

M02510G 
SO 2 5200 

PRINT  MSG»Q  »SP2 

X+l 

MQ2530G 

MAXIMAL  RELATIVE  error  IS: 

MO  2 540  0 

SET  r II  E W P I T r Q3 

’>*02  5 50  0 

r I L F WE  1 T E V , X 

M0256Q0 

00  T O 22  FOR  W = 1 (1  ) X 

EO  2 5700 
M025800 
S 0 2 5 ° 0 0 

2 2 

OlMGA  =OMFGA  SORT  ( Q 4”  )/(2  PI) 

W K vs  W 

MO 2 6 CO C 
S026  100 

22 

c I L L WRITE  CMC,  A 

W 

M 0 2 fc  2 0 0 

CO  TO  34  E OR  W = L(1)X 

MO 2 6 300 
S026400 

BETA  = — 0 M T G A Q /(CMC,  A 2 PI) 

W K .A  A 

MO 2 6 500 
SO 2 6600 

34 

PUP  WRITE  3ETA 

W 

MO 2 6700 

bCT  FILE  R f A r Q3 

M026RO0 
S02  6°0  0 

TILE  k FAD  V,X, 0 MG A T J C MG A , 3PT  A 

1 X 

MG27000 

5FT  FILE  WRITE  JC 

MO 27 100 

FILL  WRIT  f V 

MO  27 200 
SO  27  300 

FILE  WRITE  QMG  A TG  C MG  A 

I X 

MO 2 7 400 

0 C = 0 C + I d 

MO  2 7 500 

SET  FILE  WRITE  QD 

MO 2 7 600 
5027700 

FILE  WRITE  V,  6 1 T A T P B^TA 

1 X 

M027800 

0 0 = Q D + 1 5 

M027600 

12 

DUG  MR  = 0 

R 02 E 000  E 
R 0 2 “ 1 0 0 

CE  FI!  f WE'  Ilf*  START  Of  pj|_r  RE&Q  aND 

M028200 

10  TO  100  FQR  L HO P = 1(1) COUNT 

PET  A 

X 
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M02  1700 
S02  1800 

M02  1900 
SO 2 200  0 

”022100 

S02220G 

R022300 

MO 2 2 400 

”022500 
SO  2 2 80  0 

”022700 

M022P00 

MO  2 2C00 

MO  2 2000 
S022100 

M023200 

M023300 
S02  3400 

M022500 

”023600 

S023700 

M023R00 
S02  3900 

M024C00 

M024100 
SO  2 4 200 

M024300 

S02440C 

M024500 

M024600 

”024  7CC 

MO  24800 
S024900 


s y 

= -2  BETA 

CNEGA 

OMEGA  /(SQOMEG  E 

P01 

THREE  THREE 

K THFTA  PRIME 

S” 

= - S M , 

SM  =2  BETA 

OMEGA  /OMEGA 

-41 

= 01 

E 6 1 

THFTA  K THFTA 

S” 

= SQOMEG 

/ SOC  ME  G 

5 62 

K 

THETA 

!'  NIL.  CE  SEVENTH  ROW. 

CALL  PI  LY551  ( SUB)  ,C,1 ,7,2 

RESUME  N, GA INN * GA I NDC  » POL Y TO  POL  Y 

0 N 

COUNT  = COUNT  + 1 

Q A = 1000  + 25  COUNT,  G 3 = 2000  + 50  COUNT 

SFT  R 1 1 E vs R I T r-  g A 

FILE  WRI7R  V , U , POL  Y TO  °GL  Y 

0 N 

I C = N + 1 

C OE F =0 
0 

on  TO  11  FOR  A =1 1 1 ) ( N+l ) 

U r,RcF  = POLY 

A N+ l -A 

PRINT  MSG , COE F TG  COEF  ,SP2 
1 N+l 

COFcF I C IF NTS  IN  ASCENDING  ORDER  AP E: 

CALL (NONM AC, RESIDENT )GG, IC ,COFF  TO  COEF 

1 N+l 

RESIJMR  I FR » X , 0 TO  0 , EE  TO  FF  , POL  TO  POL 

1 X+l  1 X 1 X+l 

I F 1 R R NZ,  GO  TO  13,  OTHERWISE  GO  TO  2C 
13  C0UNT=C0UNT-1 
GO  TO  12 

20  PRINT  MSG,  POL  T(j  POL  ,SP2 
1 X + l 
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MO  1 1 OOO 

SQOMEG 

= KG/MS 

SOI  1100 

G 

M0 1 1200 

SOD MEG 

= KA/  C 

SOI  1300 

A 

M01 1400 

SCO  MEG 

= k r* 

SO  1 1 500 

Y 

Y 

M01 1600 

OMEGA 

= SORT (SCCMEG  ) 

SOI  1700 

Y 

Y 

MO  1 1 RQC 

PET  A 

= C / ( 2 M CMEGA  ) 

SOI  1 900 

Y 

Y Y 

MO  12000 

£ = 

E / B t T A 

SOI  2100 

PRIME 

K 

MO  1 2 200 

L 

L/PF  TA 

SCI  2300 

PR  I ME 

K 

MO  1 2400 

R 

OMEGA.  /CMFGA 

SO  1 2 500 

L 

ONE  L 

EO 12600 

2 

MO  1 2700 

SQR 

R 

SO  1 2 BOO 

L 

L 

MO  1 2900 

R 

CMr  G A /CMFGA 

SO  13000 

T 

TWO  T 

E013 100 

2 

MO  1 3200 

SQR  = R 

SOI  2300 

T T 

MO  13400 

CC=  75000 

M013500 

00=80000 

RO 13600  MATRIX 
RO 13700 

cor: ST  ANTS 

MO  13800 

AOO  = SQR 

+ SQOMEG  /SQOMEG 

SO  1 3°00 

L 

G L 

MO  1 ^tOOO 

A 15  = —SQR 

c 

S014 100 

L PRIMP 

MO  14200 

SM  = AOO, 

S M = -1,  S M = A 15, 

SM  = A 15 

SO  14300 

0 

2 C 40 

60 

E 0 144  00 

2 

MO  14500 

SM  = SM 

, SM  = -1,  SM 

-SQR  H L 

SO  1 4600 

160 

0 180  200 

L PRIME 
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M0Q7200 
SO G 7 300 

OMEGA 

TWO 

M0C7A00 
SOC7  50C 

BETA 

TWO 

MO  07  600 
S007700 

SQOMEG 

PSI 

MO  07  f 00 
S 0 C 7 c 0 0 

OMEGA 

PSI 

Moceooo 

sooeioo 

BETA 

PSI 

MO C E 2 00 
S00E300 

SQOMEG 

THET 

MO  C 8400 
SOO  ft  500 

OMEGA 

THETA 

MOC  EfcOO 
SO0P700 

BFT  A 

THETA 

E0CB800 
MO  08900 
S00°000 

SQOMEG 

THRE 

MOC9IOO 

S009200 

OMEGA 

THREE 

EQC9300 

M00c400 

S0CC50C 

EFT  A 

THREE 

M00c600 

S00970C 

SQOMEG 

FOUR 

M0C9600 

S009900 

OMEGA 

FOUR 

M010000 

S010100 

BET  A 

FOUR 

F010200 
MO  1 C 300 
SO  10400 

SQCMEG 

FIVE 

MO  10500 
S01C600 

OMEGA 

FIVE 

EO 1 0700 
MO  1 C POO 
SO  1 0°00 

BETA 

FIVE 

= S G RT ( S COMFG  ) 

TWO 

SC  / ( 2 C OMEGA  ) 

PS  I TWO 

= K /J 
PSI 

= SQP.T  ( SQCMEG  ) 

PSI 

= C / ( 2 J OMEGA  } 

PSI  PSI 

= K /I 

A THETA 

= SORT (SQCMEG  ) 

THETA 

= C / (2  I OMEGA  J 

THETA  IHFTA 

2 

= K E /I 
E Y 

= SORT { SQOMEG  ) 

THREE 

2 

= C E / ( 2 I OMEGA  ) 

Y 1 HR  r F 

= 2 SK  /J 
PSI 

= S QRT ( SQCMEG  ) 

FOUR 

= C / ( J OMEGA  ) 

PSI  FOUR 

2 2 

= K H L /J 

Y 

= SQR  T { SQCMEG  ) 

FIVE 

2 2 

= C H L / ( 2 J OMEGA  ) 

Y FIVE 
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MO  00  100 
S000200 

MO  C C 300 

M00G400 

S00C500 

MO  00600  BEGIN 

M000700 

MOOCBOO 

MO  00Q00 

MO  C 1 COO 

MO  0 1 1 0 0 

MO  01200 

MO  0 1300 

MOO  1400 

MO  C 1 500 

M001600 

M00170C 

M001800 

MOO  1500 

MO  02000 

MO  0 2 100 
S002200 

E002300 
MO  02400 
SO  02  500 

F0C2600 
MO  0 2 700 
SO  02800 

E002900 

M0C3C00 

S003100 

M003200 

S003300 


COMMON! P0LY55 1 ) ,SM 

224° 

I NO  E X N,W,A,X 

RESERVE  POLY  ,COFF  ,Q  , EE  ,POL  , OMG  A , BE  TA 

15  16  16  16  16  16  16 

PEAL)  VSTART ,V INC,VMAX 

READ  ALPHA, 6»C,E,FL, FT 

REAP  H,I, J,KA,KG,L 

PEAD  M,MS,RO 

REAP  Ci,C2,C3,Kl,K2,K3 

PRINT  HDG,  VSTART, VINC,VMAX,SP2 

V START  INCP  V END 


PRI NT  HDG, ALPHA, B,C,E,FL,FT,SP2 
ALPHA  BCE 

PRINT  HDG » H » I ,J,KA,KG,L,SP2 
H 1 J KA 

PRINT  HDG, M, MS , RO, SP2 
M MS  RC 

PRINT  HDG,C1,C2,C3,K1,K2,K3,SP2 
Cl  C2  C3  K 1 

K = 4 K2 

Y 

2 

SK  - 2 K1  B 

PSI 

2 2 

K = K HL+2SK 

PSI  Y PSI 

2 2 

K = K E+4K3B 

THETA  Y 

C = 4 C2 

Y 


FL  FT 


KG 


K 2 K3 
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EO  03400 
MO  03  500 
S003600 

2 

SC  = 2 Cl  8 

PSI 

E003700 
M0C3800 
SO  03900 

2 2 

C = C HL  + 2SC 

PSI  Y PSI 

E004000 
MO  04 100 
SO  04  200 

2 2 

C = C E + 4 C3  e 

THETA  Y 

MO  04300 
SO  04400 

PRINT  HDG,K  » SK  ,K  ,K  ,C  ,SC  ,SP2 

Y PSI  PSI  THETA  Y PSI 

M004500 

S004600 

K SK  K K C 

Y PSI  PSI  THFTA  Y 

M0C4700 

S0C4800 

PRINT  HDG  t C ,C  ,SP2 

PSI  THETA 

MO  04  900 
S005000 

C C 

PSI  THETA 

MO 05 100 
S005200 

SQRETA  = RO  L/ALPHA 
K 

M005300 
SO  05400 

BETA  = SQR  T ( SQB  FTA  ) 

K K 

M005500 

S005600 

SQDMEG  = 2 F L/ ( MS  BETA  ) 
L K 

MO  0 5700 
S0C5800 

DMEGA  = SQR  T ( SQCMEG  ) 

L L 

M005900 

S006000 

SQOMEG  = K / ( 2 MS) 

ONE  Y 

M006100 

S0C6200 

OMEGA  = SORT ( SCCMEG  ) 

ONE  ONE 

MO  06  300 
SO  06400 

BETA  = C / ( 4 MS  OMEGA  ) 

ONE  Y ONE 

F006500 
MO  C 6600 
S006700 

2 

SQOMFG  = 2 FT  L / ( C BET  A ) 
T K 

MO  C 6 800 
S006°00 

OMFGA  = SQRTtSQCMEG  ) 

T T 

M007C0G 
SO 07 100 

SQOMEG  = SK  / C 

TKO  PSI 

sc 

PSI 
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CONTROL  CARDS 


The  following  control  cards  are  required  to  invoke  POLY552 


from  storage: 

1)  //  JOB 


This  card  enables  job  accounting  and  has  the  form: 


Col  1-2: 
Col  3-8: 
Col  14-16: 
Col  18-21: 

Col  23-28: 


// 

job  name 
JOB 

job  charge  number  [to  be  obtained  from  Data 
Services] 

User's  name 


The  following  parameters  may  also  be  inserted: 

a)  Region 

b)  Time  estimate  in  minutes 


Example  JOB  card 

Card  Column  1 bEblank 


2 ) //PGPOLY552bbJOB4 37  9 , SOMERS , REGION=22  OK , TIME=6 

2)  //bbbbbbbEXECbMACOMPIL 

Indicates  that  the  MAC  compile-time  monitor  is  to 
process  a compilation. 

3)  / / SYSINbDDb* 

This  card  informs  the  Operating  System  that  the 
following  cards  will  be  processed  by  the  program 
begmg  executed. 

4)  *bbbbbbbMAC*POLY552 

Compilation  Control  Card 

5)  /* 

Last  card  of  a job  step. 

6)  //bbbbbbbbEXECbFORTGCL 

Causes  compilation  and  linking  of  Fortran  subroutines 

7 ) / /LKED . SYSLINbDDbDSN=MAC .MAC4TRAN ,DISP=SHR 

/ /bbbbbbbbbbDDbDSn=AAAALIN , DISP= (OLD , DELETE ) 

causes  new  linking  routine  to  be  included  in  the 
load  module  into  which  the  Fortran  subroutines  have 
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8 ) //dddddddEXECbMACRUN 

This  card  indicates  that  the  MAC  monitor  is  to 
process  the  run  steps  in  the  job. 

9 ) //R . EXTL IBbbDDbDSN=  &AAAAMOD , DISP=OLD , DELETE ) 
//bbbbbbbbbbDDbDSNAME=SYSl .MACEXT ,DISP=SHR 

Defines  data  set  in  which  non-MACload  modules 
are  located  at  MACRUN  collection  time. 

10)  //R.FT06F001bDDbSYSOUT=A 

Messages  produced  by  the  Fortran  error  monitor 
are  written  in  the  data  set  described  on  the 
DD  card  corresponding  to  Fortran  unit  6 . 

11 ) //R . CATLGbbDDbbDSN=SYSPLOT . CATLGB , DISP=SHR 

This  card  alerts  the  plotter  and  should  be  used 
only  if  a plot  is  desired. 

12)  *bbbbbbbbRUNPOLY552 

Triggers  POLY552  from  storage  and  transfers  control 
to  it . 

This  card  is  followed  by  the  data  cards. 
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INPUT 


The  input  parameters,  which  are  outlined  below,  are  pro- 
jected into  the  deck  in  the  order  indicated.  The  comma-delimited 
format  is  utilized;  this  format  dispenses  with  the  need  of  para- 
meter definition  by  separate  format  statements,  but  does  require 
the  parameters  to  be  separated  from  one  another  on  a data  card 
by  commas.  No  comma  should  follow  the  last  parameter  on  a data 
card.  Parameters  may  be  punched  on  any  of  the  80  columns. 

FIRST  DATA  CARD  [3  Parameters] 


Order  of 

Parameter  Name  Description  Units 


1st 


VS TART  Initial 

Velocity 


ft/sec 


2nd  VINC  Incremental 

Value  of 
Velocity 

3rd  VMAX  Final 

Velocity 


ft/sec 


ft/sec 


SECOND  DATA  CARD  [6  Parameters] 

Order  of 

Parameter  Name  Description 

1st  ALPHA  Effective  conicity 

of  a wheel,  of  a 
wheel-set  in 
question 


2nd 


B 


Half  length 
of  contact  area 
in  direction 
transverse  to  the 
direction  of  rolling 


Units 


ft/ft 


ft 
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3rd 


C 


4th  E 


5th  FL 


Moment  of 
Inertial  of  a 
wheel-set  in 
yaw 

Height  of  body 
centre  of  gravity 
above  wheel-set 
axle  center-line 

Longitudinal 
creep  ccefficient 


6th 


FT 


Transverse 
creep  coefficient 


THIRD  DATA  CARD  [6  Parameters] 


Order  of 


Parameter 

Name 

Description 

1st 

H 

Constant  which 
when  multiplied 
by  "L"  gives  "HL" 
the  semi-wheel 
base 

2nd 

I 

Moment  of  Inertia 
of  body  in  roll 

3rd 

J 

Moment  of  Inertia 
of  body  in  yaw 

4 th 

KA 

5 th 

KG 

6th 

L 

Half -distance 
between  contact 

points  in  a lateral 
direction 


slug  ft2 

ft 

lbs 

lbs 

Units 

slug  ft2 

slug  ft2 

f t-lbs/rad 

lbs/ft 

ft 
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FOURTH  DATA  CARD  [3  Parameters] 


Order  of 

Parameter 

Name 

1st 

M 

2nd 

MS 

3rd 

RO 

FIFTH  DATA  CARD  [6  Parameters] 


Order  of 

Parameter  Name 

1st  Cl 

2nd  C2 

3rd  C3 

4 th  K1 

5th  K2 

6th  K3 


Description 

Units 

Mass  of  body 

slugs 

Mass  of  wheel-set 

slugs 

Wheel-tread  circle 
radius,  wheel-set 
in  central  position 

ft 

Description 

Units 

Longitudinal 

suspension 

damping 

coefficient 

lbs/fps 

Lateral 

suspension 

damping 

coefficient 

lbs/fps 

Vertical 

suspension 

damping 

coefficient 

lbs/fps 

Longitudinal 

suspension 

stiffness 

lbs/ft 

Lateral 

suspension 

stiffness 

lbs/ft 

Vertical 

suspension 

stiffness 

lbs/ft 

The  significance  of  the  various  parameters  and  their 
diagrammatic  representation  is  given  in  the  Appendix. 
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DECK  STRUCTURE 

//JOB 

//  EXEC  MACOMPIL 
//  SYSIN  DD* 

POLY552 
MAC  PROGRAM 

/* 

//  EXEC  FORTGCL 
//  FORT.  SYSIN  DD* 

FMAC 

FORTRAN  SUBROUTINE 
/* 

//  LKED . SYSLIN  DD  DSN=MAC.MAC4  TRAN,  DISP=SHR 
//  DD  DSN=&AAAALIN , DISP  = (OLD,  DELETE) 

//RUN  EXEC  MACRUN 

//R.EXTLIB  DD  DSN=&AAAAMOD , DISP=  (OLD , DELETE ) 

DD  DSNAME=SYS1 . MACEXT , DISP=SHR 
//R.FT06F001  DD  SYSOUT  = A 
//R.CATLG  DD  DSN=SYSPLOT .CATLGB,DISP=SHR 
//SYSIN  DD* 

* RUN  POLY  552 

DATA 
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Sample  Data 

Card  Col  1 Col  9 

5,  20,  220 
0.4,  3,25,  360,  4, 
2,  4E3,  12E3,  1787 
4E2,  90,  1.75 
0.0,  0.0,  0.0,  5E3 


3E6,  3E6 
949E2,  2.5 

5E3 , 5E4 
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PLOTTING 


The  Plot  Control  Card  should  be  used  if  a plot  is  required; 
otherwise  remove  this  card. 

POLY552  will  provide  2 framed  "7.5  x 9"  plots  for 

1)  System  Frequency  (on  the  ordinate  axis) against 
velocity . 

2)  Damping  Constant  (on  the  ordinate  axis)  against 
velocity . 

Points  are  marked  by  small  dots. 

The  program  in  its  present  form  will  only  provide  for  plot 
values  of  velocity  in  ft/sec  from  0 to  220  against  corresponding 
values  of  frequency  in  hertz  from  0 to  12,  and  Damping  Constant 
from  -0.10  to  +0.10. 

For  the  sample  data,  the  following  plots  were  produced. 
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OUTPUT 


For  each  velocity  value  POLY552  will  supply  system  frequencies 
and  corresponding  damping  ratios. 

The  two  identical  sets  of  values  observed  is  due  to  the 
fact  that  complex  roots  of  the  system,  on  which  computation 
is  performed,  occur  in  pairs. 

For  the  sample  data,  the  following  is  the  output.  Only  one 
data  set,  of  the  identical  two,  is  listed. 

A later  modification  will  be  to  remove  the  second  set  of 
identical  values  from  the  print  format 
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Appendix 

This  computer  simulation  model  is  based  upon  seven  sets  of 
second-order  differential  equations.  From  the  definition 
of  the  externally-applied  forces  and  moments  Equations  1-4 
defined  below  and  the  use  of  Figure  1 showing  the  vehicle 
geometry,  the  seven  equations  of  motions  of  the  vehicles 
Equations  5-11  can  be  written. 

Externally-applied  forces  and  moments 


Q]_  ~ ~ 2fjj  (y^/V  - - kg  (y-^  - 6^) 


(1) 


2 

(Yl  " 6i)  + T-igl  + ka  ^1 


(2) 


Q6  = ' 2fL  <VV  - - kg  <*2  “ 62> 


(3) 


(4) 


lateral  motion,  fwd  wheelset: 


+ 2f L (y q/V  - + kg  yx  = kg 


(5) 


yaw  motion,  fwd  wheelset: 


(6) 


o 


body  lateral  motion: 


Mya  + Si  (2y 


a 


Yx  " Y2  + 2e0)  + ( 2y 


a “ Yl  " y2  + 2eG  =0 

(.7) 


A- 2 6 


body  yaw  motion: 


J,\  + h ~ cip  (lK  + h£  tfi  - y2> 


(8) 


K. 


+ % - kj,  '♦l  + *2)  " r-  h£  <yl  ■ y,>  ” 0 


body  roll  motion: 


h. 


10  + CQ  0 + K0  0 - -3-  e (y1  + y2  - 2yfl)-  ^ e (yx  - Y2  ~ 2y^)- 0 


lateral  motion,  aft  wheelset: 


(9) 


K, 


my  + _Z  (y2  - y + h£  ip  - e6)  + (y2  - ya  + h£  ipa  - e9) 


(10) 


+ 2fL  (y2/v  + *2)  + kg  Y2  = kg  «2 


yaw  motion,  aft  wheelset 
"2 


a£  £2 

C^2  + (^9  " V + kU,  (^2  ' ^a}  + 2fT  ( y2  + ~ i 


) 


(ID 


a£ 


- k ip  = 2frp  — 60 
a 2 T r 2 


where 


C = C h2  l2  + 2 c, 

^ y 1 

K,  = K h2  l2  + 2 k 


C0  C0  + Cy  e 
Ke  = Ke  + Ky  ^ 


The  seven  equations  of  motion  are  transformed  ^non-dimensionalized 
in  time  and  amplitude/  by  the  following  substitutions 
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82  = 
K 

-2- 

-2i- 

W2  = 
2 

0) 

3 

2 

0)  . = 


Vp^ 

a 

2f  T. 


!z, 

2m 


a>K  = 


tOrp 


V 


2f  r 


X t 


K 


- Sl 


^ ' p2 


C 

K 2 

= ye 


4mw  j 
cip 


cye 


2k, 

_Ji 

J 


2Ito. 


c i 

= _1_ 
Jco , 


5 

2 

a) , = 

2 

W0 
032 


Kyh2^2  Cyh2,2 

: , 35  = 


2 Jm5 


hi 


K, 


= s 


2 Jco 


30  “ JT^r 


K 


y 

2 

j 

g 

l 


= , 3v  = I* 


M 


2Mw, 


kCT  2 k 

= _2.  , a)  = a 


e = 


m 
--  e 


l'  = 1 , Y a = — 

K K 


Y = 


_ Yl  + Y2 


W 


Ya  = yi  - Y2 


2BK 


R" 


= OJl  = Ky  3K 

2 4 f r 

“l  l 

,2 


(jo  2 

T 


l|>  = l|»,  - 1 1<. 


= n - *2 

A 2 


(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 
(22) 

(23) 

(24) 

(25) 

(26) 
(27) 


A- 2 8 


(28) 


61  “ 61'  62  “ 62 


ip  =f(0  = 0 

a a 


(29) 


using  the  above  substitutions,  the  final  set  of  equations 
become : 


'£  = + [1  + 26l  RL  S y + 


WL 


R L + 2 I T 


R“  (y  4 e " 0)  - 28-1  RT  ^ (y  4-  e"  0) 
La  ± u ,,  a 


a)! 


- ip  = -2.  $ + 6 

,2  £. 


0) 


L 2 8 


K 


(30) 


to 


K 

2 ' A 


TA  + [i  + 26,  Rl  £]  y + R2  + (>2|  YA 


2 ^ v 

- \p-  R h l ' ip  - 2 8n  Rr  — h t! ; = 

A l a L oj  a 


^2.  / 51  " 62^ 


< l 2 8 


K 


(31) 
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2 / 2 

0)v-  ..  j.  In 

- * + I1  + 2 B2  *Sp  ^ * + I RT  - “T  i* 

U)2  WT  \ “ 

T 1 


0)-  I 


<uv  2 

= 2 6_  Rt  _K  ip  - R if,  + r = 
2 1 /.»  a la 


61  + 62 


2 e 


K 


(32) 


Rm  “k  ] $ + 


0)2 

K .. 

2 ip  + [1  + 2 3 - 

A 2 1 a 

T o)T 


41  - 4 2 


R - 
T 


ip  + r. 
VA  A 


2 3 


K 


(33) 


o)2 

_K 

2 r + 2 3,,  ~ r + r - (r  - e"  0) 
wy  a y wy  a a 


0) 

- 2 8 i (r  - e'  0)  = 0 

y wy 


(34  ) 


0), 


A *a  + 2 e*  _k  *a+  ta  - 


0). 


!i  » - 2 s4  !k  M i 

^ 'V 


h£ ' 


"2  3, 


(35) 


A- 30 


(36) 


where  the  dot 
dimensionless 


represents  differentiation  with  respect 
time  (t) . 


to 


A-31 


Wh 
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Figure  A-l.  Lateral  Dynamics  Model  for  Two  Axled  Vehicle 


APPENDIX  B 

COMPUTER  PROGRAM  TO  PREDICT  TRACK  DEFLECTION  DUE 
TO  RESPONSE  OF  RAIL  VEHICLES  TO  VERTICAL 
TRACK  IRREGULARITIES 
"HALF  CAR  MODEL" 

Prepared  by 
B.  Mackenzie 
Service  Technology  Corp . 
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PURPOSE 


The  program  computes  the  response  of  the  dynamic  model 
shown  in  Figure  1,  consisting  of  a table  of  values  of  ten 
functions.  It  then  plots  the  magnitude  of  each  function  against 
frequency . 
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Figure  B-l . The  Half  Car  Model 


EQUATIONS 

The  following  system  of  equations  was  solved  for  F^/v^, 

F2/v1'  yl/vl'  y2^Vl'  y3//vl'  y4//vl'  6]/vl'  and  <S2//vl'  aS  functions 
of  frequency.  Refer  to  Figure  1. 


where 


y.  = v,  + 6. 


y.  = v9  + 6. 


61  G11F1  + G12F2 


62  “ G11F2  + G12F1 


F1  = my4  = F2  " my3 


yl 

y 


1 + 2 j 3 — 

a)i 


2 -(^)2)+  57 

F1  + f2  - ra(y3  + y4)  + «2  V2  + 

y3  + y4’ 


8.  K,  y,  - 


+ M2y2  + M1y1  = 0 


1J1 


< 

I-1 

2t rx 

= vq  sm  — ^ — , v2  = 

V 

0 

sin  2 tt 

*•« 

yi 

2 . 0 
=-w  yi  , x = 1,  2, 

3, 

4, 

X 

= vt 

0) 

2ttV  V 0) 

X ' £ X 2 TT 

The  functions  and  G^2  are  the  dynamic  compliance  coefficients 

of  a beam  on  a visco-elastic  foundation  and  are  given  by: 
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let 


= J 

M W 


then  for  co  < co 


11 


k A fl  - ^ 

ol  u2 
a 


\~  "5~  / /i  where  A 
2 \ 3/4  o 


= 2 7T 


and 


2tt£ 


12 


G11  6 


' 2tt£  2 tt  £ \ 

sin  t + cos  r 

V A1  A1  / 


for  w > co 


where  A^  = 


2 \ 1/4 


1 - ^ 


11 


- TT  ( 1 + j ) 


/2  k XQ  ^ - i 


,3/4 


and 


12 


( - /2  3-,£ 

tt  \ e + je 


- j /2  3-,  £ 


/2  k A ( ^ - 1 

° 1 (O2 

a 


3/4 


where 


2 it  f co 

A_  ~2 

o \ co 

\ a 


\ 1/4 

1 
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INPUT 


The  program  calculations  are  in  units  of  inches,  pounds 


and 

seconds.  The  input 

variables  are 

given  below 

• 

Card 

Variable 

Units  of  Input 

Format 

1 

k 

lb/in2 

F10 .2 

2 

Ki 

lb/in 

F10 .2 

3 

k2 

lb/ in 

F10.2 

4 

Wi 

lb 

F10 .2 

5 

W 2 

lb 

F10 .2 

6 

w 

lb 

F10 .2 

7 

W 

lb/yd 

F10 .2 

8 

I 

in4 

F10 .2 

9 

E 

lb/in2 

El  0 . 7 

10 

r 

in 

F10 .2 , 

1 

in  col 

15 

lv 

in/sec 

F10 .2 , 

2 

in  col 

15 

11 

1 

in 

F10 .2 

12 

3 

- 

F10 .2 

Either  X or  V can 

be  an  input 

variable  and  the 

other  ' 

value 

is  then  calculated.  This  requires 

a 

code  number 

in 

column 

15 

of  the  data  card.  For 

X input  the 

code  is  1,  for 

V 

the  code  is 

2. 

The  values  of  w,  Wi , W2  are  read  in  units  of  pounds;  con- 
version to  mass  units  is  done  in  the  program.  The  program  also 
converts  W from  lb/yd  to  lb/in. 
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Information  about  the  desired  frequency  range  must  also  be 


read  in.  The  program  can  accept  up  to  5 consecutive  intervals 
of  frequencies  with  a different  increment  in  each  interval.  To 
do  this  the  following  cards  are  read: 

Card  Description  Format 

13  Number  of  intervals  12 

14  Increments  for  each  of 
the  intervals  (in 

order)  5F10.4 

15  Boundaries  of  the 

intervals  6F10.4 

Figure  2 shows  a sample  list  of  data  cards.  The  last  three 


cards  give  a frequency  range  that  goes  from  0.1  to  1.0  in  in- 
crements of  0.1,  from  1.0  to  2.0  in  increments  of  0.2,  from  2.0 
to  5.0  in  increments  of  0.5,  from  5.0  to  10.0  in  increments  of 
1.0,  and  from  10.0  to  100.0  in  increments  of  10.0. 

The  table  of  frequencies  can  contain  a maximum  of  100 
values . 

$DATA 

1500.0k' 

14193.00 

39856.00 

18600.00 

7000.00 

1000.00 

119.00 
312.5 

. 3 00OE+08 

1760.00  2 

82.00 
.10 
05 

.1  0.2  0.5  1.0  10.0 

.1  1.0  2.0  5.0  10.0  100.0 

Figure  2 
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OUTPUT 


A list  of  input  variables  precedes  the  tables  of  functional 
values.  The  values  printed  are  converted  values,  i.e..  Mi,  M2, 
and  m are  in  mass  units,  W is  in  lb/in.  The  printed  value  of 
each  of  the  other  variables  is  the  same  as  the  value  read  in. 

Two  computed  constants  3Q  and  k*  are  also  printed  out: 


The  values  of  the  ten  functions  (Gi 1 and  Gi2  are  printed 
and  plotted  along  with  the  eight  functions  solved  for  in 
equations  (1)  - (8))  are  printed  out  in  magnitude-phase  form. 
The  magnitudes  are  in  exponential  format;  the  phases  are  in 
floating-point  format. 

All  functions  except  G11  and  Gi2  are  normalized  by  vj  . 
This  normalization  leaves  the  output  for  y^  (i=l, 2,3,4),  61 
and  62  unitless.  Fi/vi  and  F2/vi  have  units  of  lb/in;  Gn  and 
Gi 2 are  in  in/lb. 
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OPERATING  INSTRUCTIONS 


This  program  is  written  in  FORTRAN  for  the  IBM  7094  computer. 
The  deck  is  set  up  to  run  on  the  DOT/TSC  computer. 
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SAMPLE  OUTPUT 


Figure  B-2.  Track  Deflection  Amplitude  Ratio 


.1 

I 


ACCEL ( T 1 T 2 T3  TU) 


Figure  B-3.  Car  Body  and  Truck  Accelerations  Due  to 
1-Inch  Amplitude  Track  Irregularity 
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Tl/V  T 2/ V T3/V  ru/v 


Figure  B-4 . Car  Body  and  Truck  Displacement 
Amplitude  Ratio 

1 

j 

i 

I 
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Figure  B-5.  T-Theel  Rail  Forces  Produced  by  Unit 
Track  Irregularity  (lb/in) 
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Gil  G12 


Figure  B-6 . Track  Compliance  Function 
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PROGRAM  LISTING 
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MA  I N 


EFN 


SOURCE  STATEMENT 


IFMS  ) 


05/18/71 


01  ME  NS  I CN  I RU  F { 1 024  ) 

CALL  PLGTSl IBUF,1C24) 

CALL  PLOT (0., -11., -3) 

CALL  PLCT ( J. , .5,-3  ) 

COMPLEX  G A , GB , GC  , GO 

DIMENSION  ACCV1 ( 100 ) , ACCY2( 100 ) , ACCY3( 100),ACCY4(100) 
DIMENSION  AF1  1100)  ,TF1 1100)  , AF2  1 100 ) , T F2 ( 100 ) 
DIMENSION  ADELK  ICC)  , TOELK  1001,  ADEL  2 ( 100),  TDEL2  (100) 
DIMENSION  AYIl 100 ) , TY1 ( 100)  , AY2(1C0),TY2(1C0) 
DIMENSION  AY3(10U)  , TY3(100>  , A Y4 ( 1 GO >, TY4 ( 100 ) 
DIMENSION  AG  It  100) , TGll  100) ,AG2( ICC) f TG2I1C0) 

COMPLEX  Fi  ( 100  ),  F2 ( 100  ), DELI ( 100 ), DEL2( 100 ) 

COMPLEX  Yl(lu0),Y2(100),Y3(100),Y4(100) 

COMPLEX  Bl, B2, B3,B6,B7 
COMPLEX  A1,A2,A3,A4 
COMPLEX  A 5 

COMPLEX  G 1 ( 100 ) , G2( 100 ) 

DIMENSION  OF ( 5 ) 

DIMENSION  F L ( 6 ) 

DIMENSION  F(IOO) 

COMPLEX  BL1,BL2»BL3,BL4,BL5 
G= 386. 

PI =3.14l592c5 
EPSW^:.  01 

READ  INPUT  VARIABLES 


READ 

(5,9) 

XK 

READ 

(5,9) 

XK  1 

READ 

(5,9) 

XK  2 

READ 

(5,9) 

XM  1 

R EAD 

(5,9) 

XM  2 

READ 

(5,9  ) 

XM 

READ 

(5,9) 

WT 

READ 

(5,9) 

XI 

READ 

(5,7) 

E 

READ 

(5,11) 

XX, 

R EAD 

(5,9) 

XL 

READ 

(5,9) 

BET 

» e a rj 

(5,9) 

ZETA 

7 FORMAT  (£10.7) 


9 FORMAT  (F10.2) 

11  FORMAT  ( F 10. 2 , 4X, I 1 ) 

C CONVERSIONS  ANC  COM PUT EDCON S T ANT S 
X M 1 = X M 1 /3 6c. 

XM  2=  XM  2 / 3 8 6 . 

X M = X M / 386 . 

WT=wT/36. 

WA=  SOR  T ( XK*G/wT) 

XLAMN=2 . *P  I * ( 4 ,*E*X  I/XK ) * * . 2 5 
Wi  = SORT ( XK1/XM1 ) 

XL  L = XL 

RT2  = SQRT ( 2 . ) 

C WRITE  INPUT  VARIABLES 
WRITE  (6,10) 

10  FORMAT  { 1 HO  , 1 5 F I N PUT  VARIABLES) 
WRITE  (6,12)  XK 


2 

4 

6 

8 

9 

10 

1 1 

i 2 

13 

14 

15 

16 

17 

19 

20 

21 

22 

23 

24 

25 

26 

27 


B-17 


MAIN  - EFN  SOURCE  STATEMENT  - I FN ( S ) - 05/18/71 


12  FORMAT  ( lHOt  5 X » 4 FK  = ,El5.8) 

wl R I T E (6*13)  XK1  28 

13  FORMAT  ( 1HC,4X, 5HK1  = ,E15.8) 

WRITE  (£>,14)  XK2  29 

14  FORMAT  ( 1 HO  , 4X  » 5 HK2  = ,E15.8) 

WRITE  (6,15)  XM1  30 

15  FORMAT  (1H0,4X,5FM1  = ,E15.8) 

WRITE  (6,16)  X M2  31 

16  FORMAT  ( 1 HO  » 4X  » 5HM  2 = ,E15.8) 

WRITE  (6,17)  XM  32 

17  FORMAT  ( 1H 0 » 5 X » 4H M = ,E15.8) 

WRITE  (o,18)  WT  33 

18  FORMAT  (1HU,5X,4FW  = ,E15.8) 

WRITE  (6,22)  E 34 

22  FORMAT  ( 1 HO  5 X , 4 F E = ,E15.8) 

WRITE  (6,23)  XI  35 

23  FORMAT  (1H0,5X,4HI  = ,E15.8) 

WRITE  l o , 24 ) BET  36 

24  FORMAT  ( 1HU  ,2 X ,7H6E T A = ,E15.8) 

WR  I T E ( 6 , 122)  ZETA  37 

122  FORMAT ( 1H0,2X,7HZET  A = ,El5.8) 

WRITE  (6,27)  XL  3 8 

27  FORMAT  (IH0,5X,4FL  = ,E15.8) 

IF(JA.tC.i)  GC  TQ  31 
IC(JA.EQ.2)  GO  TO  32 
IF ( (JA.NE.l ) .ANC.( JA.NE.2)  )G0  TO  6 
32  WRITE  (6,21)  XX  48 

21  FORMAT  ( 1HG, 5X  , = ,E15.3) 

GO  TG  34 

31  WRITE  (6,19)  XX  50 

19  FORMAT  ( 1H0, 9HLAMIBDA  = ,E15.8) 

GO  TC  34 

6 WRITE  (6,8)  52 

8 FORMAT  (2X,32HERR0R  ON  V/LAMBDA  IDENTIFICATION) 

34  CONTINUE 

C GENERATE  FREQUENCY  TABLE 

R E A 0 ( 5 , 20  ) NOF  53 

20  FORMAT  (12) 

READ  (5,30  IDF (II)  ,11=1,  NOF)  55 

30  FORMAT  (5F10.4) 

\F  L=NDF+1 

READ  (5,40  (FL(L)  ,L=1  ,NFL)  62 


40  FORMAT  ( 6 F 10 . 4 ) 

L = 0 
N=  0 

DU  50  11=1, NCF 
N1  = 1 
L = L + 1 
70  N = N + 1 

XF  1 =M  — 1 

XFN  = FL ( L ) +XF  i*DF III) 
XF2=XFN-FL ( L*1 ) 

XF  2 = AR  S ( XF 2 ) 
IF(XF2.GT..002)G0  TO  60 
GO  TC  49 
60  F ( N ) = X F N 
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N i = N 1 + 1 
GO  TC  7 j 

49  N = N - 1 

50  CONTINUE 
N = N + 1 

F ( N1=FL (NFL  ) 

N F = N 

CO  LCCP  TO  COMPUTE  FUNCTIONAL  VALUES 
DO  80  I = 1 » NF 
I F ( JA.EQ.l ) XL  AM=XX 
IF ( JA. EQ.2 ) XLAM  = XX/F ( I ) 

5 V=XLAM*F<  I ) 

W=2  .*P  I*F (I ) 

W W = X—  W A 
AWW  = A8  S ( W W ) 

IF  ( AWW.LT.EPSW)  GO  TO  25 
I F ( WW.GT.O.G)  GO  TC  3 
GO  TO  2 

25  WRITE  ( 6 » 26  I 

26  FORMAT (3X,5HERRGR»5X  » 4 H W = W A ) 

COMPUTE  G 1 1 AND  G12  FOR  W LESS  THAN  WA 

2 W 3 = to  / W A 

D A = ( l.-K8**2  ) 

D3=2.*ZETA*WB 
PHI =DB/OA 
PH  I =ATAN ( PH  I ) 

DC  = (CA**2+CR**2>**.125 
BET  A 1 = 2 . * PI*DC/XLAMN 
ZA=PI /( XK*XLAMN*DC**3) 

Z8=  COS  I 3 • *P  F I / 4 • ) 

Z8=ZR*ZA 

ZC=  S IN ( 3.#PHI/4. ) 

ZC= (-1 . >*ZC*ZA 
G I ( I ) =CMPLX(ZB  ,ZC ) 

GUI  ) = -G  1 ( I ) 

PHIN=<3.*PI*PHI ) / 4 . 

DE  L TA  = R T2*BFTAI*XLL 

CALL  CBLEXP(DELTA,PHIN,GA) 

PHI 3-(5.*PI+PHl ) / 4 . 

CALL  DBLEXPtDELTA, PHI3,GB) 

DE=AIMAG( CB ) 

DF=REAL(GB) 

GB=CWPLX(-DE,DF) 

PH  1 2=  P I /4 . 

DE  =CCS ( PH  I 2 ) 

DF=  S I N( PH  I 2 ) 

DF=-DF 

GC=CVPLX(CE,DF) 

G2  l I ) = G 1 ( I ) *GC  *( G4  +GB ) /RT2 
GO  TO  4 

CGMPUTE  GIL  ANC  G12  FCR  W GREATER  THAN  WA 

3 to  B = W / W A 
nA=we**2-i . 

D8=2 .*ZETA* wB 
P H I = D B / OA 

ph I = AT  AN (Phil 


- 05/18/71 


113 


115 

116 

117 

118 

122 

124 

126 

127 


132 
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133 


DC •=  ( DA**2+DB**2)**.  125 
9ETA1=2.*PI*DC/XLAMN 
GA=(-1 . ) *PI/ ( RT2*XK*XL AMN*DC**3 I 
PHI N=3.*PHl /A. 

DE=  CCS ( PH  IN  ) 

OF  = S IM  PHIN) 

G R = C MPL  X ( OE ,DF  I 
DF= -CF 

GC  -C  N P L X ( DF  » D E ) 

Gil  I ) = GA*(GB  + GC) 

GL  ( I ) = -Gl < I ) 

DELTA=RT2*BETA1*XLL 
PH  I2  = P I-PHI /A. 

P H I 3 = ( 6 .♦PI-PHI  I/A. 

CALL  DPLEXP(0ELTA,PH[2  »GC) 

CALL  DBLEXP ( DELTA, PHI  3, GDI 
DE=-A IVAG ( GC I 
OF  =RE  AL ( GD I 
GO=CMPLXIOF,DE I 
G2 ( I )=GA*GR* (GC  + GD ) 

G2(  I J =-G2  ( I I 
A 07  = 2 . *BtT*W/W 1 
D3  = l .-( W/Wl )**2 
C1={D3+G7**2)/(D3**2+D7**2) 

C2= (07*08-07) /(DR* ♦2+07**21 
Bi=CNPLX(Ci,C2) 

R2=XK2/( 2.*XK2-2.*XM2*W**2-2.*XM1*61*W**2) 

R3= (-1. ) * ( w **2 l*( XM*XM2*92+XM1*B1*B2) 

86= ( B3-XM*k**2  1/2  . 

B7=  l H 3 + X.M  * 1m  * * 2 1/2. 

Ai=Gi*R7+G2*R6 
A2  = G1*B64-G2*B7-1  . 

CS1=C0S( 2.*PI*XL/XLAM) 

CS2=SII\i(2.*PI*XL/XLAM) 

A5=CNPLX( CS1,CS2) 

Y3(  I )=(A2  — Ai*A5)/(Ai*A 1-A  2*A  2 ) 

YA (I) = ( A2*A5-A1 1/ ( A1*A1-A2*A2 ) 

Fill  ) = B 6*  Y3  ( I I + B 7*  Y A ( I ) 

F 2 ( t ) = B7*Y3 ( I I + B6  * Y A ( I ) 

Y2  ( I ) =B2 * ( Y3  ( I I ♦ Y A { I I I 
Y 1 ( I ) -B  1*  Y 2 ( I I 
DcLi  ( I ) = Y 3 ( I 1-1  . 

0EL2 ( I ) = YA ( I ) — A5 
rfRITE!  6.2C01  ) I,A5,Y3(I),YA(n»Fl(I),F2m,Yim,Y2(I),0ELl(I), 
X Dt L 2 ( I ) 

20 Cl  FORMAT! 5X  ,1 3 , ( 5 X , 2 E20 . « ) / 8X  , 

1 2(5X,2E20.8)/3X, 

2 2 (5X, 2E20.8I/8X, 

3 2 ( 5X.2E2C. 8)  /8X, 

A 2( 5X, 2E20.8I I 

BL1  = Y 3 ( I ) -DEL  1 ( I I -1  . 

8 L 2 = Y A ( I ) —DEL 2 ( I )~A5 

8L3  = F 1(  I I +F2(  I >♦(  XM*(  Y31  I > + YA(  I )I+XM2*Y2(  l H-XM1*Y1U)  l*k**2 
8LA  = Fit  I ) -F2 ( I ) -X V* ( YA I I )— Y3 ( I)  I *W **2 

8L5  = XK2*  { Y2  l I )-(  V 3 t I >*YA(  I I I /2.  >-(  XM2*Y2  ( I I *-XMl*Yl  ( I I )*W**2 
3L  = 9 F AL  ( BL  1 ) 


13A 

135 


139 

1A1 


1 A7 
1A8 


166 
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WBM  = A I MAO ( BLL) 

SL  1 1 = WBL  v*2+WBM**2 

BL1 1=S0RT ( BLi 1 > 194 

*|R  L = SE  A L ( B L 2 ) 

WlH  = £ I m A G ( RLE  ) 

9L22=WriL**2+WBM**2 

BL22=SQRT(8L22)  195 

W RL=  REAL! 3L  3 ) 

WR « = A I MAG ( RL3  ) 

8Li3=wBL*A2+WBM**2 

RL  3 3=  S ORT  ( t!L  3 3 ) 196 

WriL  = RtAL(  B L 4 ) 

W R M = a I m A G ( BL4  > 

RL  44=WBL**2+WRM**2 

BL44=SCRT  { R L44  ) 197 

W BL = k E AL ( RL  5 ) 

WRM=  A I MAG  ( L3L5  ) 

BL55=wBL**2+WRY**2 

RL 5 5= SORT ( RL 55 ) 193 

WR I Tfc ( 6 ,2  000 ) I ,RL1 1 ,8L22 » BL33,RL44, 8L55  199 

2000  FORMAT  ( 5 X , I 3 , 5 ( 5 X , E 1 7 . 3 ) ) 

80  CONTINUE 

C CONVERT  COMPLEX  NCS.  TO  MAGNITUDE-PHASE  F C R M 
DU  150  1 = 1,  IMF 

CALL  PCL AR ( PI ( I I , AF1 { I ) , TF1 ( I ) ) 209 

CALL  P0LAR(F2(n  ,AF2U)  ,TF2(I  ))  214 

CALL  POLAR ( CEL  1 ( I ) , ADEL  11  I)  , TOEL 1( I ) ) 219 

CALL  PCLARt 0EL2 II) » ACFL2 ( I ) ,TCEL2 l I ))  224 

CALL  POLAR! VI  (I  ) ,AY1(  I ) ,TY1(  I ) ) 229 

CALL  POL  A R { Y 2 ( I ) * A Y 2 ( I ) » T Y 2 ( I ) ) 234 

CALL  POLAR! Y3 ( I ) ,AY3{ I > ,TY3( M)  239 

CALL  POLAR! Y4( I ) ,AY4( I ) ,TY4( I ) > 244 

CALL  POLAR  ( Gi  ( I ) , AGi(  I) , TGK  I ))  2 49 

CALL  PGLAR(G2!I  ) ,AG2(  I ) ,TG2( I ) ) 254 

150  CJNMNUE 
C PRINTED  OUTPUT 

1 C 5 FORMAT  ( 1H0  ,F6. 2 , 3X  ,2  ( El  3 . 5 , 3X,F7  .<♦  ,5X  ) ) 

WRITE  (6,102)  257 

102  F0RVAT(lH0,2X,4hFREG,5X,9HMAC-N  F 1 / V , 5 X , 8 HP  HS  F 1/V  , 5 X , OHM  AGN  F2/V,5 
IX, 3HPHS  F 2 / V ) 

00  110  I = 1 , N F 

WRITE  (6,106)  F!  I)  ,AF1 ( I)  »TF1 ( I ),AF2(  I ),TF2(  I ) 260 

110  CONTINUE 

WRITE  (6,103)  268 

103  FORMAT  ( 1H 0 , 2 X ,4HF R E Q , 5 X , 1 1H M AGN  D r L 1 / V , 3X , 1 0 HP HS  CtLl/ V , 3X , 1 1HMAG 
IN  DEL 2/V,  3X,  10HPHS  DEL  2/ V ) 

Oj  ill  1=1, NF 

WRITE  ( 6 , 10  5)  F ( I ) , ADtLl  ( I ) ,TCELi  ( I ) , ACEL2 (I)  ,TCEL2(  I ) 271 

I i 1 continue 

WRITE  (6,107)  279 

107  FORMAT  ( 1H0,2X,4HFREQ,5X,9HMAGN  Y 1 / V , 5 X , 8H PH S Y 1 / V , 5 X , 9 HM AGN  Y 2/V, 

15X , dHPHS  Yl/V) 

DO  112  1=1 , NF 

WRITE  (6,105)  F(  I)  ,AY1(  I),TYl(I  ) *AY2(I  ) , TY2 ( I ) 282 

112  CONTINUE 

WRITE  (6,103)  290 
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108  FORMAT  ( 1 HO , 2 X , 4HFR EQ , 5 X , 9HMAGN  Y3/ V , 5X , 8HPHS  Y3 / V , 5 X , 9HM AGN  Y 4/V, 
15X.8HPHS  Y4/V) 

DO  113  1 = 1,  NF 

WRITE  (6,105)  F(  I),AY3(I)  ,TY3( I ) ,AY4(I)  ,TY4( I ) 293 

113  CONTINUE 

WRITE  (6,10.9)  301 

109  FORMAT  ( 1 HO , 2X , 4FFR EQ , 5 X , 8HMAGN  G11,6X,7HPHS  G 1 1 , 6X , 8H M AGN  G12,6X, 
17HPHS  Gil) 

DO  114  1=1, NF 

WRITE  (6,105)  F(  I),AG1( I ),TG1( I ),AG2(I),TG2( I ) 304 

114  CONTINUE 

WRITEIO^ILOJ  312 


116  FORMAT  ( 1 HO ,4HFR EQ , 6X , 10HACCEL  Yl/V, 7X, 10HACCEL  Y2 / V , 7X , 1 CH ACC E L Y 
13/V, 7X,1 OH ACCEL  Y4/V) 

DO  200  1=1, NF 
W=2.*PI*F ( I ) 

ACC Y 1 ( I ) =W**2*AY1 ( I ) /G 
ACCY2I I )= W**2*AY2 l I )/G 
AC C Y 3 ( I ) =W**2*  AY3 ( I )/G 
ACC Y4( I )=W**2*AY4l I ) /G 


WRITE(6,115)F(  I),ACCY1(  I), ACC Y 2 ( I ) , ACC Y3(  I ) ,ACCY4(  I ) 32  5 

115  FORMAT  (1H0,F6.2  ,3X  ,4( E13.6,4X) ) 

200  CONT  INUE 

SCALING  FOR  PLOT  ROUTINE 
DO  155  1=1, NF 

F(  I )=(  ALOCUIFU  ) > + l.)  *2.  338 

AF1  ( I ) = ( ALCG10(AF1 (I ) )— 1 . ) * 1 • 5 341 

AF2(I)=(AL0G1C(AF2(I) )-l.)*1.5  344 

ADELK  I )=  ( ALOGlOt  ADEL  1(  I ) ) +4. 1*1. 5 347 

A0EL2 ( I ) = ( ALCG10 1 A0EL2 ( I ) )+4.)*1.5  350 

AY1( I )=( ALQG101 AY1( I ) )+6. ) 353 

AY2(I)=(ALGC10(AY2(  I))+6.)  356 

AY3( I ) = (ALCC10(AY3(  I ) ) +6  . ) 359 

AY41  I )=(ALOG101AY4( I ) )+6. ) 362 

AG1 ( I )= ( AL0G10 ( AG1(  l ) )+7. ) *2.  365 

AG24 J ) = ( ALCG10  ( AG2 II))+7.)*2.  368 


IF ( A F 1 ( I ) .LT.O.  ) AF  1 ( I )=0. 

I F ( AF2 ( I ) .LT .0. ) AF2 1 1 ) = 0 . 

IF ( A Y 1 ( I) . LT. C. ) A Y 1 ( I ) =0. 

I F ( A Y 2 ( I ) .LT.O.  )AY2I I )=0. 

IF ( AY3 { I ) . LT.O. ) AY3 ( I )=0 . 

IF ( A Y4 ( I ) .LT.O. ) A Y4  ( I ) =0. 

IF  ( ACELK  I ) .LT.O.  ) A0EL1(  I )=0. 
IF l ADEL2 (I) .LT.O. ) ADEL2( I )=0. 
I F ( AG  1 ( I ) .LT.O. ) AG  1 ( I > = 0. 

IF ( AG2 (I) .LT .0 . ) AG2( I )=0. 

IF ( AF1 ( I ) .GT.9. ) AF1 ( I ) =9 . 
IF1AF2( I).GT.9.)AF2(I)=9. 

IF ( ACEL1 ( I).GT.9.)ACEL1(I)  = 9. 
IF ( ADEL2( I ) .GT.9. ) ADEL21I )=9. 
IF ( AG1 ( I ) . GT . 8 . > AG  1 ( I ) = 8. 

IF ( AG2 (I) .GT .8 . ) AG2  II  ) = 8 . 


ACC Y 1 ( I ) = ( ALOG 1 0 ( ACC  Y 1 ( I > )+4.  ) 435 

ACCY2I I )= ( AL0G101ACCY21 I ) ) + 4.  ) 438 

ACCY3(  I ) = < ALCG10UCCY3  ( I ) ) +4  . ) 441 

ACCY4 ( I ) =( ALUG10( ACCY4( I ) ) +4 . ) 444 
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IFCACCY1C  I ) .LT.O.  ) ACCYK  I ) = 0. 

IF ( ACCY2C I ) .LT. C. ) ACCY2C I ) = 0. 

IF ( ACCY  3 ( I ) .LT.O.)  ACCY3C  I > = 0. 

I F ( A CCY 4 ( I ) .LT.O.)  ACCY4  ( I )=G. 

155  CONTINUE 
C PLOT  ROUTINE 

CALL  LAXl S( 0. ,0. ,3 , 9 HFR EQU ENC Y , 9 , 6 .,-1, 2 ) 46  5 

CALL  LAXI  SC  C.  , 0.  , d,  19FIY1 /V  Y 2/V  Y3/V  Y4/ V , 1 9 , 8 . , -6 , 1 ) 46 7 

CALL  LL  INE ( F, AYltNF  ) 469 

CALL  LLINEl F.AY2  ,NF)  471 

CALL  LL  INEC F,AY3,NF)  473 

CALL  LL  I NEC  F,AY4»NF  ) 475 

CALL  PLOT ( 1 2 . ,0.  ,-3)  477 

CALL  LAXI SC  0. , 0. ,3 , 9HFREQUENCY, 9, 6.  ,~1 ,2 > 479 

CALL  LAX  IS  CO. »0. ,6»9HF1/V  F2/V  ,9, 9 . , +1 , 1 > 481 

CALL  LLINECF, AF1  ,NF)  433 

CALL  LL IN  E C F , AF  2 * NF  ) 485 

CALL  PLCT C12.,0.,-3)  487 

CALL  LAXI SC  0. ,0.  ,3 » 9HFREQUENCY ,9 ,6 . ,-l ,2 ) 489 

CALL  LAXISC0.*0.*6»13HDEL1/V  DEL  2 / V , 1 3 , 5 . 4 , 1 ) 491 

CALL  LLINECF, ADEL1,NFJ  493 

CALL  LLINECF, ADEL2,NF)  495 

CALL  PLOT C 12. , 0. ,-3  ) 497 

CALL  LAXI SC0.,0.,3,9HFREQUENCY,9f b.,-1,2)  499 

CALL  LAXI SCO. «C.f4,7HGll  G12  » 7 » 8.  »-7 , 1 ) 501 

CALL  LL INE ( F, AGl.NF ) 503 

CALL  LLINECF, AG2 , NF ) 505 

CALL  PLOTl  12.  , 0.  ,-3)  507 

CALL  LAXISC0.,0.,3,9HFREQUENCY,9,6. ,-1,2)  509 

CALL  LAXI SC  0. , 0. ,8 ,18HACCELCY1  Y 2 Y3  Y4 ) , 18 , 8 . ,-4, 1 ) 511 

CALL  LL  INEC  F»ACCY1,NV)  513 

CALL  LLINECF, ACCYl.NF)  515 

CALL  LLINE ( F ,ACCY2 ,NF)  517 

CALL  LL  INECF, ACCY3.NF)  519 

CALL  LLINFCF, ACCY4,NF)  521 

CALL  PLOT  C 12. ,C. »— 31  523 

CALL  PLOT  10. , 0.,  999 ) 525 

33  STOP 
END 
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INPUT  VARIABLES 


K 

= 

0.  1 50COCOOE 

CA 

K1 

= 

0.70965000E 

0 A 

K 2 

= 

C. 19928000E 

C 5 

Ml 

0.A8186528E 

02 

M2 

= 

C. 18134715E 

02 

M 

= 

0.25906736E 

01 

W 

= 

C.33055556E 

01 

E 

= 

0.30000000E 

08 

I 

= 

C. 71ACG000E 

02 

BETA 

= 

0.30000000E 

00 

ZETA 

= 

C.  CG000000E- 

-38 

L 

= 

C.820OGOOOE 

02 

V 

— 

0. 13200000 E 

05 

1 
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FREQ 

f^AGN  FI /V 

PHS  FI /V 

^AGN  F2/V 

PHS  F2/V 

0.10 

0.141584E 

02 

-3.  1398 

0. 141584E 

02 

-3.  1396 

0.2C 

C. 571735 E 

02 

-3. 1384 

0.571735E 

02 

-3.1379 

0.30 

G.130718E 

03 

-3.  1377 

0. 13C7  18E 

03 

-3.  1369 

0.40 

C. 237769E 

03 

-3.1382 

0.237768E 

03 

-3.1371 

0.50 

0.382922E 

03 

-3.  1403 

0. 38291 8E 

03 

-3. 1390 

0.60 

0.572913E 

03 

3.1384 

0.572901E 

03 

3.1399 

0.70 

0.817442E 

03 

3.  1306 

C.  81 741 OE 

03 

3.1322 

0.80 

C. 113046E 

04 

3.1181 

0 . 1 1 30  38  E 

04 

3.1199 

0 .90 

0. 153220E 

04 

3. 0990 

0. 153203E 

04 

3.1009 

1.00 

0. 205227E 

04 

3.0705 

0.205192E 

04 

3.0724 

1.20 

G.364173E 

04 

2. 9648 

0. 364C36E 

04 

2.9667 

1.40 

0. 646893E 

04 

2.7266 

0 .64642  5E 

04 

2. 7281 

1.60 

C.  102  55  5E 

05 

2. 2310 

0. 102462E 

05 

2.2320 

1.80 

0. 108 162E 

05 

1.6425 

0 . 1079  29  E 

05 

1.6426 

2.00 

C.  89959.8E 

04 

1. 3204 

0.  8 965  32 E 

04 

1.3194 

2.50 

0.6^7 06  5 E 

04 

1.1745 

0 .641422E 

04 

1. 1703 

3.00 

C.  602  180E 

04 

1.2809 

0.592023E 

04 

1.2746 

3.50 

0.657180E 

04 

1.4127 

0.640406E 

04 

1. 4067 

4.00 

0. 785126E 

04 

1.  4921 

0.759692E 

04 

1.4868 

4.50 

0.967767E 

04 

1.4899 

0 .931545E 

04 

1.4833 

5.00 

0. 117057E 

05 

1. 4047 

0.1 12184E 

05 

1.3930 

5.50 

0 . 13321  IE 

05 

1.2559 

0. 127C86E 

05 

1.2341 

6.00 

C.  1 392  63E 

05 

1. 0831 

0.132129E 

05 

1 .0453 

6.50 

0. 134557E 

05 

0.9286 

0. 126743E 

05 

C. 8670 

7. CO 

C. 123134E 

05 

0. 8152 

0.1 14837  E 

05 

0.7198 

7.50 

0.1091C8E 

05 

0.7468 

0. 10C333E 

05 

0.  6037 

8.00 

C. 947643E 

04 

0. 7201 

0.353531E 

04 

0.5092 

8.50 

0 . 8 10  4 3 4 E 

04 

0.7332 

0. 706984E 

04 

0. 4244 
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9.  GC 

C.683249E 

04 

0.7091 

0.565783E 

04 

0.3345 

9.50 

0.569284E 

04 

0.8987 

0. 43C697E 

04 

0. 2164 

10.00 

C.474584E 

04 

1.0812 

0 . 30492  3 E 

04 

0.0185 

10.50 

G.410849E 

04 

1.3553 

0.2G3465E 

04 

-0. 3873 

11-00 

C.  394026E 

04 

1 • 7031 

0. 177478E 

04 

-1.1724 

11.50 

0.43222  8E 

04 

2.  0437 

0.256331E 

04 

-1.7739 

12.00 

C.  516970E 

04 

2 . 3078 

0 .384638E 

04 

-2.0412 

14  .00 

0. 109647E 

05 

2. 7700 

0. 1 C4106E 

05 

-2.2981 

16.00 

C.  189614  = 

05 

2.9118 

0 . 186276E 

05 

-2.3178 

18  .00 

C.291 330E 

05 

2. 9798 

0.288918E 

05 

-2.2923 

20.00 

0.420899E 

05 

3.0236 

0.418950E 

05 

-2.2513 

22  .00 

C. 583329E 

05 

3.  0583 

0. 586617E 

05 

-2.2048 

24.00 

0.809629E 

05 

3.0910 

0.808012E 

05 

-2.1577 

26 .00 

0. 11120  IE 

06 

3. 1262 

0. 111036E 

06 

-2.1139 

28.00 

0. 154629E 

06 

-3.1147 

0 . 1 54448E 

06 

-2.0776 

30.00 

C. 22218  IE 

C6 

-3. 0584 

0.221959E 

06 

-2.0556 

32.00 

0.343427E 

06 

-2.9749 

0.343106E 

06 

-2.0612 

34.00 

0.642524E 

C6 

-2. 8328 

0.64 1896  E 

06 

-2.1260 

36.00 

0.345347E 

07 

-2.5457 

0 .344947E 

C 7 

-2.3438 

38.00 

C. 125329E 

07 

1.2631 

0.125192E 

07 

0.2220 

40.00 

0.173213E 

07 

2.0616 

0.173174E 

07 

-0.  5 C 03 

42.00 

C.  5 0 8 2 1 CE 

07 

-0.  6864 

0.508223E 

07 

2.3258 

44.00 

0.121489E 

07 

-0. 5060 

0. 121494E 

07 

2.2234 

46.00 

C. 760389E 

06 

-0. 4057 

0.760416E 

06 

2.2012 

48  .00 

0.583090E 

06 

-0.3409 

0. 5831C4E 

06 

2.2145 

50.00 

C. 4R8198E 

Oo 

-0.  2945 

0.488206E 

06 

2.2461 

52  .00 

C.428878E 

06 

-0.2585 

0.428983E 

06 

2.  2882 

54.00 

C.3  88480E 

06 

-0.  2290 

0.388483E 

06 

2.3367 
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56.00 

0.359041E 

Oo 

-0.2037 

0. 359043E 

06 

2.  3895 

58. 0C 

C. 336675E 

C6 

-0. 1814 

0.336676E 

06 

2.4452 

60  .00 

0. 319097E 

06 

-0.  1611 

0. 319098E 

C6 

2. 5030 

62.00 

C.3U49C2E 

06 

-0.1422 

0.304903E 

06 

2.5622 

64  .00 

0.293  1 77E 

06 

-0. 1245 

0.293177E 

06 

2. 6226 

66.00 

C.2  83298E 

06 

-0. 1077 

0 . 283299E 

06 

2.6838 

68  .00 

0.2  74830E 

C6 

-0.  0915 

0.274830E 

06 

2.  7456 

70.00 

0.267455E 

06 

-0.0758 

0.267455E 

06 

2. 8081 

72  .00 

C.260937E 

C6 

-0.0606 

0. 26C937E 

06 

2.8709 

74.00 

C.255C99E 

06 

-0.0458 

0 . 2 5 5099  E 

06 

2.9342 

76  .00 

0 . 2498C2E 

06 

-0. 0313 

0. 249802E 

06 

2.9978 

78.00 

G.244940E 

06 

-0.0172 

0 .244940E 

06 

3.0617 

80.00 

C.240425E 

06 

-0. 0033 

C.240425E 

06 

3.1259 

82.00 

0.236191E 

06 

0.0102 

0 .236191E 

06 

-3.0928 

84.00 

G. 232182E 

C6 

0. 0235 

0. 232182E 

06 

-3.0280 

36.00 

0. 228353E 

06 

0.0364 

0.2 28353E 

06 

-2.9628 

88.00 

0. 224668E 

C6 

0. 0490 

0.224668E 

06 

-2.8974 

90.00 

0. 221096 E 

06 

0.0612 

0.221096E 

06 

-2.8316 

92.00 

C. 2 17614E 

C6 

0. 0731 

0.217614E 

06 

-2.7654 

94 .00 

0.214202E 

05 

0.0846 

0.214202E 

C6 

-2.6988 

96.00 

C.  2 1 0 8 4 4 E 

C 6 

0. 0957 

0.21 0844E 

06 

-2.6318 

98  .00 

0.207526E 

06 

0. 1063 

0.207526E 

06 

-2.5643 

1 GO . 00 

C. 2G424CE 

C6 

0. 1164 

0 • 2 04240  E 

06 

-2.4964 

FREO 

M AGN  CEL1/V 

PHS  DEL  1/ V 

MAGN  0EL2/V 

PHS  QEL2/V 

C.  1C 

0.  112  564E- 

■03 

0. 0018 

0.1 12549  E 

-03 

0.0020 

0.20 

0 .4  54  5 3 3 E- 

03 

0.0032 

0.454518E 

-03 

0. 0C36 

0.30 

C.  1C3920E- 

■02 

0. 0040 

0.1 03920  E 

-02 

0.0046 

0.40 

0.189024E- 

02 

0.0036 

0. 189022F 

-02 

0. 0C44 
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0.50 

C.304419E-02 

0.0015 

0.3 04-4 16  E-02 

0.0024 

0 .60 

0.455457E-02 

-0.0030 

0. 455451E-02 

-0. 0C1 9 

0. 70 

0. 649852E-02 

-0.0108 

0 .649832  E-02 

-0.0096 

0.80 

0 .893695E— 02 

-0.0232 

0. 898649E-02 

-0. 0220 

0. 90 

0. 121 806E-01 

-0.0423 

0.121796E-01 

-0.0409 

1 .00 

0. 163149E-01 

-0.0708 

0. 163129E-01 

-0. 0694 

1 .20 

C.289497E-01 

-0.1765 

0.289419E-01 

-0.1752 

1 .40 

0. 51421 9E -01 

-0. 4148 

0. 513953E-01 

-0.4137 

1. 60 

C. 815464E-01 

-0.9104 

0.8 147 1 IE-01 

-0.9098 

1 .80 

0.85961  IE-01 

-1.  4991 

0.  8582  86E-01 

-1.4990 

2.00 

C. 714823E-01 

-1.8213 

0.713079E— 01 

-1.8220 

2.50 

0. 5137b9E-01 

-1.  9677 

0. 510561E-01 

-1.9707 

3.00 

0. 47758^ E—Ol 

-1.8616 

0 .47 180 1E-0 1 

-1.8661 

3.50 

C.520551E-01 

-1.  7297 

0. 511012E-01 

-1.7340 

4.00 

0.621283E-01 

-1.6503 

0 .606823E-01 

-1.6640 

4.50 

0. 765262E-01 

- 1.  6526 

0.  7 44664E— 01 

-1.6573 

5 .00 

0 • 92  5 06  9 E— 0 1 

-1.7385 

0 . 89  73  56  E-0 1 

-1.  7469 

5.50 

0. 1052C5E  CC 

-1. 8887 

0.101722E  00 

-1.9043 

6.00 

0 • 109896E  00 

-2.0636 

0.105839E  00 

-2.0907 

6.50 

C.106064E  00 

-2.2214 

0. 10 16 19 E 00 

-2.2654 

7.00 

0.96  8 998  E— 0 1 

-2.3392 

0. 921787 E— 01 

-2.4074 

7.50 

0.  85646  IE  — 01 

-2. 4137 

0 . 80649  IE-01 

-2.5161 

8 .00 

0.740852E-01 

-2.4488 

0.6871 86E- 01 

-2.5997 

8.50 

0. 629267E-C1 

-2. 4471 

0.57  009  3 E-0 1 

-2.6683 

9 .00 

0 . 524136E-01 

-2.4061 

0.456494E— 01 

-2.7326 

9.50 

0. 427C9CE-C1 

-2. 3149 

0 .346050  E-01 

-2.8075 

10  .00 

0.3417 10E-01 

-2.1487 

0.238692E-C1 

-2.9252 

10.50 

0. 276897E-C1 

-1. 8681 

0.138109  E~0 1 

3.0817 
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11.00 


11.50 


C.249344E-01 
C. 272639E-01 


-1.4597 
-1. 0313 


0. 78873 2E -02 
0.1 40008  E— 0 1 


2. 1183 
1.1474 


12.00 

0 . 339470E 

-01 

-0.7136 

0.247587E-01 

0. 8777 

14.00 

C.  796339E- 

-01 

-0. 2325 

0.761972E-01 

0.6903 

16  .00 

0. 140903E 

00 

-0.0994 

0. 138875E  OC 

0. 6891 

18.00 

C.217766E 

00 

-0.0342 

0.2163 08 E 00 

0.7198 

20  .00 

0. 314830E 

00 

0.0097 

0.313653E  CO 

0. 7615 

22.00 

C.439536E 

00 

0.0458 

0 .438501 E 00 

0.8070 

24.00 

0.60389  lE 

00 

0. 0799 

0.6C2911E  00 

0.  8529 

26.00 

C. 8286 1 CE 

00 

0.1159 

0.827614E  00 

0.8960 

28  .00 

0. 1 15307E 

01 

0.  1577 

0.115198E  01 

0.9329 

30. 00 

C.166354E 

01 

0.2108 

0 . 166220E  01 

0.9581 

32  .00 

0.259701E 

01 

0. 2854 

0.259509E  01 

0.9616 

34.00 

0.495579E 

01 

0.4041 

0 .4952 10E  01 

0.9201 

36  .00 

0.27382  IE 

C2 

0. 6245 

0.273592E  02 

0.7692 

38.00 

0.934146E 

01 

-2.0097 

0.933314E  01 

-2.7881 

o 

o 

o 

<r 

0. 102275E 

02 

-1.  1851 

0.1 G2245E  02 

2.7465 

42.00 

0.289589E 

02 

2.4808 

0.289599E  02 

-0.8414 

44.00 

C.704563E 

01 

2.  7139 

0.704606E  01 

-0.9965 

46.00 

0 .446593E 

01 

2.8342 

0.446613E  01 

-1.0388 

48.00 

0.  3 44  5 5 CE 

01 

2. 9068 

Q.344561E  01 

-1.0332 

50.00 

0.2B9008E 

01 

2.9554 

0.289014E  01 

- 1. 0038 

52 .00 

C. 253758E 

Cl 

2. 9905 

0.253762E  01 

-0.9609 

54.00 

G.2292o3E 

01 

3.0172 

0.229271E  01 

-C. 9C95 

56.00 

C. 2112 C4C 

Cl 

3. 0383 

0.211205E  01 

-0.8525 

58  .00 

0.1972 99 E 

01 

3.0555 

0.187300E  01 

-0. 7916 

60.00 

C.  1 8 6 2 4 8E 

Cl 

3. 0698 

0.1  86249  E 01 

-0.7278 
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62.00 

C.  1 77243E 

01 

3.0819 

0.  177243E 

01 

-0. 6619 

64 . 0 u 

0. 1 u9  7 5 5E 

01 

3.0923 

0.169755E 

01 

-0.5942 

66 . 00 

0. 1 6 3 4 2 5 E 

Cl 

3.1013 

0.  163425E 

01 

-0. 5252 

6 8 • G 0 

C.  1 57999F 

ol 

3.  1092 

0 .157999  E 

01 

-0  .4550 

70.00 

0 . 1 5329^E 

01 

3.1161 

0.  i 5 3 2 9 BE 

01 

-0.3  839 

72 . 0C 

. 149169E 

01 

3. 1223 

0.199169E 

ul 

-0.3120 

74  • 0 J 

0. 148523c 

01 

3.1277 

0. 14552  3E 

Cl 

-0. 2394 

7b  • 0 0 

0. 1 -*2  2 7 5E 

ol 

3 . 1 326 

0. 142275E 

01 

-0.1661 

74  .00 

0.13  / 3ciE 

0 1 

3.  1369 

C. 1 393clE 

01 

-0.0929 

40.  JC 

0.13o730E 

01 

3.  1407 

0 . 1 3&730E 

01 

-0.0182 

82  .00 

0. 1 34344c 

01 

-3.  1391 

C. 1 34344E 

01 

0.05o5 

84.00 

0.13216  .3  E 

01 

-3.1360 

0 .1 32168E 

01 

0.1315 

86 .00 

0 . 1 301 7oE 

01 

- 3.  1334 

C. 13C1 76E 

01 

0.2069 

88.00 

0 . 1 2 8 3 4 8 E 

01 

-3.1310 

0.128345E 

01 

0. 2 827 

90 .00 

0.  1 2 o t 5 cE 

0 1 

- 3.  1290 

0.  1 2 o 6 5 6 E 

01 

0.3587 

92.00 

0.  1 23094K 

ol 

-3.1273 

0.125094E 

01 

0. 4 350 

9 4 . 0 

C.123t4tt 

0 1 

-3.  1258 

0.  123646E 

01 

0.5116 

96  • 0 0 

0.1  222  9 i F 

ol 

-3.  1245 

0.1  22299 E 

01 

0. 5 *«4 

9H  . 00 

0.  1 2 1 04  5E 

01 

-3.1235 

0.  12  1045E 

01 

0.6&54 

1 00 .00 

C.l  1 9m75E 

ol 

-3.  1226 

0 . 1 198  75E 

01 

0.7427 

FREQ 

MAGM  Yi/V 

PH  S Yl/V 

MAGN  Y 2 / V 

PHS  Yl/V 

0.10 

0. 100412F 

01 

0.0019 

0. 100143E 

01 

0. 0C2C 

0.20 

0.  1 0 1 c 6 BE 

01 

0. 0032 

0.100577E 

01 

0.0 039 

0.30 

0.103804E 

ul 

O.0035 

0.  101322E 

01 

0. 0058 

C . 4 u 

6.  lbo922E 

Cl 

0. 0022 

0. 102408  E 

01 

0 .0077 

0.50 

0. 111153E 

0 1 

-0.0014 

0.  103885E 

01 

0. 0C94 
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0 . 6 C 

C.  1 l669oE 

01 

-0. 0083 

0.105827E 

01 

0.0109 

0.7  0 

0.  I 2381  2E 

01 

-0.0195 

0. 108336E 

01 

0. 0117 

C.  HO 

0.1 329C9E 

0 1 

-0.0  364 

0.111561E 

01 

0.0114 

0 .90 

0. 144534E 

01 

-0.0613 

0. 115713E 

01 

0. 0091 

1 . 0 u 

C.  1 594  7 lb 

01 

-0.0971 

0.121092E 

01 

0 .0031 

1 .20 

0 . 2 0 4 0 4 3 fc 

01 

-0. 2222 

0.  1373  3 IE 

01 

-0. 0334 

1 .40 

C • 2 7 7 6 9dE 

01 

-0. 4876 

0 . Id  3928E 

01 

-0.1560 

1 .60 

0. 35267CF 

01 

- 1. 0189 

0.  1 8562  BE 

01 

-0.4724 

1.80 

C.  307&3  5c 

01 

-1.6527 

0.154236E 

01 

-0.8226 

2 .00 

C. 21684 9 E 

01 

-2.  0 300 

0. 1 1 521  IE 

01 

-0. 8994 

2.  50 

C. 1 0^65  7E 

jL 

-2 . 3564 

0.891386E 

00 

-0.7901 

3.00 

0.737069c 

00 

-2.4716 

0. 512515E 

00 

-0.6634 

3.  50 

C.5  75007E 

Jo 

-2 . 5614 

0 .984575E 

00 

-0.6913 

4.00 

0.43  3 3 3 7c 

00 

- 2 • 6 6 9 5 

0. 1C7646E 

01 

-0.7823 

4 . 5 J 

0.434033E 

00 

-2.8161 

0.117266E 

01 

-0.9  301 

5 .00 

C. 389356c 

00 

- 3.  0098 

0.  1 24 55 2 E 

Cl 

-1 .1330 

5.50 

0.  ? 3 9 3 2 0 E 

00 

3 .0412 

0.125502E 

01 

-1.3775 

6.00 

0. 2 81  9 2 CE 

C J 

2.8010 

0.  11  7940E 

01 

-1.6312 

6.50 

0.223413c 

00 

2 . 5 8 7 d 

0 . 1043  19E 

01 

-1.8579 

7.00 

c»»  1 7 3 u 74E 

C 0 

2.4189 

0.  8 5 1 4 86 E 

00 

-2.0398 

7.50 

0 . 1 33876E 

00 

2 .2929 

0 .753499E 

00 

-2.1781 

8.00 

C. 104636E 

Cu 

2. 2009 

0.6  38487E 

00 

-2.2815 

8.50 

0.930261E- 

01 

2.1341 

0.545655E 

cc 

-2.3591 

9.00 

0. 669429E- 

01 

2. 0853 

0.471216E 

00 

-2.4178 

9.50 

0.548087E- 

01 

2.0494 

0.41124CE 

00 

-2.4628 

10.00 

0. 455087E- 

01 

2. 0230 

0. 362466E 

00 

-2.4977 

10.50 

0. 382682 E- 

01 

2.0037 

0.  322 3 7 IE 

CO 

-2. 5249 

B-31 


11.00 

C. 225468E-C1 

1. 9897 

0.2890546  00 

-2.5463 

11 .50 

0.2796366-01 

1.9798 

0.261C82E  CO 

-2.  5630 

12 .00 

C. 2424646-01 

1. 9730 

0.237374E  00 

-2.5761 

14.00 

0. 147767E-01 

1.9667 

0.  17127  IE  00 

-2. 6038 

16.00 

C.  <5388  106-02 

1. 9791 

0.132245E  00 

-2.6081 

18.00 

0.7099986-02 

2.0009 

0. 107536E  00 

-2. 5996 

20.00 

C. 5393566-02 

2.0279 

0.912842  E-0 1 

-2.5833 

22  .00 

0.4315696-02 

2.0582 

0. 8C5551E-C1 

-2. 5619 

24.00 

0.3617256-02 

2.0905 

0.738546E-01 

-2.5371 

26  .00 

0. 31836  16-02 

2. 1242 

0. 705654E-C1 

-2. 5098 

28.  CO 

0.2963346-02 

2.1590 

0.708535E-01 

-2.4805 

30  .00 

0.296762E-02 

2. 1944 

0. 761265E-01 

-2.4499 

32.00 

C. 333040E-02 

2.2303 

0.912286E-01 

-2.4182 

34.00 

0.4730096-02 

2. 2664 

0.1377946  00 

-2.3858 

36.00 

C. 1997636-01 

2.2989 

0 .6  16641E  00 

-2.3566 

38  .00 

0.4768286-02 

-0. 7986 

0.155468E  00 

0. 8262 

40.00 

0. 16641 9E-02 

-0.7621 

0.5714796-01 

0.8599 

42  .00 

C.  858047E-03 

-0. 7248 

0. 309531E-01 

0.8948 

44.00 

0.5110016-03 

-0.6872 

0 . 193196E— 01 

0.9302 

46  . Ou 

0. 328  5C1E-C3 

-0. 6495 

0. 1298906-01 

0.9660 

48  .00 

0.2214496-03 

-0.6115 

0.9139726-02 

1.0020 

50.00 

C.  1541516-03 

-0. 5735 

0.6629126-02 

1.0384 

52 .00 

0. 1097666-03 

-0.5354 

0 .491044E— 02 

1.0749 

54.00 

C.  7 9 4454c  — 04 

-0. 4972 

0.369152E-U2 

1.1117 

56  .00 

0.5816756-04 

-0.4589 

0.2803486-02 

1. i486 

58.00 

0.  429  1 7 46-  04 

-C. 4206 

0.2142736-02 

1.1856 
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60*00 

0.318014E-04 

-0. 3822 

0. 164276E-02 

1.2228 

62.00 

C. 23587tE-C4 

-0. 3438 

0. 125926E-02 

1 .2602 

64.00 

0.1  74  5 0 8 E - 04 

-0. 3053 

0. 561813E-03 

1.2576 

66. 00 

C.  12824  7E-04 

-0. 2668 

0.729019  E-0 3 

1.3351 

68 .00 

0.931270F-05 

-0. 2283 

0. 545482E-03 

1. 3727 

7C.  CO 

0.  to321GE-C5 

-0. 1 897 

0. 399934  E-03 

1.4104 

72.00 

0.45782  IE-05 

-0.1511 

0. 28399 3E -03 

1.4482 

74. CC 

C.  30009  3t  — 05 

-0. 1125 

0. 191339E-03 

1 .4861 

76 .00 

0.1  7888  3E-05 

-0.0739 

0.  117147E-03 

1. 5240 

78. CC 

C. 85836  3F-06 

-u. 0352 

0.57o962  E —0 4 

1.5620 

80  .00 

0. 146288 E- 06 

0.0035 

C. 10G858E-04 

1.6000 

82.00 

C. 395663F-06 

-3.0994 

0 . 2 7962  5 E-04 

-1.5035 

84  .00 

0.804593E-06 

-3.0607 

0. 532530E-04 

-1. 4654 

86.  OC 

C. 110922E-05 

-3.0220 

0.8  2 22*6  E-04 

-1 .4272 

88  .00 

0.133190E-05 

-2.  9832 

0. 101032E-03 

-1. 3890 

<30.  00 

0.  14901  IE-05 

-2.9445 

0.115608  E-03 

-1.3508 

92  .00 

0. 159762E-05 

-2. 9057 

0. 1267C9E-03 

-1.3125 

94.00 

0. 166527E-05 

-2.8669 

0 .134951E-03 

-1.2742 

96  .00 

0. 170167E-05 

-2. 8281 

0. 1 40840E -03 

-1  .2359 

98.00 

0.171 364E-05 

-2.7893 

0 . 144791E-03 

-1.1975 

100 .00 

C. 1 70662E-05 

-2. 7505 

0. 147146E-03 

-1.1591 

FREQ 

MAGN  Y 3 / V 

PHS  Y 3/ V 

MAGN  Y 4/ V 

PHS  Y4/V 

0.10 

C.lCOcllE  01 

0. 0000 

0.100011E  01 

0.0039 

0.20 

0.100045E  01 

0.0000 

0 . 100045  E 01 

0.0078 

0.30 

C.io01C4E  Cl 

0. cooo 

0.1CC1C4E  01 

0.0117 

0*40 

0.1O0189E  01 

0. 0000 

0.100189E  01 

0.0156 
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0. 5u 

0.  1 00  3 0 4t 

01 

0 . c 0 0 0 

0. 100304E 

01 

0.0195 

0 .0  0 

0. 1GG455E 

01 

-0.0000 

0.1CC455E 

01 

0.0233 

0.  70 

C.  1 006  5CE 

cl 

-0. 0001 

0. 100649E 

01 

0.0271 

0.80 

0. 100«93E 

01 

-0.0002 

0 . 10C887E 

Cl 

C. 0308 

0.90 

0.  10121 7 E 

01 

-0. 0005 

0. 1012  14E 

01 

0.0342 

1 . 0 0 

0. 101t>2  7E 

01 

-0.0011 

0. 1C1622E 

Cl 

0.0373 

1 . 20 

C.  102  85  IE 

01 

-0. 0C4  9 

0. 102825E 

01 

0.0406 

1 .40 

0 . 1 JW2  7E 

01 

-0.0198 

0. 1046  12E 

01 

0. 0325 

1 . 6 u 

C. 10515SE 

01 

-0. 0613 

0.104807E 

01 

-0  .0018 

1 .80 

0.10098  IF 

01 

-0.0850 

0.  1 003  8 IE 

01 

-0.  0154 

2 . c 0 

0.  984  71  8E 

00 

-C. 0704 

0.979267E 

00 

0 . 009  1 

2 . 5.0 

0.98128  )F 

00 

-0.0483 

0. 976665E 

oc 

0. 0516 

3 . 0 u 

C. 937  366F 

JO 

-0. 0464 

0.982042E 

00 

0.0731 

3.50 

C ,99309tF 

00 

-0.0518 

0. 9861  IbE 

00 

0. 0871 

4.00 

C. 996994E 

JO 

-0. 0622 

0.987374  E 

00 

0.0964 

4.50 

0.996672F 

00 

-0.0766 

0. 983335E 

oc 

0. 1 C24 

5. 00 

0. 98H774E 

Ou 

-0.0924 

0.971051  E 

00 

0 .1089 

5 .50 

0.972265E 

00 

-0. 1030 

C. 950964E 

00 

0. 1233 

6. 0C 

C.  952  93  6 E 

00 

-0 .1018 

0.930744E 

00 

0.1512 

6 .50 

G.9395o3E 

00 

-0. 0900 

0.91 9355E 

00 

0. 1 392 

7.  CC 

C.935255E 

u u 

-0 . U 74b 

0.91H358E 

cc 

0.2286 

7 .50 

0.937792E 

00 

-0. 0608 

C. 9241 50E 

oc 

0.2642 

8.  Cu 

C.S4-1 83E 

\J  s J 

- J . 0 5 0 1 

0 .9  3 32  16E 

00 

0. 2955 

8.50 

0.952502c 

00 

-0. 0423 

C.  5 4 3 5 9 5 E 

cc 

0. 3232 

9.00 

0 . 0 0 1 770E 

JU 

-0.0366 

0 .9544  30E 

GO 

0 • 3485 

9 .50 

0.971 5 8 1 F 

C 0 

-0. 0323 

0.  So 54  1 9 E 

CO 

0.3721 
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10.00 

0.631749E 

00 

-0.0292 

0 .976500E 

00 

0. 3945 

10.50 

C.992241E 

oc 

-0. 0267 

0. 5877C5E 

00 

0.4162 

11  .00 

0. 1 0J307E 

ol 

-0.02*7 

0 .999  10  IE 

00 

0. 4372 

11. 5o 

0.  lul42PF 

01 

-0.  0231 

0.  10 10  76E 

Cl 

0.4578 

12.00 

0. 10259 OE 

Ol 

-0.0217 

0.102276E 

Cl 

0.4780 

14.00 

0.  1 u 7 7 6 5e 

ol 

-0.  01  7C 

0. 1U7547E 

01 

0.5566 

16 .00 

u.  1 14  029F 

01 

-0.0123 

0. 1 13862E 

ol 

0. 6324 

13.00 

0. 121 7fctE 

ol 

- o . 006 1 

0. 12 1628E 

01 

0 . 7 0 5 6 

20.00 

0.1 .>  1 4 3 2 E 

01 

0.0023 

0.  1 3136  IE 

01 

0. 7761 

22 .00 

C.  1 4 3 9?  2E 

01 

0. 0140 

0. 1 4 3 8 0 9 E 

ol 

0.8430 

24.00 

0 . 1 6 02  ^3 

31 

0.0301 

0.  160  169E 

01 

C.  9 C 52 

2 6.00 

0.  i 32  5 5 7C 

01 

0. u525 

0. 1 8 24  4 2 E 

ol 

0.9610 

23  .00 

J • t 1 6 4 2 F 

01 

0.0845 

0. 2 145  12E 

Cl 

1. 0072 

30.00 

c.2c49c9E 

Ol 

J.  1 31  7 

0 .264307  E 

01 

1 .o)330 

32.00 

0. 5 5 o 7 7 C F 

01 

0 . 2064 

0.  3 565  35E 

C 1 

1. 0412 

34  . Ou 

0.5^33330 

ol 

J . 3373 

0.5  3 8 3 76  E 

ol 

0 .9875 

36  .00 

0.231994c 

C 2 

0.6037 

0.  2H17C0E 

C 2 

0. 7903 

38.00 

0. 89t>2  J 5E 

01 

-1 . 9085 

0.8952O2E 

01 

-2.3 893 

40.00 

0.  1 0 o 4 4 1 1 

02 

-1.0930 

0.  1 C64  10E 

02 

2. 6593 

4 2 • u 0 

C . 2 H 1 7 & 1 E 

02 

2.4590 

0 . 2ri  1 7 71  E 

02 

-0.8 196 

44 . 0 0 

o .6  149  72t 

0 1 

2 . 6464 

0.615015E 

01 

-0. 9290 

46.0c 

0.  362 5 79 E 

01 

2 . 7483 

0 . 352t>uCF 

01 

-0.9529 

48  .00 

0.2  43  386E 

01 

2. 3 1 30 

0.  248397F 

01 

-0.9394 

50.00 

C . 1 9 1 6 3 1 l 

01 

2 .8537 

0 . 1916  38E 

01 

-0. 9071 

52  .00 

0 . 1 55fc  2 7E 

01 

2.8937 

0. 155631E 

01 

-0.8640 

5 4 . 0 0 

0. 1 30folF 

0)1 

2.9221 

O.  1 306  34  E 

01 

-0. 8144 

56  .Jo 

0. 1122UE 

0 1 

2. 9463 

0. 1122 1 3E 

Cl 

~0. 7 605 
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58.00 

0 . 9 8048  OE 

00 

2.9676 

0 .980490E 

00 

-0.7038 

60 .00 

C.868033E 

00 

2. 9870 

0. 868040E 

00 

-0.6451 

62.00 

0. 776507E 

00 

3.0049 

0.776512E 

00 

-0. 5849 

64.00 

C. 7005C3E 

CO 

3. 0218 

0. 7C0507E 

00 

-0.5237 

66 .00 

0.636341E 

00 

3.0379 

0 .6  36343E 

00 

-0.4618 

68.00 

0. 581422E 

00 

3. 0534 

0.581424E 

00 

-0.3992 

70  .00 

0 . 533862  E 

00 

3.0684 

0 . 533863E 

oc 

-0.3361 

72.00 

C. 492256E 

CO 

3. 0830 

0.492257E 

00 

-0.2727 

74.00 

0.455538E 

00 

3.0973 

0.455538E 

00 

-0.2C89 

76.00 

0. 422883E 

cc 

3.1112 

0 . 422  884  E 

00 

-0.1448 

78.00 

0.393 64 5E 

00 

3.1249 

0. 393645E 

OG 

-0.0804 

80.00 

0. 3673G6E 

CO 

3.  1384 

0 . 367306  E 

00 

-0.0158 

82  .00 

0.343452E 

00 

-3.  1317 

0. 343452E 

CO 

0. 0451 

84. 0C 

0. 321744E 

00 

-3. 1188 

0. 321744 E 

00 

0.1143 

86.00 

0. 301903E 

00 

-3. 1062 

0. 3019C3E 

CO 

0.  1797 

88.00 

0 . 2 33  69  BE 

00 

-3.0938 

0.283698E 

00 

0 .2455 

90  .00 

0 . 26  69  3 6 E 

00 

-3.0819 

0. 266936E 

00 

0. 3115 

52.00 

C.251453E 

00 

-3.0702 

0.251452E 

00 

0.3780 

94  .00 

0.237110E 

00 

-3.0589 

0. 23711 OE 

00 

0.  4448 

56.00 

C. 223791E 

00 

-3.0481 

0.223791 E 

00 

0.5120 

98  .00 

0.211 394E 

00 

-3.0377 

0.211394E 

00 

0. 5796 

100.00 

C. 1 9983 1 E 

00 

-3.0277 

0 . 1 998  30  E 

00 

0.6477 

FREQ 

MAGN  Gil 

PH S G 1 1 

MAGN  G 12 

PH S Gil 

0.10 

0. 681840 E- 

•05 

3.1416 

0 . 1 1 3 148  E- 

•05 

-3.1416 

0.20 

0.631843E- 

05 

3. 1416 

0. 11 3149E- 

•05 

-3.1416 

0.30 

C.681849E- 

•05 

3.1416 

0.113151E- 

■05 

-3. 1416 

0.40 

0.681857E- 

05 

3.  1416 

0. 113154E- 

•05 

-3.1416 

0.50 

0.681363E- 

•05 

3.1416 

0 .113158E- 

05 

-3.  1416 

0.60 

0.681 88CE- 

05 

3. 1416 

0. 113163E- 

■05 

-3.1416 
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0.70 

0.68189  5 E-05 

3. 1A16 

0.1131 69 E-05 

-3.  1416 

0.  BO 

C.681913E—  05 

3.  1A 16 

0. 1131 75E-05 

-3.1416 

0.90 

0.681932  E-05 

3. 1 A 16 

0.1131 8 2E— 0 5 

-3.1416 

1.00 

C.  6S195AE-C5 

3. 1 A 1 6 

0. 113191E-05 

-3.1416 

1 .20 

0.68200  5E-05 

3. 1 A 16 

0.113210E-05 

-3.  1416 

1 . AO 

C.  68206  5E-05 

3. 1A16 

0. 113232E-05 

-3.1416 

1 .60 

0 .6  82 1 3 A E-05 

3 . 1 A 1 6 

0. 113258E-05 

-3.  1416 

1.80 

0.  6 8 2 21 2E- 0 5 

3. 1A16 

0. 113287E-05 

-3.1416 

2 .00 

0.682300E-05 

3 . 1 A 16 

0. 113320E-05 

-3.  1416 

2.50 

0. 682560E-05 

3. 1A16 

0.113417  E-05 

-3.1416 

3 .00 

0.6 82 878 E-05 

3 . 1 A 16 

0. 113536E-C5 

-3.  1416 

3.5C 

C.68325AE-C5 

3. 1A16 

0.113677  E-05 

-3.1416 

« 

o 

o 

0.6 83 689 E-05 

3 . 1 A 1 6 

0. 113840E-C5 

-3. 141b 

A.  50 

C.68A182E-05 

3.1416 

0.1 14Q24E-05 

-3.1416 

5 .00 

0.68A735E-05 

3. 1416 

0.  11423  IE -05 

-3.  1416 

5.  50 

C.  6853A6E-05 

3. 1416 

0 . 114461 E-05 

-3.1416 

5 .00 

0.636016E-05 

3. 1 A 1 6 

0.  11471 3E -05 

-3.  1416 

6.  50 

C.6867A9E-05 

3. 1416 

0.114987E-05 

-3.1416 

7.00 

0.6875A2E-05 

3. 1416 

0.  1 15285E-05 

-3. 1416 

7.  50 

0.688395 E-05 

3 . 1416 

0. 115606E-05 

-3.1416 

8 .00 

0.68931 0E-05 

3. 1416 

0. 115950E-05 

-3.1416 

8.50 

0.6902  87E-05 

3.1416 

0 . 116319E-05 

-3.1416 

9.00 

0.691 32  6E- 05 

3.  1416 

0. 116711E-05 

-3.1416 

9.50 

0.692A29E-05 

3.1416 

0.117  127E-05 

-3.  1416 

10  .00 

C.693597E-05 

3.  1 A 1 6 

0. 11 7569E-05 

-3.1416 

10.50 

0.69A829E-05 

3.1416 

0.1 18035E— 05 

-3.  1416 

11.00 

C.696126E-C5 

3. 1416 

0. 11 8527E-05 

-3.1416 

11.50 

0.697A90E-05 

3.1416 

0.1 19044E-0  5 

-3. 1416 
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12.00 

C.698922E-05 

3.  1416 

0. 11 9588E-05 

-3.1416 

14.00 

0 .705  339E-05 

3.1416 

0 . 122037E-05 

-3.  1416 

16.00 

0. 71291 8E- 05 

3. 1416 

0. 124949E-05 

-3.1416 

13.00 

0.721 740E-05 

3.1416 

0.128365E-05 

-3. 1416 

20.00 

0. 7319C7E-05 

3. 1416 

0. 132337E-05 

-3.1416 

22  .00 

0.743 54  IE-05 

3.1416 

0. 136929E-05 

-3.  1416 

24.00 

0.  7 5 57  5 CE- 05 

3. 1416 

0.  1422 16  E— 0 5 
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0 • o44  8 7 3 E- 0 5 

0.7854 
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2.3533 

98.00 

C.608204E-05 

0.7854 

0.402209 E-05 

2.4072 

100  .00 

0. 575436 E- 05 

0.7854 

0. 3807C7E-C5 

2.  4588 

ACCEL 

Yl/V  ACCEL  Y2/V 

ACCEL  Y3/V 

ACCEL 

FREQ 

0.10  0.102697E-02 
0.20  C . 4 1 5 90  7E -02 
0.30  0.95549  BE -02 
0.40  0.1 74969E— 01 
0.50  0.284205E-01 
0.60  C . 42°644E-01 
0.70  0.620484E-01 
0.80  C. 859974E— 01 
0.90  0.119737E  00 
1.00  0.1631C0E  00 


0. 102422E-02 
0.41 1465  E-02 
0. 932647E-02 
0.167581 E— 01 
0. 265623E-C1 
0 . 389646E— 01 
0. 542926E-C1 
0 .730238E-01 
0. 9586C3E-C1 
0.123 84 7E  00 


0. 102287E-02 
0.409289 E-02 
0.  921 43  8E-02 
0. 163950E-01 
0. 2 56468E-01 
0. 369869 E-01 
0. 504407E-01 
0.660445  E— 0 1 
0. 83851 5E-01 
0 . 103940E  00 


0. 102287E-02 
0.409289E-02 
0.921 43 8E-02 
0. 163950E-01 
0 . 2 5646  7E-0 1 
0. 369869E-01 
0. 504405  E-01 
0.66043  8E— 01 
0 • 8 38  494  E-01 
0. 103934E  00 
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APPENDIX  C 


COMPUTER  PROGRAM  TO  PREDICT  VEHICLE 
RESPONSE  TO  VERTICAL  TRACK  IRREGULARITIES 


FULL  CAR  MODEL 


EQUATIONS  OF  MOTION  FOR  VERTICAL  DYNAMICS 


For  an  evaluation  of  the  vertical  response  to  track  profile 
irregularities,  the  model  in  Figure  C-l  which  includes  damping 
in  the  secondary  suspension  can  be  used  to  describe  the  vehicle. 
Since  the  parameters  of  interest  are  the  motion  of  the  car  body 
and  displacements  of  the  wheels,  a convenient  set  of  coordinates 
to  describe  the  system  is  y2  and  02,  the  vertical  and  angular 
displacements  of  the  car  body,  and  v-^ , v2 , v^ , v^ , the  displace- 
ments of  the  trucks  at  their  connections  to  the  equalizer  springs. 
Since  the  motion  described  by  this  model  is  not  coupled  to  roll 
or  lateral  response  of  the  vehicle,  the  wheel  inputs,  v^q,  v2q, 
v30,  v^q  represent  the  average  of  the  two  rail  profiles. 


m2^2  + 2c2  ^2 
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(C-l) 


+ 2ak 
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(C-2) 
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As  discussed  previously,  the  normal  modes  containing  pitch  and 
bounch  motions  of  the  car  body  can  be  decoupled  from  the  pitch 
motion  of  the  trucks.  In  addition,  the  bounce  can  be  uncoupled 
from  unsymmetric  translation  of  the  trucks  (Figures  C-2d  and 
A-2e)  and  the  pitch  from  symmetric  translation  (Figures  C-2b  and 
A-2c)  so  that  the  dynamics  of  the  car  can  be  interpreted  in  terms 
of  two  simpler  equivalent  two  degrees  of  freedom  systems. 

Symmetric  translation  of  the  trucks  can  be  represented  by 
v^=V2=V2=v^=y^ . Substitution  in  the  equations  of  motion  reduces 
the  six  equations  to, 

raAi + c2(yi-y2)  + kiyi + k2  (a~^2)  = 0 (c-7> 


C-3 


C-4 


Figure  C-l.  Full  Car  Dynamic  Model 


0 


(C-8 ) 


m2y2  + + 2k2(y2'yl)  = 

Thus  the  vertical  motion  of  the  car  body  can  be  interpreted  in 
terms  of  the  equivalent  system  of  Figure  C-2a  which  also  is 
described  by  Equations  C-7  and  C-8.  Similarly,  the  unsymmetric 
translation  of  the  trucks  can  be  represented  by  v^=V2=y3=-V3= 

-v^  . In  this  case,  the  six  equations  of  motion  reduce  to, 

mly3  + c2(y3+a®2)  + kly3  + k2(y3+a02)  = 0 <C"9) 

I2 9 + 2a  c2(^3+a®2)  + 2a  ^2(^3+a02)  = ® (C-10) 

which  describe  the  equivalent  system  of  Figure  C-2b. 

These  equivalent  descriptions  provide  a description  of 
inputs  at  the  vehicle  wheels  in  terms  of  simple  base  motions  of 
the  equivalent  systems.  For  example,  a unit  movement  of  each  of 
the  wheel  displacements,  v^,  v^q/  v^q,  v4o'  ec3uivalent  to  a 
unit  base  displacement  y , where,  v = 1/4  (v]_o+v2 0+V3 0+V4 0 ^ ’ 
Similarly,  angular  body  motion  is  excited  by  an  equivalent  base 
motion, 

VA  = I (V10+V20-V30-V4o)  • 

Figures  C-2a  and  C-2b  indicate  that  these  base  motions  are  inputs 
to  the  system  represented  by  forcing  functions  k(y-^g)  and  kCy^g) 
on  the  right  hand  sides  of  Equations  C-7  and  C-9,  respectively. 

The  transfer  functions  for  these  inputs  are, 
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Figure  C-2.  Simplified  Models  for  Vertical  Dynamics 
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^ = U±  (s)  = 


v 
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where, 


.2  _ kl  .2  k2  m2  ^ _ a2  n C2 

co  — f co  — — ■ f y — f G — ^ 

1 m , 2 m-  m,  2'  2m-oj- 

1 2 1 p 2 2 

The  vertical  response  of  any  point  in  the  car  is  a linear  com- 
bination of  the  y 2 and  02  responses. 

A program  to  solve  the  above  equations  for  a sinusoidal 
irregularity  has  been  prepared.  The  program  solves  the  following 
equations : 


v 


1 + 


23s 

co- 

2z 


1 + 


2 3s 

CO- 

2z 


\ lz  + (2+y)u2z 


s + 


2co 
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co 
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2z  lz 


2co 


1 LU 

lz  2z 


1 + 


2G£s 
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“lz  + <2G+y)“2z 
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4 


Goj0oj,  „ 
2 lz 


given  that 

v = 1/4  (v1+v2+v3+v4) 


^2  = (x+Jl ) 

v ^ (x+L) 

v4  = v 2 (x+L) 

as  functions  of  to  with  3 as  a field  variable  and  to  plot  the 
curves  of  these  values  against  given  frequencies  on  a log-log 
axis.  The  method  by  which  the  program  solves  the  equations,  as 
well  as  the  required  input  and  sample  output  are  discussed  on 
the  following  pages. 
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INPUT 


The  input  quantities  for  the  program  are  described  below 


in  the 

order  and  units 

they  must 

be  put  into  the  data 

deck . 

Column 
tity . 

4 shows  the  equivalent  in 

Figure  1 for  the  input  quan- 

Data  Card  Quantity 

Format 

Description 

Units 

1 

V 

F7 .3 

Velocity 

MPH 

2 

XL 

F7 .3 

L in  Fig . 1 

feet 

3 

XLL 

F7 .3 

LL  in  Fig.  1 

feet 

4 

K1 

F12 .4 

Spring  Constant 
K1  in  Fig.  1 

lb  ./in . 

5 

W1 

F12 .4 

Ml  in  Fig.  1 

lb. 

6 

K2 

F12 .4 

Spring  Constant 
K2  in  Fig . 1 

lb ./in . 

7 

W2 

F12 .4 

M2  in  Fig.  1 

lb. 

8 

NDF 

12 

No  . of  F ' s 

See  dis- 
cussion 

9 

DF 

6 F10.4 

No . of  A F ' s 

of  Input 
f requen- 

10 

FL 

7 F10.4 

No.  of  F limits 

cies,  P. 

11 

Nl 

12 

No.  of  8 values 
input 

12 

B 

7 F10.4 

8 values  (field 
variable)  more 
than  7 values  may 
be  input,  but  only 
7 are  allowed  per 
data  card . 

M2a  2 

13 

G 

F12 .4 

- in  Fig.  1 

i2 

Sample  Input  Data,  found  on  page  14,  is  a list  of  the  data 


cards  in  the  order  and  format  they  should  appear  in  the  program. 


C-9 


Figure  C-3.  Full  Car  Model  Assumed  in  Program 
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INPUT  FREQUENCIES 


The  user  may  determine  the  frequency  values  against  which 
he  wants  the  equations  plotted  and  input  these  himself,  the  only 
restriction  being  that  the  plot  routine  is  equipped  to  handle  a 
range  of  0 to  100. 

The  actual  frequencies  are  computed  within  the  program  from 
the  information  input  on  the  9th  and  10  data  cards. 

The  smallest  frequency  value  computed  in  the  first  fre- 
quency limit  (FL)  and  is  punched  beginning  in  column  1 of  the 
10th  data  card.  This  quantity  is  repeatedly  incremented  by  the 
first  AF  (DF)  on  the  9th  data  card,  until  it  reaches  the  second 
F limit,  or  the  second  value  on  the  10th  data  card.  Upon 
reaching  the  second  limit,  the  program  increments  this  value  by 
the  second  AF  or  the  second  quantity  on  the  9th  data  card  until 
it  reaches  the  next  limit.  This  process  continues  until  it 
arrives  at  the  last  limit. 

As  there  will  always  be  one  more  F limit  than  AF,  there 
will  always  be  one  more  quantity  entered  on  the  10th  data  card 
than  on  the  9th. 
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The  first  equation  solved  by  the  program  is: 


v 


1 + 


20s 

0). 

2 z 


1 + 


2 Bs  , wlz 


co  n ^ + ( 2+y)  (0 


2z  2 (2+y ) B ,,3 


2z 


_ 2 2 
2“lz  w2z 


s + 


s + 


GO.  GO.. 

2z  lz 


4 

s 

2 2 

20),  co. 
lz  2z 


The  value  for  z2/v  ^ s comPute<3  for  all  co,  which  are  functions  of 
the  input  frequency  values  - co  = 2irf. 

In  the  above  equation  and  the  one  which  follows: 
s = jco  where  j = /-I 

K1 

coiz  = , referring  to  Figure  1 


k2 

J2z  M2  ' 


M, 


referring  to  Figure  1 
referring  to  Figure  1 


Let 


2 B 


2f 


C2  = “lz  + (2+u)“2z 


C3  = (2+u)6  , C4  = 2wl2  W22 


C.  - GO.  GO, 

5 2z  lz 


Then  let 
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The  denominator  can  now  be  separated  into  its  real  and 
imaginary  part.  Since  s2((joo)2)  and  s^  ( ( j oo ) 4 ) are  real  numbers 
the  real  terms  in  the  denominator  are 


DR  = 1 


A3S 


s(jaj)  and  s2  ((jco)2)  are  imaginary,  the  imaginary  terms  are 
written 


DI  = C±s  - (a2s3) 

The  denominator  may  be  expressed  as  DR  + DI  and  both  numer 
ator  and  denominator  may  be  multiplied  by  DR  - DI. 

(1  + C1s)  (DR  - DI) 

(DR  + DI)  X (DR  - DI) 

The  terms  resulting  from  this  multiplication  may  be  grouped  by 

real  and  imaginary,  so  that  the  real  terms  are  written  as 

DR  + (C,  s)  (DI) 

z2vR  = ~ ~ 

DR  + DI 


and  the  imaginary  part  is 

DI  + (C, s)  (DR) 

z2vl  = - ~ 5 

DR  + DI 


|-^-|  = ^/z2vR2  + z2vl2 

This  equation  is  solved  for  every  value  of  s for  all  $'s. 
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The  same  approach  is  used  in  solving  the  equation 

2G(3s  \ 


a 

c 

V 


1 + 


w2z  / 


1 + + 
W2z 


2 v 2 

u,  + (2G+]i )(!),  o 

2 2~~  + 

2Gwlz  W2z 


(2g2+pi|.  s3  + 
Ga)2z  W1z 


2 

2“lz 


Where 


G = 


M2a‘ 


Let 


xn  = 


x„  = 


X-  = 


x „ = 


Xr  = 


2G3 

2 

CO.  + 

lz 

2 

2 

2Gcon 

w . 

lz 

2z 

(2G2  + 

u) 

G co  ^ 00 

2z 

2 

lz 

Then 


x. 


Yi  = 


2z 


Yo  = 


x. 


X. 


X4 

Yo  = ~ 


3 
1/4 


xr 

1_ 

X. 
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The  next  step  is  separating  the  real  from  the  imaginary 

2 4 

terms  in  the  denominator,  since  s and  s are  real 

yR  = 1 - (y2s2)  + (y4s4) 
s and  s2  are  imaginary,  so 

r 3 

yl  = yxs  - y3s 


Both  numerator  and  denominator  are  multiplied  by  YR  - yl 

(1  + 2Ggs)  (yR  - yl) 

(yR  +yi)  (yR  - yl 


After  multiplication,  the  real  part  is  written  as 


APVR 


yR  - (yl  x yls) 

2 2 
yR  + yl 


and  the  imaginary  part  is 


APVI 


yl  - (yR  x s) 

2 2 
yR  + yl 


Therefore , 


r 


APVR2  + APVI2 


The  next  step  in  the  program  is  to  solve  for  z2^v0  f°r 
every  lambda  (A  = ^-)  , where  v is  velocity  in  mph  and  f is  fre- 
quency . 

Given  that 


v = 4 <VV2+  v3+  V4> 


V 

V 


1 

2 


= V 

o 


sin 


2ttx 

A 


= v^ (x+£) 
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V3  = V1 ^X+L^ 
v4  = v2 (x+L) 

Where  x - vt,  L = L in  feet  (see  Figure  l) and  1 = LL  in  feet 
(see  Figure  1) 

Let 


j {(sin  (fx)  [l  + cos  (f£)  + cos  (fL)  + cos  (fL)  cos  (f£) 
o 

- sin  (f£)  sin  (fL)  cos  (fx) [sin  (f£)  + sin  (fL) 

+ sin  (f£)  cos  (fL)  + cos  (f£)  sin  (fL)J 
If 

= 1 + cos  (f£)  + cos  (fL)  + cos  (f£)  cos  (fL) 

- sin  (f£)  sin  (fL) 

and 

= sin  (f£)  + sin  (fL)  + sin  (f£)  cos  (fL)  + cos  (f£)  sin  (fL) 


then 


v_ 

v 


2 2 
Ai  + Bi 


z„  z„ 

2 , . , . 2 v 

— may  be  written  as  — x — . so 
v 2 v > 

o vo 


z? 

(l 

r 2 2~| 

2\ 

Z 2 

A + B. 

V 

\ 4 

1 1 

— 

o 

\ 

/ 

V 

Likewise  — may  be  written  as 

v v 

o o 


V 


V 


so  in  order  to  compute 


o 
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. v a 

— for  all  small  A the  program  solves  first  for  — 

o Vo 

VA  ” I (vi+v2-v3-v4 ) 
so 

V 

— = j sin  (fx)  1 + cos  (f£)  - cos  (fL)  - cos  (fZ)  cos  (fL) 
o 

+ sin  (fZ)  sin  + cos  (fx) ^sin  (f£)  - sin  (fL) 

- sin  (f£)  cos  (fL)  - cos  (fZ)  sin  (fL) 

Let 


A^  = 1 + cos  (f Z)  - cos  (fL)  - cos  (f£)  cos  (fL) 

+ sin  ( f £ ) sin  (fL) 

and 

B2  = sin  (f£)  - sin  (fL)  - sin  (f Z)  cos  (fL)  - cos  (f£)  sin  (fL) 


then 


so 
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OUTPUT 


The  output  for  this  program  consists  of  a list  of  the  input 
variables  and  one  set  of  tables  for  each  value  of  3.  Each  of 
these  tables  contains  a list  of  all  frequencies  and  the  corre- 
sponding to  and  A,  as  well  as  the  values  of  z2/v,  a$/vQ,  z2/vq, 
and  ac^/vQ  for  that  frequency. 

In  addition,  all  curves  for  each  of  the  above  four  values 
are  plotted  by  the  CALCOMP  plotter  against  the  frequency  on  a 
log-log  axis. 

The  following  pages  contain  a listing  of  the  program, 
followed  by  sample  data  cards  and  sample  output.  The  subroutine 
which  draws  the  log  axes  (LAXIS)  is  listed  in  binary. 
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Z2/VB 


IRREGULARITY  WAVELENGTH  (FT) 

880  88  8.8  0 88 


Figure  C-4.  Car  Body  Response  Characteristics  to  Track 
Irregularity  for  Different  Damping  Ratios 
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Z2/V0 


IRREGULARITY  WAVELENGTH  (FT) 


880  88  8.8  0 88 


Figure  C-5.  Car  Body  Response  to  Sinusoidal  Track 
Irregularity 
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RP/VD 


IRREGULARITY  WAVELENGTH  (FT) 

880  88  8.8  0.88 
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11/1  5/6-) 


I BJOB 
SI BJOB 
$ I BFTC 


TIME...  0781.24  MIN 
SEXECUTE  I BJOB 

VERS  I CN  5 HAS  CONTROL. 
FIOCS 

MAIN 
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MAIN 


- EFN  SOURCE  STATEMENT  - I FN { S ) - 


c 


c 


c 


c 

c 

c 

c 


c 


c 


c 


c 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENS  ION 
DIMENSION 

dimension 

DIMENS  ION 

dimension 

DIMENSION 
DIMENSION 
DI MENS  ION 
DIMENSION 
DIMENSION 
REAL  M 


B ( 10) 

F (150  » 

PF ( 1 50 ) 

L AM( 150  ) 

DF(1U) 

FL (10) 

Z2CDV( 150, 10) 

S1150) 

APDDVl 150, 10) 

Z 2 A ( 1 50 , 10) 

APA( 150,10) 

I BUF ( 1024) 

ZADDV1150, 10) ,Z2APA( 150,10) 


REAL  LAM 

RE  AL  Kl,  K 2 

CALL  PLOTS ( I BUF ,1024 ) 

CALL  PLOT(0. ,-li. ,-3) 

CALL  PLOT (0., .5,-3) 

PI  =3.  14159265 

FIRST  DATA  CARD  IS  VELOCITY  IN  MPH 
RE  AD ( 5 , 1 ) VI 
1 FORMAT ( F7. 3 ) 


V=V 1*5280. 0/3600.0 
SECOND  DATA  CARC  IS  L IN  FEET 
R EA  D ( 5,2)  XL 

2 FORM  AT ( F 7 . 3 ) 

THIRD  DATA  CARD  IS  LL  IN  FEET 
R E A D ( 5 , 3 ) XL  L 

3 FORMAT (F7. 3) 

FOURTH  DATA  CARD  IS  Kl 

R EA  D ( 5,4)  Kl 
FIFTH  DATA  CARD  IS  W1 
R E A D ( 5,4)  W1 
SIXTH  DATA  CARD  IS  K2 
READ ( 5 ,4 ) K2 
SEVENTH  DATA  CARD  IS  W2 
READ (5,4)  W 2 

4 FORMAT { F12 .4 ) 
W1Z=SQRT(K1*386. /Wl) 
W2Z=SQRT(K2*386./W2 ) 


M=W2/Wl 

EIGHTH  DATA  CARD  IS  NUMBER  OF  DELTA  FS 
READ(5,10)  NDF 
10  FORMAT! 12 ) 

NINTH  DATA  CARD  IS  DELTA  FS 

R EAO ( 5 » 20 ) (DF{ 11 ),  1 1 = 1 , NDF ) 

20  FORMAT! 6F 10 . 4 ) 

NFL=  NDF  + 1 

TENTH  DATA  CARD  IS  F LIMITS 

READ ( 5,30)  (FL(L)»L=1,NFL) 

30  FORMAT (7F10. 4 ) 

ELEVENTH  DATA  CARD  IS  NUMBER  OF  BETAS 
READ( 5,35)  N 1 
35  FORM  AT ( 12 ) 
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2 

4 

6 


8 


9 


10 


11 

12 

13 

14 

15 

16 


17 


19 


26 


33 
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MAIN  - EFN  SOURCE  STATEMENT  - IFNIS)  - 


C TWELFTH  DATA  CARD  IS  VALUES  FOR  BETA 

R E AD ( 5 » 40 ) ( B ( J ) v Js 1 , N 1 ) 35 

40  FORMAT  I7F10.4I 

C THIRTEENTH  DATA  CARD  IS  G 

READ ( 5,41  G 42 

WR I T E ( 6 » 41 ) V 43 

41  FORMAT! 1H1, 5X,4HV  = ,F7.3) 

WRITE(6,42)  XL  44 

42  F0RMAT(1H0»5X,5HXL  = ,F7.3) 

WRITE(6,43)  XLL  45 

43  FORMAT !1H0*5X#6HXLL  = ,F7.3) 

WRI TE ( 6 » 44 ) M 46 

44  FORMAT! 1H0, 5X,4HM  = ,F10.5» 

WRITE<6,45>  W1Z  47 

45  FORMAT(1HO,5X,6HW1Z  = ,F10.5) 

WRITE(6,46I  W2Z  48 

46  FORMAT (1 HO  »5X*6HW2Z  = »F10.5) 

WR I TE ! 6 » 47 ) N1  49 

47  FORMATdHO,  5Xf  18HNUMBER  OF  BETAS  = tI2» 

C COMPUTE  FREQUENCIES 

DO  55  J=1 » N 1 

WR I T E ( 6 1 48 ) B ( J ) 52 

48  FORMAT! 1H0  *5X»7H8ETA  = ,F10.4) 


55  CONTINUE 
L=0 
N=0 

DO  50  11=1, NDF 
N2=  1 
L=L+ 1 
70  N=N  + 1 
XF1=N2-1 

XFN=FL!  L)  + XFi*DF!  ID 
XF2=XFN-FLIL+1) 

XF2=ABS(XF2I 
IF! XF2.GT. .002) GO  TO  60 
GO  TO  49 
60  F ( N) =XFN 
N2=N  2+ 1 
GO  TO  70 

49  N=N-1 

50  CONTINUE 
N=N  + 1 
L=NF  L 

F ( N ) =FL ( L ) 

NF  = N 

C COMPUTE  LAMBCAS 
DO  80  N=1 , NF 
I = N 

LAM! I )=V/F!N) 

80  CONTINUE 

C COMPUTE  Z2/VM  FOR  EVERY  S FOR  ALL  BETAS 
DO  100  J = 1 » N 1 
DO  100  1=1, NF 

C COMPUTE  S FOR  EVERY  LAMBDA 
S ( I I =2 .0*P I* F ( I ) 

C 1 = 2 ♦ 0*B ( J ) /W2  Z 
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MAIN  - EFN  SOURCE  STATEMENT  - IFN(S>  - 

C2=WlZ**2+( 2.0+M)*W2Z**2 
C3=(2.0+M>*B(J) 

C 4=  2 • 0*( W 1 Z **2 ) *( W2Z**2) 

C5  = W2  Z*W 1Z**2 
A1=C2/C4 
A2=C  3/C  5 
A3=l .0/C4 

C REAL  TERMS  IN  DENOMINATOR 

DR=1»0-(A1*S( I )**2)+A3*S( I >**4 
C IMAGINARY  TERMS  IN  DENOMINATOR 
DI=C1*S(I)-(A2*S(I)**3) 

DM=DR**2+DI**2 
C REAL  TERMS 

Z2VR= ( DR  + C 1*  S ( I )*DI  I /DM 
C IMAGINARY  TERMS 

Z2VI=(-DI+C1*S(  I ) *DR  )/ DM 
C Z2/VM  EQUALS 

Z2VM=SQRT< Z2VR**2+Z2VI**2» 

C COMPUTE  APHI/V  DELTA  FOR  EVERY  S FOR  ALL  BETAS 
X1=2.0*G*B( J) 

X2=WlZ**2+( ( 2.0*G+M ) *W2Z**2> 

X3=2 • 0*G* ( W1 Z**2 ) *( W2Z**2 ) 

X4=  ( 2»0*G**2+M)*B(J) 

X5=G*W2Z*(W1Z**2) 

Y1=X1/W2Z 
Y2=X2/X3 
Y3=X4/X5 
Y4=l • 0/X3 

C REAL  TERMS  IN  DENOMINATOR 

YR=1.0-(Y2*S ( I )**2 )+(Y4*S( I )**4> 

C IMAGINARY  TERMS  IN  DENOMINATOR 
YI=(Y1*S(  I ) )-<Y3*S(  I )*  * 3 ) 

YM=(YR**2)  + ( Y 1**2  ) 

C REAL  TERMS 

APVR= ( -YR- ( Y I*Y1*S(  III) /YM 
C IMAGINARY  TERMS 

APVI={ YI-YR*(Y1*S( Illl/YM 
C APHI/V  DELTA  EQUALS 

APVD=SQRT( APVR**2+APVI**2  > 

FB=2.0*PI/LAM( I \ 

C COMPUTE  Z2/VO  FOR  EVERY  LAMBDA  FOR  ALL  BETAS 
FBLL=FB*XLL 
FBL=FB*XL 
CFBLL  = COS( F8XLL  > 

CFB  L = COS ( F B XL  > 

SFBLL  = S I N ( F BLL ) 

SF  BL  = S I N ( F 8L  ) 

A1=1.0+CFBLL+CFBL+CFBLL*CFBL-SFBLL*SFBL 
BI=SF  BLL+SF8L+SFBLL*CFBL+CFBLL*SFBL 
C Z2/V0  EQUALS 

Z 2 VN= ( SQRT( A 1** 2+ B 1**2) *Z2VM) /4.0 
C COMPUTE  APHI/VO  FOR  EVERY. LAMBDA  FOR  ALL  BETAS 
A2=l • O+CF  BLL-CFBL-C  FBLL*CFBL+SFBLL*SFBL 
B2=SFBLL-SFBL-SFBLL*CFBL-CFBLL*SFBL 
C APHI/VO  EQUALS 

A P VN= ( SORT ( A2**2  + B2**2 ) *APVD)/4.0 
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111 


122 


124 

125 

126 
127 


128 


025 


MAIN 


EFN  SOURCE  STATEMENT 


I FN ( S ) 
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C COMPUTE  Z2/VO  ♦ A PHI/VO  129 

P1  = Z2VR*A1 
P2=Z2VI*B1 
P3=  Z 2V  I *A 1 
P4=Z2VR*B1 
ZAR 1=  P 1-P2 
ZAI 1=P3+P4 
P5=APVR*A2 
P6=APVI*B2 
P7= APVR*B2 
P8=APVI*A2 
ZAR  2 = P5-P6 
ZAI 2=P7+P8 
ZAR3=ZAR1+ZAR2 
ZA  I 3=ZAI 1 + ZA  12 

ZAVN=SQRT( ZAR3**2  + ZA  13**2 )/4.0  130 

PF(I)=( (2. 0*PI*F( I ) )**2)/386. 

Z2DDV ( I t J )=PF( I ) *Z2VN 
APODV ( I » J ) = PF ( I ) * APVN 
ZAODV( I , J) = Z AVN*  PF ( I ) 

Z2A( I » J ) = Z 2DDV l I, J )*LAM( I ) 

APA ( I , J)=APDDV( I » J ) *L  AM ( I ) 

Z2APAU  , J)  =ZADDV(  I » J ) *L  AM ( I ) 

100  CONTINUE 

DO  120  J=1,N1 

94  FORMAT (1H1, 3X»5HBETA=»F6.4) 


WRITE!6,94)  B(J)  154 

WRITE(6*110)  156 

110  FORMAT!  1 HO,  5X,4HFREQ,  14X , 5HZDDV0 , 1 1 X , 6H APDD VO , 1 1 X ,4HZ2 / A , 14X , 4HA P/ 

1A) 

DO  119  1=1, NF  i 

WRITE  (6,112)  F( I ) ,Z2DDV{ I , J) ,APDDV{ I, J) ,Z2A(I , J) ,APA(I , J)  158 

112  FORMAT (1H0,5(2X,E13«6,2X)  ) 

119  CONTINUE 

120  CONTINUE 

DO  140  J= 1 , N 1 

WR I TE ( 6 , 94 ) B ( J ) 171 

WRITE16, 125)  173 

125  FORMAT !1 HO , 5X ,4 HFR EQ , 14X , 7HZ +A DD/V , 1 IX , 7HZ 2+ AP/ A ) 

DO  130  1=1, NF 

WRITE(6,126)  F(  I ) , Z ADDV ( I , J » , Z 2APA ( I , J ) 175 

12  o F0RMAT(1H0,3(2X,E13.6,2X» ) 

130  CONTINUE 
140  CONTINUE 

DO  200  1=1 , NF 

F( I )=AL0G10(F( I ) ) 188 

DO  300  J= 1 , N 1 

Z2DDV! I , J) = A LOG  10  C Z2DDV ( I , J ) ) 19  3 

IF! APDDV! I, J) .LT. . 0000001 ) APDDV! I , J)=. 0000001 

APDDVd,  J)  = ALCGi0(  APODV!  I,J)  ) 200 

Z2A! I ,J)=AL0G10(Z2A( I ,J) ) 203 

IF (APA! I , J) .LT. .000001) APA! I , J)=. CO 0001 

APA! I , J)=ALCG10( APA! I,J  ) ) 210 

Z2APA! I , J)=AL0G10(Z2APA( I ,J) ) 213 

ZADOV! I, J )=AL0G10( ZADDV!  I, J)  ) 216 

300  CONTINUE 
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200  CONTINUE 

F(NF+1)=-1. 

F ( NF+2 ) =• 5 
00  500  J= 1 » N 1 
Z20DV(NF+1, J»=-5. 

Z2DDV(NF+2, J)=l. 

APDDV(NF+l,J>=-7. 

AP0DV!NF+2,J»=1. 

Z 2A ( NF+1 , J ) =-5 • 

Z2A(NF+2, J )=1. 

APA(NF+1, J>=-6. 

APA(NF+2,J)=1. 

Z2APAINF+1 , J ) = -4. 

Z2APA(NF+2, J I =1 . 

ZADDV1NF+1 ,J)=-5. 

ZADDV!NF  + 2, J ) = 1. 

500  CONTINUE 

CALL  LAXI SI  0. ,0. , 3 , 9HFR EQUE NC Y ,9,6. ,-l ,2)  243 

CALL  LAXIS(0.,0., 7,6HZ0D/V0,6,7. ,-5, 1>  245 

CALL  SYMBOL ( 1 • ,3 . , . 1 4 ,3 HV  =»0.,3)  247 

CALL  NUMBER(999. , 999. ,. 14, VI, 0,0. ,0)  249 

CALL  SYMB0L(999.,999.,.14,3HMPH,0.»3)  251 

CALL  SYMBOL! 1.0,2. 5,. 14, 3HG  =,0.,3)  253 

CALL  NUMBER(999.,999.,.14,G,0.,2)  255 

DO  1010  1=1, N1 

CALL  LINE(F,Z2DDV(1 ,1 ) ,NF,1  ,0)  261 

XP= (F(NF)-F( NF+1)) /F( NF+2) 

YP= ( Z2DDV! NF,  I )-Z2DDV(NF+l , I ) ) / Z2DDV ( NF  + 2 , I ) 

CALL  SYMBOL ( XP , YP  ,.07, 3HB  = ,0.  ,31  269 

CALL  NUM8ER(999.,999.,.07,B(  I),0.,2)  272 

1010  CONTINUE 

CALL  PLOT ( 12. ♦ 0. »-3)  276 

CALL  LAXIS(0.,0.,3,9HFREQUENCY,9,6.,-1,2)  278 

CALL  LAXI S(0. ,0. ,9,7HAPDD/V0,7 ,9. »-7 ,1 ) 280 

CALL  SYMB0L(1.,3.,.14,3HV  =,0.,3>  282 

CALL  NUMBER(999.,999., . 1 4 , V 1 , 0 , 0 . , 0 ) 284 

CALL  SYMBOL ( 999. ,999. ,.  14,3HMPH  ,0. ,3>  286 

CALL  SYMBOL ( 1.0, 2.5, .14, 3HG  =*0.,3)  288 

CALL  NUMBER(999.,999.,.14,G,0.,2 J 290 

DO  1020  1 = 1 , N 1 

CALL  LINE! F, APDDV!  1,  I ) ,NF, 1, 0)  296 

XP=! F(NF)-F( NF+1) )/F (NF+2) 

YP= ( APDDV! NF  » I )-APDDV( NF+1 , I ) ) /APDDVINF  + 2 , I ) 

CALL  SYMB0L!XP,YP, .07,3HB  =,0.,3I  304 

CALL  NUMBER! 999. ,999. ,. 07, Bl I ) ,0.  ,2 ) 307 

1020  CONTINUE 

CALL  PLOT ( 12 . , 0 . , -3  ) 311 

CALL  LAXI S( 0. ,0. ,3 ,9HFREQUENCY ,9 ,6. ,-l ,2 ) 313 

CALL  LAXIS(0.,0.,7,4HZ2/A,4, 7., -5,1)  315 

CALL  SYMBOLd.  ,3.  ,.14,3HV  =,0.,3>  317 

CALL  NUMBER(999. ,999. ,. 14, VI, 0,0. ,0)  319 

CALL  SYMBOL (999. , 999. , . 14, 3HMPH,0., 3 ) 321 

CALL  SYMBOL! 1.0,2. 5,. 14, 3HG  =,0.,3)  323 

CALL  NUMBER ( 999. , 999. ,. 14, G, 0. ,2 ) 325 

DO  1030  1=1, N1 
XP=(F{NF)-F(NF+1 ) ) /F(NF+2) 
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11/15/69 

MAIN  - EFN  SOURCE  STATEMENT  - IFN(S)  - 

YP= ( Z2  A ( NF » I ) -Z2A ( NF  + 1,  I ) ) / Z 2A ( N F+ 2 , I ) 

CALL  SYMB0L(XP,YP,.07,3HB  =,0.,3>  336 

CALL  NUMBER(999.,999.f.07,B<n»0.,2)  33q 

CALL  LINE(F,Z2A(i, I),NF,1,0)  342 

1030  continue 

call  PLOT ( 12 • » 0 • » -3 ) 346 

CALL  LAXIS(0.,0.,3,9HFREQUENCY,9,6.,-1,2)  348 

CALL  LAXIS(0.,0.,9,6HAPHl/A,6,9.,-6,l»  350 

CALL  SYMBOL! l.,3. , .14,3HV  =,0. ,3)  352 

CALL  NUMBER(999.,999.,.14,V1,0,0.,0)  354 

CALL  SYMBOL (999. ,999. ,. 14.3HMPH  ,0. ,3)  356 

CALL  SYMBOL ( 1.0, 2.5, .14,3HG  = ,0.,3)  358 

CALL  NUM8ER(999. ,999. ».14,G»0.»2)  360 

DO  1040  1=1, N1 
XP=(F(NF)-F(NF+l) ) / F (NF  + 2 ) 

YP= ( APA (NF ,1 )— APACNF+l, I) )/APA(NF+2,IJ 

CALL  SYMBOL(XP,YP,.07,3HB  =,0.,3>  37l 

CALL  NUMBER! 99 9. ,999., .07, B( I 1,0. ,2)  374 

CALL  LINE(F,APA(1,I) ,NF,1,0)  377 

1040  CONTINUE 

CALL  PLOT (12., 0., -3)  381 

CALL  LAXI S( 0. ,0. , 3 , 9HFR EQUENC Y , 9 ,6. ,-l ,2 ) 383 

CALL  LAXIS(0.,0.,6, 10HZ2+APDD/V0, 10,6. »-5, 1)  385 

CALL  SYMBOLd.  ,3.  ,.14,3HV  =,0.,3)  387 

CALL  NUMBER(999. ,999. ,. 14, VI, 0,0. ,0)  389 

CALL  SYMB0L(999.,999.,.14,3HMPH,0.,3)  391 

CALL  SYMBOL ( 1 . 0 ,2 . 5 , . 1 4 , 3HG  =,0.,3)  393 

CALL  NUMBER ( 999. , 999. ,.  14,G,0.  ,2)  395 

DO  1050  1=1, N1 
XP=<  F(NF )-F(NF  + l ) )/F(NF  + 2) 

Y P= ( ZAODV! NF, I ) -Z ADOV ( N F + 1 , I ) ) /ZADDV( NF  + 2 , I ) 

CALL  SYMBOL(XP,YP ,.07,3HB  =,0.,3)  406 

CALL  NUMBER(999. ,999. ,.07,B( I)  ,0.  ,2>  40g 

CALL  LINE(F,ZADDV( 1, I »,NF, 1,0)  4l2 

1050  CONTINUE 

CALL  PLOT! 12., 0., -31  416 

CALL  LAXIS(0.,0.,3,9HFR EQU ENCY ,9, 6., -I, 2)  418 

CALL  LAXI S( 0. ,0. , 7 , 7HZ2+AP/A.7 ,7 • »-4» 1 ) 420 

CALL  SYMBOL! 1.0, 3.0,. 14,3HV  =,0. ,3)  422 

CALL  NUMBER (999., 999., . 14 , VI , 0 , 0 . , 0 ) 424 

CALL  SYMR0L(999. ,999. , . 14 , 3 HM PH , 0 . , 3 ) 426 

CALL  SYMBOL! 1 . 0 , 2 . 5 , . 1 4 , 3HG  =,0. ,3)  428 

CALL  NUMBER(999.,999.,.14,G,0.,2)  430 

DO  1060  1=1, N1 

CALL  L INE! F, Z2APA! 1, I),NF,1,0>  436 

XP=(F(NF)-F(NF+1) )/F(NF+2l 

YP=( Z2APA(NF,  I )-Z2APA( NF  + 1, I ) ) / Z2 AP A ( NF  + 2 , I ) 

CALL  SYMB0L(XP,YP,.07,3HB  =,0.,3)  444 

CALL  NUMBER(999. ,999. ,.07,B( I ) , 0 . , 2 ) 447 

1060  CONTINUE 

CALL  PLGT(0. ,0. ,999)  451 

STOP 

END 
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V 


88.000 


XL  = 60.500 


XLL  = 

6.833 

M = 

4.15642 

W1Z  = 

29.31442 

W2Z  = 

8. 58324 

NUMBER 

OF  BETAS  = 7 

beta  = 

0.0000 

BETA  = 

0.0500 

BETA  = 

0. 1000 

BETA  = 

0.2000 

BETA  = 

0. 3000 

BETA  = 

0.5000 

BETA  = 

0. 7000 
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APDDVO 
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0. 641051E 

02 

0. 106629E 

01 

0. 607460E  00 
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02 
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00 

0.101108E  01 

0.187557E 

02 

0 • 6 87  947E 

00 

0.844081E  00 

0. 178057E 

02 
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0.360000E  01 
0.380000E  01 
0.400000E  01 


0.662018E  00 
0.734712E  00 
0.840986E  00 


0. 150072 E 00 
0.622237E  00 
0.981795E  00 
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ZDDVO 

APDOVO 

Z2/A 

0. 104952E 

-02 

0.214 594E-03 

0.923573E 

00 

0.445376E 

-02 

0. 156940E-02 

0. 195965E 

01 

0. 106207E 

-01 

0.451406E-02 

0.3115 39 E 

01 

0.195432E 

-01 

0. 833707E-02 

0.429950E 

01 

0.307683E 

-01 

0.11 1196E-01 

0.541523E 

01 

0.441805E 

-01 

0. 102152E-01 

0.647981E 

01 

0.612744E 

-01 

0.317879  E-02 

0.770307E 

01 

0.854611E 

-01 

0. 111060E-01 

0. 940073E 

01 

0.119900E 

00 

0 • 314509  E-01 

0.117235E 

02 

0. 164649E 

00 

0. 538986E-01 

0. 144891 E 

02 

0.216245E 

00 

0.720370E-01 

0.17  299  6E 

02 

0.27009  OE 

00 

0. 781 837E-01 

0.19  8066E 

02 

0 • 325032E 

00 

0.65  3047E-0 1 

0.220022E 

02 

0. 388387E 

00 

0. 293000E-01 

0.244129E 

02 

0.475203E 

00 

0 • 288870E-01 

0.278786E 

02 

0. 595686E 

00 

0.  101  947E  00 

0.327627E 

02 

0.741284E 

00 

0.176207E  00 

0.383723E 

02 

0. 886061E 

00 

0.233279E  00 

0.43  3185E 

02 

0.100168E 

01 

0.253341E  00 

0.463936E 

02 
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01 

0.219625E  00 

0.473636E 

02 

0 • 1 19375E 

01 

0.3201 87E-01 

0. 477501E 

02 
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01 

0.431778E  00 

0.51&622E 

02 

0.151o6lE 

01 

0.679284E  00 

0.513313E 

02 

0. 142061E 

01 

0.418093E  00 

0 . 446478  E 

02 

0. 141068E 

01 

0 . 47 1 58  IE  00 

0.413799E 

02 

0. 152072E 

01 

0. 160867E  01 

0.418199E 

02 

0. 147130E 

01 

0.205938E  01 

0 • 3 80807E 

02 

0. 133226E 

01 

0.563153E  00 

0.325663E 

02 

0 • 1 36575E 

01 

0.356840E  01 

0.316279E 

02 

0. 141619E 

01 

0. 632 592 E 01 

0.311563E 

02 
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uETA=0. 1 000 


FREQ 

ZODVO 

APOD VO 

Z2/A 

AP/A 

"•a  1 OOQOOE 

00 

0 . 104953E 

-02 

0.2146Q2E-03 

0.923586E 

00 

0. 1 88850E 

00 

0 « 20Q000E 

00 

0. 445474E 

-02 

0. 157024E-02 

0.1 96009  E 

01 

0 . 690907  E 

00 

0. 300000E 

00 

0 . 106325E 

-01 

0.452509E-0 2 

0. 311886E 

01 

0.132736E 

01 

0.4C0000E 

00 

0. 196118E 

-01 

0. 839335E-02 

0 .43 1459  E 

01 

0.184654E 

01 

0.500000E 

00 

0.310322E 

-01 

0 • 1 127  74E-0 1 

0. 546167E 

01 

0.1 98483E 

01 

0.600000E 

00 

0.449698E 

-01 

0.104729E-01 

0.659558E 

01 

0.153602E 

01 

0 » 7 GOO  00  E 

00 

0.633238E- 

-01 

0 . 330608E-02 

0.  79607  IE 

01 

0.415621E 

00 
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00 

0. 904300E 

-01 

0. 117586E-01 

0.994731 E 

01 

0.129344E 

01 
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00 

0.131411E 

00 

0.34012  IE-01 

0. 128491E 

02 

0.  332  562 E 

01 
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01 
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00 

0. 597260E-01 

0.167088E 

02 

0.525589E 

01 

0.11 0000  E 

01 

0.268085E 

00 

0.820399E-01 

0.21 446  8E 

02 

0. 65632  OE 

01 

0. 120000E 

01 

0.  370850E 

00 

0.917628E-01 

0.271957E 

02 

0.672927E 

01 

0.130000E 

01 

0.515487E 

00 
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0. 348945E 

02 
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01 
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01 
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02 
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01 

0.240000E 

01 

0.779208E 

00 

0 . 674708  E 00 
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02 

0.247393E 

02 
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BETA=0.  2000 


FREQ 
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0.177710E  00 
0 • 240826E  00 
0.  312406E  00 
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0.786577E  00 
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0.100364E  01 
0.111677E  01 
0.112574E  01 


APDDVO 

0.214599E-03 
0. 156994E-02 
0.452088E-02 
0. 036995E -02 
0 • 1 12Q56E-01 
0. 103449E-01 
0 • 323737E-02 
0. 113786E-01 
0.324165E-01 
0.558722E-01 
0.750695E-01 
0. 81 8595E-01 
0.686548E-01 
0.309088E-01 
0.305574E-01 
0 • 108068 E 00 
0.187057E  00 
0.247841E  00 
0.269193E  00 
0 • 233246E  00 
0.338986E-01 
0.454264E  00 
0. 707313E  00 
0.428376 E 00 
0.471068E  00 
0. 154146E  01 
0.183616E  01 


Z2/A 


0.923582E 

00 

0. 195994E 

01 

0. 311767 E 

01 

0.430927E 

01 

0 • 544478E 

01 

0.655181E 

01 

0.785862E 

01 

0.971831E 

01 

0. 123458E 

02 

0.156385E 

02 

0.192661E 

02 

0 • 2 29098  E 

02 

0. 265616E 

02 

0 .307708  E 

02 

0. 363700E 

02 

0.432617E 

02 

0.4948 1 3E 

02 

0 . 524365  E 

02 

0.512228E 

02 

0.47  2006  E 

02 

0. 397585E 

02 

0.384514E 

02 

0. 362966E 

02 

0.312617E 

02 

0. 294401 E 

02 

0 • 307 1 1 IE 

02 

0.291369E 

02 

AP/ A 

0.188847E  00 
0. 690774E  00 
0 « 132612E  01 
0.1841 39 E 01 
0.197219E  01 
0.151725E  01 
0.406984E  00 
0.125165E  01 
0.316961E  01 
0.491675E  01 
0.600556E  01 
0.600303E  01 
0. 464740E  01 
0 • 194284E  01 
0. 179270E  01 
0 • 594376E  01 
0.  968295E  01 
0.121167E  02 
0.  124679E  02 
0.102623E  02 
0.135594E  01 
0 . 1 66564E  02 
0.239398E  02 
0.134632E  02 
0.138180E  02 
0.423902E  02 
0.475242E  02 


0. 360000E 

01 

0. 106746E 

01 

0.445724E 

00 

0.260935E 

02 

0. 108955E 

0 . 380000E 

01 

0.114809E 

01 

0. 23271 IE 

01 

0. 265872E 

02 

0.538910E 

0.400000E 

01 

0.124810E 

01 

0 . 395401 E 

01 

0.274582E 

02 

0. 869883E 

038 


B ETA=0 • 0500 


FREQ 

ZODVO 

APDDVO 

Z2/A 

AP/A 

0.100000E 

00 

0.  104953E- 

-02 

0.214602E-03 

0.923587E 

00 

0. 188850E 

00 

uu 

o 

o 

o 

o 

o 

• 

o 

00 

0.445477E- 

-02 

0.157027E-02 

0. 196010E 

01 

0. 69092  0E 

00 

0 • 300000E 

00 

0. 106329E 

-01 

0. 452  552E-02 

0.311897E 

01 

0 .132749E 

01 

0.400000E 

00 

0 . 1 96 1 41 E 

-01 

0.339589E-02 

0.43151 IE 

01 

0. 1 8471 0E 

01 

; 0.500000E 

00 

0.310417E 

-01 

0. 112859E-01 

0.546334E 

01 

0.198631E 

01 

1 0. 6 000 00 E 

00 

0 .450000  E 

-01 

0. 104893E-01 

0. 66000  IE 

01 

0. 153843E 

01 

0. 700000E 

00 

0.  634087E- 

-01 

0.3315 82 E-02 

0.797138  E 

01 

0 .416846E 

00 

0.300000E 

00 

0.906572E 

-01 

0.118188E-01 

0. 997229E 

01 

0. 130007E 

01 

0.900000E 

00 

0.132005E 

00 

0 . 342974E -01 

0 . 129071 E 

02 

0.335353E 

01 

0.100000E 

01 

0. 191386E 

00 

0.605101E-01 

0.168419E 

02 

0. 532489E 

01 

0. 110000E 

01 

0.271 846E 

00 

0.836656E-01 

0.21 7477  E 

02 

0 .66  93  2 5 E 

01 

0 • 1 20000  E 

01 

0.380212E 

00 

0 • 944  309E- 0 1 

0.278822E 

02 

0.692493E 

01 

0. 130000E 

01 

0. 540207E 

00 

0. 824894 E-01 

0.365679E 

02 

0.558390E 

01 

i 

0.140000E 

01 

0.331280E 

00 

0.  3894  50E- 01 

0. 522519E 

02 

0.244797E 

01 

0. 150000E 

01 

0. 153380E 

01 

0 • 4067  77  E -0 1 

0.899832E 

02 

0 .238642E 

01 

0. 160000 E 

01 

0 • 357047E 

01 

0.153224E  00 

0. 196376E 

03 

0.842732E 

01 

0. 170000E 

01 

0.  3 1 81 92  E 

01 

0.284963E  00 

0.16  47  HE 

03 

0.147510E 

02 

0.180000E 

01 

0. 190906E 

01 

0.  409407E  00 

0.933319E 

02 

0 • 2001 5 5 E 

02 

0. 190000E 

01 

0. 136632E 

01 

0.486505E  00 

0.632822E 

02 

0. 225328E 

02 

0.200000E 

01 

0 • 106659E 

01 

0. 464403E  00 

0. 469301 E 

02 

0 . 204337E 

02 

0.220000E 

01 

0.788B76E 

00 

0.811715  E-0 1 

0. 315550 E 

02 

0. 324686E 

01 

0 • 240000E 

01 

0.74772  IE 

00 

0. 113701E  01 

0. 274164E 

02 

0.416903E 

02 

i 

, 0.260000E 

01 

0.719783E 

00 

0 • 142247E  01 

0.243619E 

02 

0.48145  IE 

02 

| 0.280000E 

01 

0.643216E 

00 

0.605845E  00 

0. 202154E 

02 

0. 190403E 

02 

0.300000E 

01 

0.634735E 

00 

0.465456E  00 

0.186189E 

'02 

0.136534E 

02 

0 « 320000E 

01 

0.699424E 

00 

0.108791E  01 

0. 192342E 

02 

0.299175E 

02 

0.340000E 

01 

0. 706912E 

00 

0. 946652  E 00 

0 • 1 829  66  E 

02 

0 . 24  50 16  E 

02 

C-39 


0.360000E  01 
0.380000E  01 
0.400000E  01 


0.681351E  00 
0. 756845E  00 
0 • 866068E  00 


0.173156E  00 
0 • 73  32  52  E 00 
0.1 17531 E 01 


0.166552E  S2 
0.175269E  02 
0. 190535E  02 


0.  42  32  71 E 01 
0.169806E  02 
0.258568E  02 


C-4Q 


00 

00 

00 

00 

00 

00 

00 

00 

00 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 


00 

00 

01 

01 

01 

01 

00 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

02 

02 

01 

01 

02 

02 

02 

02 

02 

02 

02 

03 

03 


ZD0V0 

APDDVO 

Z 2/A 

0. 104950E 

-02 

0.214587E-03 

0.923561E 

00 

0.445282E- 

-02 

0. 156870E— 02 

0.195924E 

01 

0.  106099E- 

-01 

0. 450657E-02 

0.311223E 

01 

0.  194848  E• 

-01 

0.830661E-02 

0.42  8666E 

01 

0.305629E• 

-01 

0.110520E-01 

0.537908E 

01 

0.436290E- 

-01 

0 • 101278E-0 1 

0. 63989  IE 

01 

0.600129E 

-01 

0.314430E-02 

0.754448E 

01 

0.828203E- 

-01 

0 • 109632E-01 

0.911023E 

01 

0.114733E 

00 

0.309938E-01 

0.112184E 

02 

0.155331E 

00 

0.530417E-01 

0. 136691E 

02 

0.200966E 

00 

0.708137E-01 

0.160773E 

02 

0.247314E 

00 

0. 767908E-01 

0.181364E 

02 

0.2936  74E 

00 

0.64101 IE-01 

0.198795E 

02 

0 • 347324E 

00 

0.  2874  76E-01 

0.218318E 

02 

0.422708E 

00 

0.283350E-01 

0.247989E 

02 

0.530875E 

00 

0. 999885E— 01 

0. 291981 E 

02 

0.668182E 

00 

0.172829E  00 

0.345882E 

02 

0.817143E 

00 

0. 228850E  00 

0.399492E 

02 

0. 957060E 

00 

0.248610E  00 

0.443270E 

02 

0.107835E 

01 

0.215622E  00 

0. 474475E 

02 

0. 134387E 

01 

0.314780E-01 

0.537550E 

02 

0.176355E 

01 

0. 425373E  00 

0.646634E 

02 

0.2002  99E 

01 

0.671283E  00 

0.677934E 

02 

0.1872  36E 

01 

0.415109E  00 

0. 588455E 

02 

0. 179603E 

01 

0.471735E  00 

0.526835E 

02 

0. 185007E 

01 

0. 163043E  01 

0.508770E 

02 

0.171160E 

01 

0.214253E  01 

0.443004E 

02 

0. 148986E 

01 

0. 62861 OE  00 

0.364187E 

02 

0.147773E 

01 

0.458233E  01 

0.342211E 

02 

0. 149195E 

01 

0. 854440E  01 

0. 328230E 

02 

C-41 


BET A=0. 0000 


FREQ 

Z+ADD/V 

Z2+AP/A 

0 • 100000E 

00 

0.11 1349E- 

-02 

0. 979872E 

00 

0.200000E 

00 

0.522330E- 

-02 

0.229825E 

01 

0.300000E 

00 

0. 132392E 

-01 

0. 388349E 

01 

0.400000E 

00 

0.245175E 

-01 

0.53938 5E 

01 

0.5  OOOOOE 

00 

0. 368821E 

-01 

0.649124E 

01 

0.600000E 

00 

0.487326E- 

-01 

0 • 714745E 

01 

0.700000E 

00 

0.637183E- 

-01 

0.  801  03 1 E 

01 

0.800000E 

00 

0.932792E- 

-01 

0 . 102607E 

02 

0 . 900000E 

00 

0. 146922E 

00 

0.143657E 

02 

0. 100000E 

01 

0 . 2243  08  E 

00 

0.197391E 

02 

0. 110000E 

01 

0. 319865E 

00 

0. 255892E 

02 

0. 120000E 

01 

0.431386E 

00 

0.3163 50E 

02 

0.130000E 

01 

0. 579988E 

00 

0. 392607E 

02 

0.1 40000  E 

01 

0.869469E 

00 

0.546523E 

02 

0.  150000E 

01 

0. 173009E 

01 

0.  10 1 498  E 

03 

0 • 160000E 

01 

0.8884 19 E 

01 

0.488630E 

03 

0. 170000E 

01 

0.417141E 

01 

0.215932E 

03 

0. 180000E 

01 

0 • 185429E 

01 

0. 906541E 

02 

0. 190000E 

01 

0. 125310E 

01 

0. 580384E 

02 

0.200000E 

01 

0. 103641F 

01 

0. 456023E 

02 

0.  220000E 

01 

0.  792  853 E 

00 

0. 317141 E 

02 

0.240000E 

01 

0.198113E 

02 

0. 726414E 

03 

0.260000E 

01 

0.282569E 

01 

0.956387E 

02 

0.280000E 

01 

0.103908E 

01 

0. 326567E 

02 

0.3  OOOOOE 

01 

0. 841565E 

00 

0 . 2468  59  E 

02 

0 • 320000E 

01 

0.144212E 

01 

0. 396584E 

02 

0. 340000E 

01 

0. 128109E 

01 

0.3315 77 E 

02 

C-42 


B ET  A=  0 . 0 500 


FREQ 

0.100000E 

00 

0.200000E 

00 

0. 300000E 

00 

0.400000E 

00 

0.5  OOOOOE 

00 

0 .6  OOOOOE 

00 

0. 700000E 

00 

0.800000E 

00 

0. 900000E 

00 

0.1 OOOOOE 

01 

0. 110000E 

01 

0 • 120000E 

01 

0. 130000E 

01 

0.140000E 

01 

0. 150000E 

01 

0.160000E 

01 

0. 170000E 

01 

0. 180000E 

01 

0.190000E 

01 

0.2  OOOOOE 

01 

0. 220000E 

01 

0.240000E 

01 

0. 260000E 

01 

0.230000E 

01 

0.3  OOOOOE 

01 

0.320000E 

01 

0.340000E 

01 

0.360000E 

01 

0. 380000E 

01 

Z+ ADD/ V 
0. 111349E-02 
0 • 5 22332E-02 
0. 132393E-01 
0.2451 74  E-01 
0.368803E-01 
0.487246E-01 
0.636910E-01 
0.931947E-01 
0.146694E  00 
0.223771E  00 
0.31 8728E  00 
0 . 428975E  00 
0.573377E  00 
0.840  5 54E  00 
0.  152773E  01 
0.348594E  01 
0.289836E  01 
0.153435E  01 
0.951007E  00 
0.696415E  00 
0. 856383E  00 
0.182903E  01 
0.213998E  01 
0.121916E  01 
0.561289E  00 
0.11 8475E  01 
0. 109555E  01 
0. 630875E  00 
0.137284E  01 


Z2+AP/A 
0. 979873E  00 
0.229826E  01 
0. 388352E  01 
0 . 539  382  E 01 
0. 649093 E 01 
0.714628E  01 
0.800686E  01 
0.102514E  02 
0.143434E  02 
0.  196918E  02 
0 • 254982  E 02 
0. 3 14582E  02 
0.388132E  02 
0. 528348E  02 
0 . 8962  68  E 02 
0. 191726E  03 
0.150033E  03 
0. 750370E  02 
0 .440466  E 02 
0. 306423E  02 
0 • 342  553  E 02 
0.670643E  02 
0 • 724  300  E 02 
0.383163E  02 
0.184645E  02 
0.325807E  02 
0.28 3553E  02 
0.15<t2l4E  02 
0.317920E  02 


C-43 


0.360000E 

01 

0.69013  3E- 

00 

0. 168699E 

02 

0.380000E 

01 

0.109971E 

01 

0. 254671 E 

02 

0.400000E 

01 

0. 155232E 

01 

0.  341 51 0E 

02 

C-44 


B ETA=0 • 1000 


FREQ 

0. 100000E 

00 

0.200000E 

00 

0. 300000E 

00 

0.400000E 

00 

0. 500000E 

00 

0 .6  00000  E 

00 

0.700000E 

00 

0.800000E 

00 

0. 900000E 

00 

0.1 00000  E 

01 

0.1 1 OOOOE 

01 

0.120000E 

01 

0.130000E 

01 

0.1 40000E 

01 

0.150000E 

01 

0 • 160000E 

01 

0. 1 70000E 

01 

0.1 80000E 

01 

0 . 1 90000E 

01 

0.2  OOOOOE 

01 

0.220000E 

01 

0.240000E 

01 

0 • 260000  E 

01 

0.230000E 

01 

0.3  OOOOOE 

01 

0. 320000E 

01 

Z+ACD/V 
0. 111349E-02 
0 • 522331E-02 
0. 132389E-01 
0.245144E-01 
0. 36  86  81 E -01 
0.486899E-01 
0. 636049E— 01 
0.929630E-01 
0.146046E  00 
0.222114E  00 
0.314867E  00 
0.420164E  00 
0.550348E  00 
0 • 766634E  00 
0.120127E  01 
0.186693E  01 
0.189264E  01 
0.1 34983E  01 
0.939715E  00 
0.734320E  00 
0.8579 79 E 00 
0.133502E  01 
0.167938E  01 
0.119237E  01 
0. 370571E  00 
0.101179E  01 


Z2+AP/A 
0.979873E  00 
0 . 229826E  01 
0. 388341 E 01 
0. 539318E  01 
0.648879E  01 
0.714118E  01 
0 • 799604 E 01 
0. 102259E  02 
0.142801E  02 
0.  19546  IE  02 
0.251893E  02 
0. 308120E  02 
0.3723 82 E 02 
0.481884E  02 
0 . 704742  E 02 
0. 102681E  03 
0.979721E  02 
0.659917E  02 
0.435236E  02 
0.323101E  02 
0.343192E  02 
0. 507840E  02 
0 • 568576E  02 
0. 374745E  02 
0.108701E  02 
0.278242E  02 


3400006 

01 

0. 102484E 

01 

0.265252E 

02 

360000E 

01 

0.604854E 

00 

0 • 1478  53  E 

02 

380000E 

01 

0. 174584E 

01 

0. 404299E 

02 

40000OE 

01 

0. 248841E 

01 

0 . 547450  E 

02 

046 


BETA=0. 2000 


FREQ 

Z+AOD/V 

Z2+AP/A 

0. 1000 JOE 

00 

0. 111349E- 

02 

0.979372E 

00 

0.2  OOOOOE 

00 

0.522318E- 

02 

0. 229820E 

01 

0. 300000E 

GO 

0. 132367E- 

■01 

0.338276E 

01 

0.400000E 

00 

0. 244998E- 

01 

0. 538996E 

01 

0. 500000E 

00 

0. 368124E- 

•01 

0.6478 9 8E 

01 

0.6  OOOOOE 

00 

0.4854C7E- 

01 

0.711931E 

01 

0. 7 OOOOOE 

00 

0.632681E- 

•01 

0.795371E 

01 

0.8  OOOOOE 

00 

0.921484E- 

01 

0. 101363E 

02 

0. 900000E 

00 

0. 143903E 

00 

0 . 140705  E 

02 

0 • 100000E 

01 

0. 216759E 

00 

0. 190748E 

02 

0. 1 1 OOOOE 

01 

0.302534E 

00 

0 . 242027  E 

02 

0.120000E 

01 

0. 392795E 

00 

0. 288050E 

02 

0. 13O0OOE 

01 

0.487986E 

00 

0.330329E 

02 

0. 140000E 

01 

0. 609492E 

00 

0. 383109E 

02 

0.1 50000  E 

01 

0.791 46 5E 

UO 

0 • 464326E 

02 

0.160000E 

01 

0. 101170E 

01 

0. E56433E 

02 

0. 170000E 

01 

0. 113969 E 

01 

0 . 5 899  59 E 

02 

0. 180000E 

01 

0. 108991F 

01 

0. 532344E 

02 

0.190000E 

01 

0. 942241 E 

00 

0.4  36406E 

02 

0. 200000E 

01 

0. 821187E 

00 

0. 361322E 

02 

0.220000E 

01 

0.927  31 IE 

00 

0. 370924E 

02 

0. 240000E 

01 

0. 132521E 

01 

0.485911E 

02 

0.2 600 DOE 

01 

0.155331E 

01 

0. 527429E 

02 

0. 280000E 

01 

0. 124034E 

01 

0. 389973E 

02 

0.300000E 

01 

0.358183E 

00 

0. 105067E 

02 

0. 320000E 

01 

0. 789646E 

00 

0.217153E 

02 

0.340000E 

01 

0.980905E 

00 

0. 253881E 

02 

C-47 


0.360J00E 

01 

0. fc0^160E 

00 

0.147684E 

02 

0.390000E 

01 

0.261304E 

01 

0. 605125E 

02 

G.4000G0E 

01 

0. 381540E 

01 

0 • 83 9389  E 

02 

C-48 


BETA=0. 3000 


FREQ 

z+aoo/v 

Z2+AP/A 

0.1  OOOOOF. 

00 

0.1113 49 E 

-02 

0.979870E 

00 

0.200000E 

00 

0.52229 2E 

-02 

0. 229808E 

01 

0.300000E 

00 

0.132326E 

-01 

0.388155E 

01 

0 • 400000E 

00 

0. 2447A3E 

0 

1 

0. 538435E 

01 

0. 5 00 00 OE 

00 

0 • 3 67 1 95  E 

-01 

0.546263E 

01 

0.6  000  JOE 

00 

0.433015E 

-01 

0.  7 0 842  2 E 

01 

0. 700000E 

00 

0.627  543  E 

-ol 

0.788912E 

01 
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APPENDIX  D 


PRELIMINARY  ESTIMATES  OF  THE  ACCURACIES  REQUIRED 


FOR  SIMULATION  OF  LATERAL  DYNAMICS  AND  TRACTION 


SUMMARY 


In  order  to  accurately  simulate  the  behavior  of  a rail 
vehicle  operating  on  a track,  it  is  necessary  for  the  wheel/rail 
simulator  modules  to  maintain  an  accurate  geometric  interrela- 
tionship between  the  roller  motions. 

To  provide  a good  simulation  of  lateral  dynamics,  the 
calculations  given  below  indicate  that  it  is  desirable  to  control 
the  actuator  motions  to  the  following  accuracies: 

Roller  - 0.1  milliradian  (20  arc  seconds) 

Lateral  Vibration  - 0.01%  of  forward  velocity 

Longitudinal  Vibration  - 0.01%  of  forward  velocity 

(0.025  inch/second  at  15  mph) 

A deterioration  of  performance  by  a factor  of  ten  would  tend  to 
invalidate  the  simulation. 

Stray  longitudinal  vibrations  may  be  theoretically  compen- 
sated for  _by  changes  in  roller  speed  to  maintain  the  relation: 

V = r A = x 
a a a 

to  within  0.01%  of  the  simulated  forward  velocity  (see  page  D-3 
for  definitions) . 

The  evaluation  testing  of  the  prototype  module  should  there- 
fore include  accurate  measurements  of  the  stray  lateral  and 
longitudinal  vibrations  to  an  accuracy  of  0.025  inch  per  second 
under  the  full  range  of  loads,  amplitudes  and  frequencies  con- 
templated for  the  machine.  The  measurements  in  the  tests  should 
also  include  measurements  of  the  controlled  yaw  angle  of  the 
roller  to  an  accuracy  of  20  arc  seconds. 
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Evaluation  testing  should  also  demonstrate  decrowning  of 
less  than  0.5  inch. 

DISCUSSION 

The  actuator  modules  of  the  wheel  rail  simulator  are  in- 
tended to  provide  motions  and  forces  on  a wheel  of  a rail  vehicle 
that  are  identical  to  those  that  would  be  experienced  by  a rail 
vehicle  operating  on  a real  track.  As  shown  in  the  following 
paragraphs,  a simulation  of  the  forces  requires  that  the  roller 
a gular  velocity,  horizontal  motions  and  yaw  angle  satisfy  the 
following  relations 


V = Simulated  forward  velocity 

= Roller  angular  velocity 

r = Roller  radius 
a 

ip^  = Actuator  yaw  angle 

ip  = Wheelset  yaw  angle 

6^  = Actuator  lateral  position 

= Actuator  longitudinal  position 

x = Wheelset  longitudinal  displacement. 

If  these  relations  are  not  satisfied,  forces  will  be  exerted  on 
the  wheel  which  will  not  exist  in  practice.  Equation  1 requires 
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that  the  apparent  forward  velocity  of  the  roller  be  equal  to  the 
apparent  velocity  of  the  rail  as  seen  by  an  observer  travelling 
on  a vehicle  having  velocity  V. 

Equation  2 requires  that  the  instantaneous  lateral  velocity 
of  the  roller  be  zero  since  the  rail  will  not  have  an  instantan- 
eous lateral  velocity  except  for  the  effects  of  rail  compliance. 

Equation  3 requires  that  the  roller  center  must  remain 
directly  below  the  axle  center  to  avoid  destabilizing  forces 
and  moments. 

The  error  in  satisfying  Equation  1 will  result  in  a longi- 
tudinal creep  force: 


For  most  vehicles  of  interest  the  creep  coefficient,  f,  will  be 
about  150  times  the  normal  force. 


A stray  longitudinal  vibration  of  0.1%  of  the  forward  velocity 
would  result  in  a creep  force  of  15%  of  the  normal  force.  Since 
the  available  adhesion  is  typically  15%  to  30%  of  the  normal 
force,  this  could  result  in  a total  loss  of  adhesion  or  a re- 
duction in  the  apparent  braking  or  acceleration  capability  of 
the  vehicle  of  50  to  100%.  It  is  therefore  desirable  that  the 
velocity  associated  with  stray  longitudinal  vibrations  be  limited 
to  less  than  0.01%  of  the  forward  velocity.  At  a forward 
velocity  of  15  mph  this  represents  a velocity  of  about  0.025 
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inch  per  second.  At  1 Hz  this  would  limit  stray  longitudinal 

-4 

vibrations  to  an  acceleration  of  4.05  x 10  g.  At  10  Hz  this 

-3 

requirement  would  be  relaxed  to  4.05  x 10  g. 

Similarly  an  error  in  satisfying  Equation  2 would  result  in 
a lateral  creep  force  of 


6,  r A 

F = 150  N — - ■ a \h 
c L V V ^a 

-3 

Therefore  an  error  m ip  of  10  radian  would  result  in  a creep 
force  of  15%  of  the  normal  force  resulting  in  a loss  of  lateral 
adhesion  of  50  to  100%.  Therefore,  it  is  desired  that  the  yaw 
angle  be  controlled  to  better  than  0.1  milliradian  at  all  fre- 
quencies. Similarly  the  stray  lateral  vibrations  should  be 
controlled  to  better  than  0.01%  of  the  simulated  forward  veloc- 
ity. At  15  mph  this  implied  0.025  inch  per  second. 


vaT  = (V  ~ (ro  + “(y-«i)) 

vbT  = (v  + s.i)  - (ro  - a(y-q)) 


VaL  = ^ 

V 

r 

o 

(rQ  + ^(y-Sj) 

^ ~ y 

- Vip 

vbL = y 

V 

r 

o 

(ro  - a(y-S2)) 

\p  ~ y 

- V 

Referring  to  Figures  3 and  4 the  corresponding  creep 
velocities  on  the  simulator  are: 
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v T = r A i|j  - 6.  + 
aL  a ara  1 

(y 

- Vi|>) 

vUT  = r.  A,  ip,  - S + 
bL  b b b 

(y 

- Vijj) 

For  the  creep  velocities  to  be  equal  in  both  cases  we  require: 

e.  = V - ( r A - x i = 0 
1 \ a a a. ) 

e2  = (raVa  ‘ Sl)  “ 0 

el  = V - (rbAb  - *b)  = ° 

e2  = (rbVb  - h)  = 0 

Creep  coefficients  used  by  British  Rail  (PB  192718)  for 
analysis  of  stability  of  the  LIM  test  vehicle  range  from  100  to 
200  times  the  normal  force.  If  the  above  requirements  are  not 
satisfied  there  will  be  undesired  forces  on  wheel  "a"  of 
(Figure  1)  . 

A - x 1 
l a a 

V 

,-hl 

V 


Ft  = f ^ = 150  N 
e2 

Ft  = f .-==■  = 15  0 N 
1j  V 


The  existence  of  these  forces  will  reduce  the  observed  adhesion 
and  because  of  the  non-linearity  of  the  creep  phenomenon  produce 
a change  in  the  apparent  creep  coefficient. 

In  addition  to  the  above  creep  force  errors,  the  curvature 
of  the  roller  produces  a destabilizing  force  on  the  wheel  as 
shown  in  Figure  5.  The  destabilizing  forces  on  the  two  rollers 
corresponding  to  a wheelset  are: 


D-6 


x - Zijj  - x 


a 


F 


aTN 


N 


F. 


x + Hip  - x^ 


N 


bTN 


for  a 40  inch  diameter  roller  and  30  inch  diameter  wheel  an 
error  of  0.35  inch  would  produce  a force  of  1%  of  the  normal 
force.  An  error  of  3.5  inches  would  produce  a force  of  10%  of 
the  normal  force.  A positioning  accuracy  of  0.51  inches  would 
produce  a force  error  of  1.5%  which  would  represent  5%  of  the 
available  adhesion  with  an  adhesion  coefficient  of  0.3. 
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Figure  D-l.  Wheelset  Travelling  on  Nominally  Straight  Track 
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Figure  D-2.  Wheelset  Travelling  on  Nominally  Straight  Track 


roller 


x 
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Figure  D-3 . Wheelset  Rolling  on  Track  Simulated  by  Roller 


Figure  D-4 . Wheelset  rolling  on  Track  Simulated  by  Rollers 
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Figure  D-5 . Force  Produced  by  "Decrowning" 


D-12 


! 


-n  T1 

0 o 

1 i 

5 o 

C Q 


n 

0 

1 


□ 034 


